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The paper is a report on a tabulation of units performed on the EDVAC at the U. 8. 
Army Proving Grounds in Aberdeen, Maryland. The algebraic number fields involved 
were 45 cyclie cubic fields of discriminant /?, where / is one of the primes of form 3m-+1 
from 7 through 499. The object of the search was the discovery of units through a specific 
method, an algorithm based on Minkowski’s geometric number theory but particularly 
amenable to a stored-program computer. In the computation, combined use was made of 
integral arithmetic and decimal arithmetic but with careful error analysis. 


1. Introduction 


An essential part of the lure of algebraic number theory consists of tabulations of almost 
unpredictable irregularity. The advanced theory has somehow surpassed the scope of numeri- 
cal data without perceptibly improving the power to predict such data. As a result, for a 
small calculation an algebraic number theorist would not forego the personal satisfaction that 
results from the many ingenious, fortuitous, and deeply meaningful combinations of integers 
that he would encounter. Yet to really appreciate these vital irregularities he needs longer 
tabulations, with the inevitable requirements of uniformity and efficiency. Here the modern 
electronic digital computer can be used profitably. 

The discussion of the problem that follows will be primarily from the point of view of the 
computer program rather than from its theoretical development. Thus, for instance, the field 
will have to be regarded not as a set of elements satisfying certain axioms, etc., in the manner 
of Dedekind and Steinitz, but as an algebra with a definite basis and definite “structure con- 
stants” for multiplications. As a matter of historical fact the exhibiting of a basis was gen- 
erally far from trivial and lagged considerably behind existence theorems [/, 10] # 


2. Description of Fields 


For any prime /(=3m-+-1) a eyelic cubic field will be defined. Our interest centers around 
algebraie integers of such a field, written as follows (owing to the presence of the so-called normal 
basis): 


aw, + bw. t+ CW (a,b,c). (1) 


Here a, 6, ¢ are rational integers (or coordinates), and @),, @, #3 (or the normal basis) are the 
three roots of the defining (normal) equation {8}. 


F (we) ='+ w—mwao—n=0. 


Here » is defined as follows: First, 4/ is decomposed (uniquely) into the combination 
4l= A*®+-279’, 
where the signs of .1 and g are specified by the conditions (in integral /) 


A=3k- q>0. (4a) 


Phe computing time was provided under contract DA-028-OR D-12332 with the United States Army Office of Ordnance Research. The 
first named author was then at Wayne University, Detroit, Mich., and the second named author was at the Ballistics Research Laboratory, 
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Then, finally an integer n always exists, for which 
27n=(A+3)/—1. (4b) 


All algebraic integers under discussion are read. 

The most general basis of the field would be given by an integral unimodular combination 
{=D Jai; 1,j=1,2,3, detla,; tl. There is essentially no other basis [9]. 

The particular designation of roots is standardized as 


Ww) W» Ws. \e)) 


Lad 


Thus we can distinguish among conjugates & (given above), and &’, £’’ given by 


t’ Cw) + da. bw. (c,a,b) (Oa) 
&” bw, t CW + Aad, (b.e,a) (6b) 
(formed by cyclic permutation of w,, w:, w;), so that &’’’=£, ete. As a matter of notation, the 
subscripts denote conjugates only in reference to the basis w,. Later on &; shall denote (unre- 


lated) integers of the field. 
It was mentioned earlier that the discriminant of the field or the determinant 


Ww) Wey WwW 
Ww Wes Ww é) 
Ww WwW) W» 


has the value /?. This fact (which can be verified directly from eq (1)) is not used directly 
in the machine computation but uniquely identifies the field in question |6]. 
The multiplication of any two numbers in the field can be expressed in integral arithmetic; 


&, (a,.b,.e, &, (d,b.,€2), (Ss, 


then one can find the explicit representation of ££ in terms of the normal basis by means of 
the structure constants @,, B,, ¥;, @, Bo, Y2, defined by 


W) a), @) Bw» - V1 @ 
and 
Wo G3 QoW Bow Y2@ 
The explicit formula is, then, e. g., 
£&—(A,B,C), Qa) 
where 
A= a,a,a, vb be t+ BiC\C2-7 B.(a,b. ' yb, ) t a (byes ; h.e,)  Yo(C o> Co)), 
B 6,a,a, Tt a,b, - ¥1C C2 (a,b. 1 dsb, ) t 8. (bye. T be, ) Ato (Co C2M,), (9b) 


C= 7102+ By by by + ayeyes + ay (Ay by + dob) + yo(byeo + bre) + Bo(e,d2+ eA J 


This product can therefore be computed in the following tabular form: 


a). a B yi 
bb, 71 Oy ey 
C102 pB 71 a 
ayb, t ab B, ¥2 Oy 
bi bic Oy B ¥2 
ea ay yo Oy By 
| B ( 


These structure constants are given explicitly by the formulas: 


m+k+1 g+m— ay g>m— a 
a 2 ) B» 9 ’ V2 5 ’ 
(9c) 
q mm a q Mt a 
a a Mi l> B, 9 ’ Yi > i 


They are listed in table 2 (see p. 166). 


In particular, the norm of & is defined for (1) by the integer 


N(&)=&t’E” n(a+b+e)'—li (k—g+ 1)Xa*b/2+ (k4+-qg+1)Zab?/2+ abe}, (10a) 
and the trace is defined by 
t(é)=£+ + &” (a+b-+e), (10b) 
so that the defining equation of & &, or &” is 
3 ___ #(£)e*-+4 #(Et’)E—N(E)—0, (10) 
with the further designation 
#(£) t(éé’) m(a+b+e)*?+/ Zab. (10e) 


Note that if a+6+c¢=—0, then N(é) is divisible by /. Thus, applving this information to 


- -/ . } . > 
eed - \\a—C) \O—C),\C a)), 
we find N(g—£’) = (g—£’)(&’ — £”’) (£”’ —£) is an integer divisible by /. 


3. Units 


Units are algebraic integers of norm equal to +1. Thev are written as 
| ; 


(11) 


Ui] U1 @, 7 UsWe-> UsWs (U,,Uo,Us), 


, ) 


The most general unit has the form + 9%n{’, where a 


where, as before, u,, u., Us are integers. 
and 6 are (positive, negative, or zero) integers, whereas mo is a specially designated fundamental 


unit |6, p. 19). The only fundamental units are +n%', +(n{)*', +(n2)*'. In what follows, 
a. 


t 
units will be normalized (by a change of sign) to have, conveniently, only norm 
(U4) ,U2,Us), We obtain a cubic equation in 
these three latter unknowns, which we are in effect solving. Now such an equation (unlike, 


say, the Pell equation), has very little intrinsic interest. The interest in units is derived en- 


If we take eq (10a) and set N(é)=1 and (a,b,c) 


tirely from their role in unique factorization. 


To give an example, when /=19, one finds (as a byproduct of the main calculation) that 


N(@,) = N(2@,—.) =7. (12) 
This does not mean that 7 has two different factorizations, although seemingly 
7 W) WoW, (2a; Ww») (2a ws) (2ws @)) (12a) 
The reason,’of course, is that the factors can pair off as associates, or 
(2. Ww Ww) Wy 2w»), (12b) 
a unit of the field. To see this, in (12a) calling one set of factors & &, &/ and the other 
f, £5, £2, respectively, we must verify that one of the ratios §/&,, &/&, &/& is a unit. Thus 
for example, in accordance with the rules of section 2, we compute 
g geiey (A*,B*,0%) P 
~ ——— = (13) 
1 18181 N(E,) 
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and find that A*, B*, and C* are now divisible by N(é,). In practice, the whole computation 

would be performed modulo N(é,)(=7 here), using the structure constants of table 2. 

Thus the preliminary step in understanding factorization in fields becomes the recognition 

of units. Traditionally, one creates, by trial and error, combinations such as (12a), and one 
hunts for units, usually with better luck than one can completely explain [4]. 

It is an algorithm for creating a 


in the process, a unit, indeed 


The method used here will be very systematic. surfeit 
of algebraic integers with small norms (including norm unity). 
a fundamental unit, is inevitably produced. (Unfortunately, machine limitations forced the 


weakening of the method, as is explained in section 6.) 


4. The Sign-Discrimination Algorithm 


$ matrices of algebraic integers |2, section 4}: 


One operates with so-called reduced 3 
- i ie 
s %. so Si 
® &) g S &, \ l 421) 
SS sk 
s s So s4 
The columns represent conjugates of algebraic integers of sum zero 
14b) 


ier 
Str 
Sty 


c 
82 


form a basis. The further, and most vital, condition is 


tr 


and such that any three of the four & 


that the sign pattern be 


sgn @ 


either as @ now stands or under some rearrangement of columns. 
The algorithm is a method of generating such matrices in chains. A matrix @, is said to 
be a neighbor of ¢ if it is formed by adding one column of ¢ to another column of @ and at the 

Thus 


same time subtracting it from a third column of @ (the three columns being different). 
there are 4! possibilities to consider corresponding to the distinet triples (),j.j,), where 1 <j <4. 
For instance, the operation 


O|341\o Lda 


can be understood to mean that the ‘third column of ri) is Increased by the fourth column, which 


in turn diminishes the first,”’ or, writing just the first rows, 
Q (1, §, &3, & } 
L5b) 
?) (g Es, &2, € Es, &4) ) 


(The relation (15a) can be read backward as an operation on ¢,.) Of the 4! possible neighbors, 
either 3, 4, 5, or 6 will be reduced, depending on inequalities satisfied by the elements 
10). 


12. section 


We next define the conjugate ¢’ of a matrix ¢@ to be the matrix formed by replacing the 
first row & by its conjugate £; and forming the remaining rows from £; instead of &. The 
conjugate ¢’ will not be reduced as it stands, but a rearrangement of the columns (in this case 
the interchange of first and second) will render it reduced. Likewise, we define ¢’’. Note 
@ egain. Finally, two matrices ¢, and @ ere said to be proportional to one another if 


srr 


0 
for 7 


o,. >, or @;, and numbers 7, 7’, 7” 
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7 O 0 


@=|0 » O ip. (16) 


It is easily seen the. n, n’, n’’ are necessarily the conjugates of a unit, possibly 1, having 
norm +1. (In sone cases ¢, is proportional to itself, leading to interesting possibilities related 
to ramified primes but bevond the scope of the present discussion {2, section 22].) In all cases, 
proportional matrices will have the same set of norms (possibly under rearrangement). 
The algorithm operates es follows: The initial matrix is taken for convenience (writing 


only the first row) as 


dy = (w1, Ws, We, 1), n>0O, (17a) 
with norms respectively 
(n, @, ®@, 2), (17b) 
Or 2s 
Py = (W) + We, W3+@, w,, 1) n<O0 (17¢) 
with norms respectively 
n+m, n+m, n. 1). (17d) 


Branching out from @), one forms all reduced neighbors of matrices present, terminating e 
branch when the matrix is proportional to one alreedy present. Indeed, the number of reduced 
matrices is finite to within proportionalities [2, section 8]. The factors of proportionality 
ny are units among which can always be found the funde rental unit |2, section 18}. 


5. Modified Algorithm 


The originel algorithm, as just described, is e branching algorithm in which each @ mey 
have more than one successor, and each @ is compared with all those preceding. Because 
each @ contains 12 components and there could be more than 50 produced, the internal storage 
of the EDVAC would be taxed before the algorithm had gone very far. The branching 
algorithm was therefore modified to a forn in which only one successor is chosen for each ¢, 
and the comparison (16) is elwavs between @) and the most recent ¢. 

The modification consisted in trying only one-third of the 4! possible neighbors, nanely, 


1123], [132], [423], [432], [214], [241], [814], [841] (18) 


characterized by the fact that the second and third indices are either 2,3 or 1,4 in some order. 
Either one or two neighbors of this type occur, but the EDVAC takes only the first that occurs 
in the program, thus we have a chain-type algorithm rather than a branching algorithm, with 


1 
the further property that >) &; decreases at each step [2, section 17], whereas the corresponding 
l 
sums for the second and third row (conjugates) increase or remain the same. This assures 
us that when the inevitable proportionality occurs, the unit will not be 1 (but it need not be 
a fundamental unit either). 
As an illustration (see table 1)‘ the matrix sequence @, is reproduced for /=19. It is 


[241] [132] [314] 3/314] [314] 5[423] [341] o{123] [123] [123] [432] Oy(214] dy. 
The EDVAC discovered that @¢ is proportional to (some rearrangement of) @». The specific 


nuperical values are of some interest. We therefore calculate 


A full explanation of table | is given in section 7 
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2.507 O18 645 
Ww» 1.221 S876 162 
) 


2.285 142 481, 


which (except for minor modifications in section 7) is the accuracy presented by the EDVAC. 


ry . 
Thus 


dy ((1.0.0) (O,1,0) (0.0.1) ( 1.—1.—1)) 
or 
2.507 O18 643 2.285 142 481 1.221 876 162 l 
dv 1.221 876 162 2.507 O18 643 2.285 142 481 l 
2.285 142 481 1.221 876 162 2.507 O18 645 | 
o)» ((7,5,5) (4,—3.6) (2.6,—1) (—13,—8,—10)) 
or 
0.014 O37 286 0.017 151 S28 0.032 O77 205 0.035 191 747 
\» 7.445 752 324 17.010 034 653 18.661 625 853 9.095 343 524 
9.570 284 962 23.992 S82 825 11.693 703 O58 21.869 464 729 


The EDVAC then “recognized” that the first column would serve as a proportionality factor, 
i. e., If we rearrange the columns of ¢, to form 


do}, ( ( 13,—8, 10), (2,6, 1), (4,—3,6), (7,5,5)), 


then 
n 0 0 
0} 0 » O Ido, 
> @ ae 
where 7=0.014 037 286, 7’ 7.443 752 34, and 7” 9.570 284 962. In this case the 


decimal accuracy would be easily sufficient for “recognition,” but integral arithmetic was used 
(see section 8). 

The unit produced, 9=(7,5,5), is not fundamental, but, early in the course of the compu- 
tation, for @, & happened to come out as m= (—1,—1,0), a fundamental unit (see section 10). 
Actually, 


n= (no)*/no. 
This can be easily verified from the multiplication scheme (9), e. g., (mo)? 6,—3,—5), 
(no) (2.3.8) 
6. Numerical Accuracy 


As the algorithm was first conceived and applied in experiment [{2, section 11], decimal 
accuracy seemed sufficient because the svmmetric functions could be correctly evaluated by 
the nearest integer when necessary. In a long computation, however, the algorithm can go 
astray in many ways: 

(a) The algorithm may incorrectly discriminate the reduced neighbors from the others. 
This is bound to happen because the subset of the reduced neighbors chosen by the modified 


é, approaching zero monotonically. The error in the quanti- 





1 
‘| j _ , j » “=. 
algorithm just mentioned has : 


ties —, can soon be bigger than the quantities themselves, making a simple discrimination in 


sign untrustworthy. 
(b) It may make uncertain the decision as to whether @ and $* are proportional; for 


1; 


. . Pa 
l 


example, it may have to be decided that (using an obvious notation), &/é —s - 


wt, 
Str 


. 
3 


Str 
Str 
Str 


y S34 4 


in which all the quantities involved are in error. 
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(c) It may of course make the approximations to the all-important resulting units untrust- 
worthy; for these are the ratios discussed in (b). 

A vital and desirable modification, to control the accuracy of the computations, consisted 
in carrying out all arithmetic operations exactly by using integral arithmetic; (explicit formulas 
for the integral arithmetic operations on two such integers are given in section 2). Furthermore, 
all quotients were removed from the algorithm (see section 8). 

It would also have been quite easy in principle to have handled the sign discriminations 
in integral arithmetic. For instance, to tell if two algebraic numbers & and & agree in sign, 
all we need have done is test to see if the three conjugates £&, &, and &7& are each positive 
(total positiveness), which, by Descartes’ law of signs, is a matter of seeing that the defining 
eq (10) for (&&) has (integral) coefficients with alternating sign. At the time, however, the 
EDVAC internal storage was not sufficient to handle the enormous integers, so that an auto- 
matic error analysis was used instead. This was the only departure from integral arithmetic. 


7. Input-Output and Error Controls 


With the exception of the following operations, the problem was run completely internally, 
with no human intervention: 

(a) For each /, the strueture constants in table 2 were computed in advance by hand 
(merely to save space in the memory) through the a, column. (The entire table is given for 
convenience of hand checking.) 

(b) For each /, the zeros (and maximum errors) of the polynomial (2) were computed in 
advance by machine. The polynomial was scaled by dividing its zeros by m and carrying 
through the computation in fixed point. The root-finding method for the scaled polynomial, 


. I ay Ff ) n 
fw)=(2) +5 (2)-a($)-3° a9) 
m m m m m m 


is the bisection method (Horner’s method in binary), for appropriate isolating intervals. The 
maximum error in the coefficients of the sealed polynomial, for EDVAC, is 2~“, whence the 
maximum error in evaluating the scaled polynomial is 72>“; if, then, we use e=2~", the 
machine will provide an automatic error analysis in solving for w, by finding the zeros rj of 


f(w)+2e, ry of f(w)—2e, and setting 


w hie ot ri r; 


m yA é 


Because the numerical values of the zeros w,; are only used in the main algorithm to discriminate 
the signs of linear homogeneous expressions, d@,+bw:+¢w;, the scaled values w,/m and their 
error bounds e, were used directly. 

(c) The input and output of the algorithm itself was by means of IBM cards, in decimal. 
(This required internal conversion on input and output because the EDVAC arithmetic is 
binary. For each / two IBM ecards were read in. One contained the integers /, m, g, a, n. 
The other contained the integer / together with the (10-digit) decimals @;/m, w:/m, w;/m, and 
the errors in @), @2, ws. 

The output cards were as follows: First, there were copies of the two input cards; 
second, for each @,, as soon as it was produced, four cards were punched with 7=1,2,3,4, 
respectively, indicating the column of ¢, (prior to possible rearrangement) being described. 
The algebraic integer &, was described by the integers /, i, j, a), ;, ¢;, N(é), t(é), t(€), where N, 


t, ? refer to the defining eq (10) and 
E = jw, +b ja + ¢ jw. 


s 


Third, there would be a card containing the integers (see (11)) /, 7, (5), uw, Ue, uy, N(E)(=1), 


t(é), #(£) if the elgorithm runs to completion to produce a ¢, proportional to ¢o with factor ». 
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TABLE | Input-outpul cards 


lhe material in this table is presented in the form in which it is printed out by EDVAC 


INPUT CARDS 


Struc- l m q h om " 
ture 
con- 
stants 19 6 


Roots l w, m wo m ws m Error w)/m Error w: ” Error w 
and 

errors 

(x 10>" 19 1178364406 2036460270 3808570802 l 


OUTPUT CARDS 





< WwW Ci@ 
Matrix Matrix Column Viz i(g 1(é 
cards number number 
a h c 
(£, 19 0 ! 0 0 7 l 
(£. 19 0 0 0 | 7 I 
é 19 0 0 l 0) 7 l 
(&s 19 0 | | | l 3 
£1 1 l 2 l l 1] } 
g, 10 l l | 0 | 2 
ny 19 | 0 | 0 7 l 
(£4 19 l | l | ] 3 
(g 1 y 2 y | | 5) 
g, 10 2 | 2 0 | 3 
\g 19 2 0 ] 0 7 | 
(£4 19 2 | | | l 3 
( 19 3 2 2 | ) 
(£2 10 3 | 2 0 ] 3 
g 19 3 2 3 | 11 6 
&4 AD 3 3 3 2 7 S 
(g iv ! 2 2 l | i) 
(£ 19 | | 2 0 | 3 
£ 10 } | a) 2 7 11 
(£4 19 ! a) 7) 3 11 13 
(£, 14 5) 2 2 | l y 
(£» 19 a) | 2 0 | 3 
(§ 19 D 6 7 3 ! 16 
(&s 19 5 7 7 | 7 Is 
(g 1 6 2 2 l l 7) 
(Es 19 6 | 2 0 l 3 
(s a 6 7 y 3 19 19 
(4 19 6 Dy Q | | >| 
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TABLE lI. Input-output cards—Continued 


OUTPUT CARDS—continued 


The algorithm is then applied to the next case automatically. A complete set of input-output 
cards is listed in table 1 for /=19. 

On the other hand, error controls could stop a case (or the current value of /) before it has 
run to completion in the following ways: 

(a) The first error control computes the numerical approximations to &, £, &, & from their 
components a,, 6,, ¢, and values of w; Thus the maximum error in £;=d,ja;+6,@.+¢,a; 1s 
e— a, 6+ |b, & + \¢,\6, so that the (sealed) condition €,/m<e would mark the sign discrimination 
on gas untrustworthy. In this case an extra card was printed containing /, 7, 5, £4, &3, &*, &, with 
<* being the numerical approximation. The computation proceeded to the next / automatically. 
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(g Aw, bw + Ci@s 
Matrix Matrix Column N(é;) t(é,) T(é;) 
ecards number number 
a h ‘ 
‘ 19 7 l 10 11 5 7 26 29 
£ 10 | 2 l 2 0 l 3 16 
(g 19 7 3 | 0 l l 2 5 
(£4 19 7 } Ss 9 i l 21 14 
Be) S l +) 9 5 11 2 75 
é 19 Ss 2 l 2 0 l 3 16 
: 19 8 3 0 2 7 14 
bs 19 S 1 Ss 9 } l 21 14 
£ 1 ” l Ss 7 5 7 20 SY 
£, 1” 9 2 l 2 0 l 3 16 
g 1? o 3 l } ] 11 4 115 
&4 1% 9 } Ss 9 } l 21 14 
c Be) 10 | 7 5 5 l 17 71 
g, 19 10 2 l 2 0 l 3 16 
© 19 10 3 2 6 l 7 7 218 
(£4 10 10 } S 9 } l 21 14 
g, 10 1] ] 7 5 5 | 17 71 
fg 10 1] y 3 S | 11 10 353 
£ 19 11 3 2 6 | 7 7 218 
&4 19 11 ! 6 3 a) l 14 21 
£) aL 12 l 7 5 5 l 17 71 
(£, 10 12 2 } 3 6 7 7 408 
é 19 12 3 2 6 l 7 7 218 
ba 19 12 } 13 Ss 10 7 31 200 
” iW) foW Maw 
Unit Matrix Ladicator N(») t(») T(n) 
eard 
i Mo { 
10 12 4D 5 7 5 l 17 71 








(b) The second error control checked each computed norm to stop the machine if any neg- 
ative norms appeared or if any ¢, failed to have one of the eight neighbors (18). (None did.) 
(c) The third error control checked the imminence of overflow by testing to see if the 
largest term appearing in the norm calculation, n(a+6+¢)*, becomes too large (>2") to 
fit in the registers. In case of this eventuality, the machine proceeded to the next case, 


printing no signal card. These controls required very little storage. 
8. Neighbor Formation and Comparison 


What remains is now to describe the portion of the program that is strictly internal, and 
which, by itself, would not come to the attention of the machine operator. The two main 
steps of the algorithm are neighbor formation and comparison. 

The neighbor formation consists of the formation of ¢,;,, from ¢, and, ultimately, the replace- 
ment of ¢, by ¢,,,. Here, essentially, each neighbor operation of type (18) was tried in succes- 
sion on @,, and the sign discriminations were made to test if the pattern (1l4c) were valid. 
Here the numerical estimates were kept in double precision, the integral and fractional parts 
being stored separately (each consisting of at most about 13 decimal digits). 

The comparison of ¢ and ¢, is more complex. First of all, the cases /=7 and 13 were 
excluded by virtue of the fact that /=1 (so that +, or —@, isa unit anvway). In the other 
cases (/=19), the set of norms in wp (see eq (17b,d)) contains exactly one norm of value 1. 
Thus a ¢, is not conceivably proportional to ¢, unless some column in ¢,; is a unit (say) 7. A 
further necessary condition is that all the norms in @, match those of ¢@,. When this wes satisfied 
the machine multiplied the top row of @, by n, 7’, and n’’, in turn (using (9a,b)) testing at each 
stage to see if ¢, or a conjugate was obtained. (In this last step, the possible column rearrange- 
ments are fewer if there is only one norm in ¢, equal to 1.) 

It is seen that the machine part of the algorithm at no time involved division (even such 
minor calculations as division by 2 were performed by shifting the binary numbers) 


9. Supplementary Hand Computations 


The restriction that the whole process be contained in a memory of 1,024 cells proved to 
be very great. As a result, in the EDVAC program provision was made only for double 
precision numerical estimates (the integral part and the fractional part being stored separately ) 
and single precision integers. Thus only 2 few cases ran to completion, but only 18 cases 
failed to produce at least one unit in some ¢, (with or without giving rise to proportionezlity). 
In these cases a hand computation produced units from the algebraic numbers automaticelly 
computed. The general scheme for these hand computations was the following: Suppose the 


machine has produced two nonconjugate numbers &, & with same norm, preferably 2 small 


must be associated with a conjugate 


prime p. If p is uniquely factorable in the field, then 
of &, as we can test in the manner described earlier in section 3, completely in integral erith- 
metic. If the norm is not prime, more complicated combinations must be tried, but they are 
no different in principle. 

At least half the hand-computed units were obtained quickly in this fashion. The cases 
that did not yield to this method directly did so when more than two integers were used in 


proper combination. 


10. Fundamental Units 


In each case those units 7 were noted for which 
max t(),t\n) (20) 


over all units produced in the algorithm was a minimum. In accordance with arguments given 
by Hasse [6], this would be the fundamental unit if one were present in the list. 

By a remarkable coincidence, almost simultaneously with the completion of the runs 
(Aug. 1954) Peter Swinnerton-Dyer kindly communicated an independent table giving the 
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traces f, ? of the fundamental units for 7</<300. A comparison of the tables showed perfect 
agreement except in the three cases, /=373, 379, 463. In just these cases the (hand) computed 
units had positive traces for ¢ and 7, which should have suggested that the units computed 
from the EDVAC output are squares. 

In the listing in table 2, the components and traces of both a unit and its reciprocal are 
shcwn in each case (here t(n)=t(1/n)). In those cases where the units were computed by 
hand, (*) is shown. In general, regardless of how the units were produced or derived, the 
symbol “f” indicates the fundamental unit (and reciprocal) and “ff” is reserved for cases where 
both fundamental unit and reciprocal were produced by machine. In the three exceptional 
cases, the units, denoted by “‘s’’, are squares of fundamental units. 

In order, however, to give some kind of reasonable impression of the amount of labor, 
particularly in view of 13-digit coordinetes (!), we should point out that the number of trials 
in testing fundamental units is easily 


0 | logimax ¢t(n),t(m))/log />*. (21) 
To see this very quickly, note that if 7 is not a fundamental unit, we can write (recalling 
nn’n’’=1), 
n=no(ny)’ 
n’ (m,)“(m)” No ”(no)* ® 9 (22) 
n” = (no)*0 = 10 °~ * (m0) ~* J 


where ny is some (as vet unknown) fundamental unit, and @ and 6 are rational integers. By 
taking reciprocels, or conjugates if necessary, we can assure ourselves that a>b>0. Then 


a(e’)*=®,", (22a) 
where 
M=a?—ab-+ 6. (22b) 
Letting max (9, 9, 9” )=—H, mex (|, \, int) = Mo, we find 
>i. (23a) 
or since (a+6)/AMf <2a/M <2a/(4a7/3)=3/(2a), 
H,<H°’. (23b) 


But, by the concluding remarks of section 2, 


l< |no— no! !no— no! |\No— 20 | << SH? 


and 
H?< 9° + (9’)?+ (n”)?=t(n)?—2t(n) 
thus 
; 3 log H leg 2 [max ft(n) . 7a 
0<b<e<—, -—_— <8 — | . » |, 
leit § ia log [//8] 


This verifies the order of magnitude (21).° 

The problem of showing the units 7 to be fundamental now resolves itself to checking on 
whether any of the units n(n’) is a perfect MW power, where a and 6 come from the finite set (24). 
But this check is somewhat laborious, even to program for the EDVAC, and owing to the 
convenience of Swinnerton-Dver’s table, it was never carried out. 


Much stronger inequalities are given by Hasse [6] using more specialized notation and information 
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11. Conclusion 


The EDVAC produced 176 pages of tabulations, too long to be reproduced here, including 
some 3,500 algebraic numbers with small norms, in less than 6 hours. With the precision 
provided for in the coding, only four cases ran to completion, but 25 of the 43 cases run produced 
at least one unit in the process (for /=97, 13 units were produced). The remaining 18 cases 
provided enough material to permit the hand computation of at least one unit. (For /—499, the 
unit required 43 decimal places to check by hand). 

In every case where g=1 the explicitly known fundamental unit (—(k—g+1)/2, 

(k—g+-1)/2, —(k—g—1)/2 was among those produced by machine, as well as the (non- 
fundamental) unit (—1,—2,0), except when /=37 and 313. 

In addition to units, mentioned earlier, the number of matrices produced after @» is shown 
in table 2, with the following symbols denoting the outcome of the computation (see section 7): 
“comp” meaning completed, ‘disc’? meaning incomplete (by failure of sign discrimination), 
“over”? meaning incomplete (by overflow in norm formation). 

It is clear, however, that to obtain the fullest benefit of this algorithm one would require 
considerably greater storage, both to perform the unmodified algorithm (guaranteeing a 
fundamental unit) and to permit sufficient decimal (or preferably integral) precision. 
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Estimation of the Frequencies of Thin Elastic Plates 


With Free Edges 


Tosio Kato,' Hiroshi Fujita,’ Yoshimoto Nakata,’ and Morris Newman 


A variational method is proposed for calculating the frequencies of thin elastic plates 
with free edges, with rigorous error estimates. As a numerical example, the fundamental 
frequency of a square plate with the Poisson ratio 0.225 is calculated with a satisfactory 
result, the possible relative error being less than 1/2000. Generalization to more complicated 
boundary conditions is straightforward. 


1. Introduction 


1.1. The present paper is concerned with the study of the vibration of a thin elastic 
plate with free edges. Let us consider a plate that occupies in its natural condition a plane 
domain J) in the sy-plane bounded by a sufficiently regular contour C. It is known [1,2] * 
that the problem of determining the frequencies of this plate may be reduced to the following 
eigenvalue problem, denoted by (Pr.), for the differential system consisting of the differential 
equation 


ey . 


A*w ( 52+ ay:) w=w in D, (1) 


with the boundary conditions 


OAw ; d Ow 


) 0 on C’, 
on + Mds Onodt 
>? (2) 
Zap 
wAw +(] Ma) ~ 0 on C’, 
on’ 


where uv is an elastic constant, called Poisson’s ratio, such that OSy< 4 (the mathematical 
theory applies for 0 Su< 1), and where s is the are length parameter of C and 0/0n, 0/Ot repre- 
sent respectively the derivatives in the directions of the outer normal vector n and the tangential 
vector ¢ to (at the point under consideration. In the case of angular points on C, where the 
direction of n varies discontinuously with a jump, w is subjected to an additional boundary 
condition {1}: 0%w/Ondf at a point P? on C tends to the same value as ? approaches any one of 
the angular points from either side. 

The main object is to propose a variational method that enables us to calculate approxi- 
mate values of the eigenvalues with a rigorous estimate of the error; in other words, to calculate 
upper and lower bounds for the eigenvalues. Such a method seems interesting and important 
from the viewpoint of application, inasmuch as it appears, even for most elementary shapes of 
(’, almost impossible to obtain an explicit solution of (Pr.). On the other hand, usual approxi- 
mating methods applied to (Pr.) seem to suffer from the lack of error estimation. For example, 
the well-known method of Ritz is applicable and was actually applied to the case of a square 
plate {2}. Although we know empirically that his method provides us with approximate 
values whose accuracy is satisfactory in most cases dealt with in practice, it should be noted 
that they are only known to be upper bounds for the quantities in question and nothing can 
be said about the error bounds. See [10] for a comparative discussion of existing methods, 


some of which furnish lower bounds. 


University of Tokyo, American University, and National Bureau of Standards 
Department of Physics, Uaiversity of Tokyo 
Figures in brackets indicate the literature references at the end of this paper 
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1.2. Our method is based on the following theorem, proved in [3], concerning a_ self- 
adjoint operator in a Hilbert space ‘ of the type 7*7, which reads, in a form slightly modified 
for our later convenience, 

THEOREM i* Let T be a closed linear 0 pe rator with its domain dense in a Hilbert Space 
§ and its range in a second Hilbert space §' and let T* be the adjoint of T. Let be&rX< 8 
(OSa< B) be an open interval containing at most one nondege nerate eige nevalue of the operator 
T*T but no other points of its spectrum. Tf u and v are elements be longing to the domains of T 


and T*, respectively, such that 


u <0. r|\ x0. (Tuwr)=(u,T*:) 20 (; 


and 


\” a (8 ! e- (4) 


(Tur) it Tu? T*r ? - 
n ’ =s he i — 7", ») 
uoeP 2 u|* v|!|? 


then there exists certainly an eigenvalue X' of the operator T*T in the interval (cé 8°) and the following 


hold, where 


ine quality is valid: ®* 


Ls’ su, 6) 
where 
é . . € . - 
L Liu, (1 8 -). / l (u,v) (n+, = 5 ‘) 
In order to make this theorem applicable to (Pr.), it is required to construct the spaces 


), 9’ and define the operators 7,7* in such a manner that the eigenvalue problem 7*7Tw= dw 
is equivalent to (Pr.) and the quantities in (5) are practically caleulable, at least for the w and 
rchosen. (Pr.) is nothing but the eigenvalue problem for a self-adjoint operator //7 (in a real 


Hilbert space 9 of square integrable functions on 1), defined as 
Hw A?w (8 


for w subjected to the same boundary conditions as those of (Pr.), in addition to certain regu- 
larity conditions making A’w well-defined Therefore, our problem is to reduce this pre- 


assigned operator // to the form 


H=I°T 9) 


with a suitable choice of 9’, 7, and 7*. Here, wu and ¢ may be called “trial functions,” or 
possibly “trial vector functions,” in estimating \’, a better choice of which will vield a better 
estimate. Finally, a and 8 are, loosely speaking, a rough upper bound of the next lower 
eigenvalue to \’ and a rough lower bound of the next upper eigenvalue, respectively. Such 
rough estimates will be obtained by other methods (for instance, the Ritz method for a’), or 
occasionally by solving an appropriate auxiliary problem, as in the case for 8? concerning a 
square plate, treated later. The decomposition of /7 into the form 7*7T is carried out in 


section 2 for a general shape of (. This is followed by some supplementary remarks on theorem 
‘ Concerning the theory of Hilbert space, we mostly follow the notations and terminology of Stone {4 For operators between two Hilbert 
spaces, see Murray 


An elementary proof of this theorem will be given elsewhere 


Any eigenvalue of 7* 7 is nonnegative, because ( 7" Tw.» rn 0 


e2 
Oand n+ 0 follow from (4 


170 


> 


| for improvement of the approximation. In section 3, confining ourselves to the case of the 
square plate bounded by «= +1, y= +1, we illustrate in detail the method of actual calculation 
along the lines mentioned above and give a numerical estimate of the smallest positive eigen- 
value of (Pr.), taking »=0.225, with a relative error bound less than 1/1000. Section 4 is 
devoted to concluding remarks, with a brief reference to other methods with the same object 
that are applicable to (Pr.). 


2. Decomposition of the Operator and Remarks on Theorem 1 


The first part of this section is devoted to decomposing //, the differential operator of (Pr.) 
including the boundary conditions, into the form 7*7 in the sense stated in the introduction, 
namely, to constructing 7 and 7* so that //=7*T7, i. e., T*Tw=de, implies (Pr.). Strictly 
speaking, the differential operators appearing here as well as in the sequel must be replaced by 
their closed extensions with certain extended domains, or in other words, differentiations should 
be interpreted in a certain generalized sense. Nevertheless, detailed discussions of this situa- 
tion will not be given, for they seem to be unnecessary for practical application, which is the 
ultimate object. 

2.1. We begin with the following heuristic considerations.” Let 7 be a closed operator 
from 9 into 9’ and 7* its adjoint. Then the eigenvalue problem 7*7Tw= nw is equivalent, at 
least under a suitable assumption on the spectrum of 7*7, to the variational problem 


56 Tu ?*=0, u!|?—1, 


where 6 is the first variation. Thus, 7*7 is the “gradient operator” of the quadratic functional 
Tu * in the sense that 6 Tu *=2(7*Tu,éu). On the other hand, as explained in textbooks 
11, S|, (Pr.) is derived from the variational problem 


beJ|u] —0, | urdrdy= 1, 
e D 


where J\u|, which represents (apart from a constant factor) the elastic energy of the plate, is 


Oru Oru ou \* 
J |u| 2( -- rdy, (10 
Ju! | | — i {2 oy” ( Or OY ) }] ate 


that is, /7 is the gradient operator of J{uj. Accordingly, © being the same as that introduced 


expressed by 


in defining //, it is reasonable to adopt 9’ and 7 such that 


Tu ?=J\ul (11) 





holds. In faet, we proceed as follows. Let L*(/):p) denote the real Hilbert space consisting 


of functions integrable (in square) on J) with norm defined by 


ul p widardy, 
D 


where p is a positive constant. 


In view of (11) and the facet that Jiu! can be written as 


J lu] AL (Au)?+(1—yp) {( a ) 4 ( se ) 4 2( cae )}] drdy, (12) 


we define §’ by 


0D = L*(Diyp) * L2(D31— w) * L2(D 31 — pw) & £7(D32—2yp), (13) 
Some general devic peftul in suc decomposition are proposed in [6 
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where represents the Cartesian product of Hilbert spaces. It should be noticed that an 


element v of ©’ is a four-component vector function r={r;,7,73.74; defined on DD. Now the 


requirement J[u!| Tu * is evidently fulfilled by the operator 7, defined by 
mm mn Ou ou OV 
Tu=Tu(r,y) Au, — =>? EH’. (14) 
R Or OY OLOY ' 


Elements uv belonging to the domain 2 of 7 are subjected to certain regularity conditions, but 
are otherwise arbitrary, i. e., free from any boundary conditions. Though these regularity 
conditions are too complicated to be specified explicitly in elementary terminology, it will be 
A 
enough for our purposes to know that the class C of functions uv, such that u, 0u/Or, Ou/Oy are 
continuous and 0*u/dz*, 0?u/Ordy, Ou/dy? are piecewise continuous on J)+C, is contained in 
A 
); namely, CC®. 
T* is defined (together with its domain D*C’) by the condition that 


(Tur) =(u,T*r) (15 


be valid for any w€ D. To determine the form of 7*, we note that 


~-. 


oi F Oru Oru Oru 
(Tur) (Au)r,+ 1 bh) t+ (1 bh) -U3t+2(1 ) r lrdy. (16) 
{[- ox? oy" u or dy ‘ Y 
In reducing the right side of (16) to the form u-T*y drdy, we can resort to devices similar to 
a/ D 
those usually employed in deriving the differential equation and the boundary conditions of 
(Pr.) from the variational problem concerning J|u|. Before doing this, we introduce some 


notations, some of which appeared in section 1. 
Let 
s=are length parameter, 
n=n(s)=(n,,.n,)—unit outer normal vector to C, 
t=t(s)=(t,,t,)=unit tangential vector to C in the positive direction, 


and note that 


° n ° nh, e , ° t 2 +t, ° ’ li 
on Ou OY eli or OY 
and also 
re) 0 O re) Oo 0 
: . ( 1S) 


or 7On’ 7dt dy “On 


Hereafter it is assumed that the components of » and ¢ are differentiable except at possible 


angular points P,,P., . . ..P, with s=s),%, . . .s,. Also it is known that, for any suffi- 


ciently regular functions Fi(z,y), Gr,y), the identities 


+ 


“OF 2 . OG ! . — 
k= Gdrdy =" | 2 or drdy+ | I ‘(in,ds, 


> (19) 


: G - 


e D 





4 


are true. 

Provided the components of ¢ are sufficiently regular that all procedures used below may be 
justified, the expression for (7u,r) in (16) becomes, by some small calculations that make use of 
(18) and (19), 
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a ’ Or, , Oo O*r; oO 
(Tur) fof nan 1. (] {opt oy t +2 roy y fete J u- Past | 5! - ‘ Qds+ | Se = Rds, 





(20) 
where 
or or 7 
P P(v,.s) 1 +(1 vl +n, +n 
{oem Tw Or wT ov Zz 
> (21) 
Q=Q (r,s) =r + (l—p){ rent t+rgnZ+2ern.n, 
P= R(v,s) = (1—p){ vont, +r gn,t,+r,(n,t,+n,t,) }. , 


Furthermore, because at a point ? with s=o we have 
ou du(s) 
(a7), “Cay ),.. 


and consequently 


"Ou ,, du 
Sp ids [Gt Ras, 


we can rewrite the last term of (20) by partial integration, taking account of the fact that C is 
closed and that 2, as well as the components of n and ¢, is differentiable with respect to s except 


at s=s,, .. ., &,, so that 
"Ou * aR 
R Xx — s—>) (Ss i, 22) 
cot 7” Jv ds a i =" dD. 
where D,=lim R—lim PR, (e==-1,3. .. ..%). 


The substitution of (22) into (20) results in 


re F Or, , Ors , Os ,, aR ou 
(Tur) [uf war -(] ee tat . } Jtcdy—f u( ij ~ 7, )ds+ Sn Qds— Suis, )D,. 
(23) 


This expression should be equal to (u,7*r). Because uw and Ou/dn on C are arbitrary, we 
immediately conclude that this is true if and only if ¢ satisfies the boundary conditions 





- ’ 
> 0 
/ = hi } 
Q—0, , (24) 
dD), QO, (7 ee * 
4 


Then 7* is given by 


——— Or, , OV O°) _ 
T'*r—T* 01,0250 3,04 wAr,+(1 i SS +2 de}. (25) 


Though an element rv belonging to D* must obey certain regularity conditions besides the 
boundary conditions (24), we will not discuss these in detail. Again it is sufficient to know 
that anv ¢ with components belonging to 0 and satisfying the boundary conditions (24) is 
contained in D* 

The construction given above shows that the operators T and 7* are certainly the desired 
ones, but it may be worthwhile to verify directly the equality /7=7*T, i. e., the equivalence 
of the eigenvalue problem 7*7w=)w to (Pr.). 
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Let w belong to the domain of 7*7. Then Tw belongs to D* and consequently is subjected 


to (24). But 


PT OAw ow Ow Ow Ow 
{ ue) + i + t. vi _i 
eon MY or oy YT Or2dy "Ody 


OAW 0 (ow 0 (ow 0 (ow 0 (ow 
- On rl ‘ ({2(% )n y dy (oz Jn} ' {2(3 )n + Oy (ay )n¥] 


OoAw 0/ow . Ow 
yu +(1—n) — (55+. ) 
On On \ Ou Ov 
Oo 
Aw 
On 
by virtue of (17 Similarly, we obtain 
(XT w) = wdw-+ (1 OW 2 OW 24 4 Ow 
uw we i°-+ ie ol i j 
v U Dr! } Sy? n; See’ } 
Awl ~se 
" 
Mt Li yn 
_ “Ww orn ow 
Ri Tw l u tind + n,t,- (nt nf. 
Ou OY Oroy 
ow 
. Oron 


Thus, w satisfies all boundary conditions of // and therefore belongs to the domain of 7/7. Con- 
versely, by a similar argument it is seen that anv element in the domain of //7 belongs to the 


domain of 7*7 Finally, we observe that 


Ten... . O Ow oO (Ow a & ow 
T*7T w= pA\ Aw) l— sp {2 ( )+ dy? (dye )4 27 (22)} 
Avw, 


as it should be 
2.2. The following remarks on theorem | include some devices that would be of use in 


obtaining a good estimate, particularly in constructing a desirable pair w and 
Remark 2.1. The estimates given by theorem 1 are sharp when & is small But it is 


known [3! that & can be written in an alternative form, 


. { | Tu " all ot | T'*, n i '); 


and & =0 occurs if, and only if, 


Tu n/ and T*r nu. 
i u 
This implies that 7*7u= yu, so that 7° is equal to an exact eigenvalue \’ of T*7 and uw is the 
corresponding exact eigenfunction w. Thus we may expect a good result when w and ¢ approxi- 
mate, apart from constant factors, wand Tw closely. This interpretation of wand ¢ should be 
taken into account, in particular when some knowledge concerning the behavior of w is given 
in some way or other, for instance from physical considerations or experimental results. As 
for precise discussions of the relations between u, 7, and w, reference is made to the paper by 


Kato [3 
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Remark 2.2. If wand ¢ are not completely preassigned but admit of some free parameters, 
it would be intrinsically the best way to choose those parameters first to minimize UV and then 
maximize 1. In practice, however, it appears sufficient as well as convenient for calculation to 
determine them by minimizing &. A standard form for such adjustable trial functions is 


u=au+au%+ ... +a,U ay 
(26) 

r= +pe+ ... +80, ) 
where uw? u®, 2. au”? are fixed functions belonging to D and vr», . . 0” are fixed vector 
functions belonging to D*, and where aj,a2, . . .,a@m, 6,82, . . .,8, are free parameters. If we 
calculate the quantities in (5) with these wand ¢, then | u)?, | Tu)? become quadratic forms in 
‘a;}, 0 *, T*r *® quadratic forms in {8}, and (7u,v) a bilinear form in {a;} and {6}. Thus é 


becomes an algebraic expression in those parameters whose minimum is desired. There is 
little difheulty in the necessary computations, especially when we can resort to an automatic 
computer, even if m and n in (26) are not very small. 

Remark 2.3. Suppose that a closed subspace I of S reduces the operator H=T*T, 
that is, let //we? for any w belonging to the intersection of Mand the domain of 7. Then we 
can regard //7 as an operator in the Hilbert space 2 denoted by H7y. It is known that eigen- 
values and eigenvectors of //y, are also eigenvalues and eigenvectors of /7, so that we may speak 
of an eigenvalue of /7in WM. If the eigenvalue ’ of /7 in question turns out to be an eigenvalue 
of //y,, we can derive an estimate for \’ from theorem 1 more conveniently by restricting 7 to 
WM, and 7* to a suitable closed subspace Y* of H’, with the properties 


TMOD) CM, T*MIOD*)\ CM. (27) 


This procedure is particularly advantageous when a degenerate eigenvalue of the 7 in question 
is converted into a nondegenerate one of HZy by a suitable choice of WM so that applica- 
tion of theorem 1 becomes possible, or when the neighboring eigenvalues of 7g, »re more widely 
separated than those of /7, thus permitting us to take smaller a and larger 8 so as to make the 
estimate better. Moreover, though the smallest eigenvalue of 77 in our problem is equal to 
zero, the associated eigenfunctions being of the type ar+by+¢, we can sometimes find I such 
that the eigenvalue \’, originally intermediate, becomes the smallest eigenvalue of H/g, and 


hence its upper bounds can be obtained by 


~ 
= 
; 
(-) 


—-2nr’, (we M ), (28) 


the so-called Rayleigh principle. This seems more convenient than the estimation by the right 
inequality of (6), for it does not involve v, an element of D*, for which a good choice is not easy 
because of the complication introduced by the boundary conditions. Therefore, we may 
prefer [’, to (as an upper bound for such an eigenvalue. 

Nontrivial 92, W* as mentioned above are often obtained by considering the symmetry 
properties of the operators 77, 7, and 7* reflecting those of C, as we shall illustrate in the 


succeeding section in the case of a square plate. 
3. Vibration of a Square Plate 


In this section the general results obtained so far are applied to the case of a square plate. 
Let ( be the square with vertices (1,1), (—1, 1), (—1,—1), and (1,—1), which is the same as 
that treated by Ritz [2]. Bounds will be calculated for the smallest positive eigenvalue of = T7*T 
according to theorem 1 and the Remarks in section 2, i. e., lower bounds by means of (6) and 


upper bounds by means of (6) or (28). 
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3.1. We begin with a preliminary consideration. — Let 


J {u] (l1—yp) (uz. uz \dady. (29) 


/ D 


Let the self-adjoint operator 7 in & be the gradient operator of J[u|. Explicitly, any element 
# in the domain of // is subjected to the boundary conditions 


o_o 


, 0, (@@=+41,-1<y<1)| 

Or Ox” 
30) 

ow Ow 0. (—1<r<l,y 1). | 

oy" Oy" 
and we observe that 
: O'w . Ow _ 

Tlu l lu ( dri? oy' ) 31 


The eigenvalue problem // #@= dw is reduced, by separation of variables to the eigenvalue prob- 
lem 
du 


dr ku, r 1),| 


to 


ee | 


and a similar one in the variable y. Equation (32) is nothing but the eigenvalue problem treated 
in the analysis of vibrations of a free bar of length 2, and is completely solved in textbooks 


[7,8]. If we denote the eigenvalues by &} (n=0,1,2, . . .) and note that 0 is the unique degen- 


erate one with eigenfunctions a+ br, we have 


k= k,=0, 0< kh.<k<k, 


where ® 





>> ) 
kt. =2.3650204, k. =3.9266023. 
kk, = 5.4978039, k.—7.0685828. 
The normalized eigenfunction u,—u,(2) associated with is given by 
7 ‘ 
u , Mu 2 


H, cosh k, ent hyd a k ‘cosh k, 4 for eVeh + 4 > 34) 
\ cosh? k, +cos? k 


sinh &,-sin kir+sin k,-sinh kis 


u, ; - , for odd n>=3. 
ysinh? k&, —sin® k, J 





It should be noticed that w, is an even or odd function, according as 7 is even or odd. 
The eigenvalues and eigenfunctions of // ere given by 


Amn (1— yw) (ki, 4-4), , Un (a )u,(y), fe ES 3 > eee | 35) 


We note that @,,,. (m,n 0,1,2... .) constitutes a complete orthonormal! system of eigenfunctions 
of the self-adjoint operator /7/. 

3.2. Comparing two quadratic functionals Jlu} and J[u|, or two operators 7/7 and 71, 
we can not only determine the constant 8 in theorem 1 in a rigorous fashion, but also the sym- 


k, is obtained as 1 root of equation tan / tanh & or tar tant wcording as s odd or ever 
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metry of the eigenfunction associated with the smallest positive eigenvalue. We note that 
J\uj sJ\u). (36) 


From (36) it follows, by virtue of the properties of /7 stated above, that the spectrum of H 
consists only of discrete eigenvalues, and also that, if we number all eigenvalues of HW and H in 
ascending order, counting any degenerate one repeatedly by its multiplicity, we have 


lA 


pm (37) 


where Ay and Ay are the Nth eigenvalue of //7 and the Nth eigenvalue of //, respectively [7]. 
Thus the eigenvalues (35) of 7/7 are the lower bounds of the corresponding eigenvalues of H, 


and hence we can take 6 as 
B- ,, " (1 u) (4, 1 he R (38) 


if the eigenvalue \’ of 77 in question is the rth eigenvalue from below and X,,, is the (r+1)th 
eigenvalue of 77. 

We turn to the study of the symmetry properties of 17 and 7/, i. e., of closed subspaces 
specified by the symmetry properties of their elements, and to the reduction of each of these 
operators. For the square plate under consideration, the boundary conditions of /7 follow 
from (2) and the remark there and take the form 











cass Se? 5 
; ort H) >, 0, (#=+1,-1<y<b), 
J te Oy" | 
T Oem 2a) 
O ow 19 P Oru 0. (y 11 —1lor<}). 
Ov L Oy" Ou | , 
Ww orw 
— at a r= +1,—1<y<), > (39) 
Oru Oru 
5,2 + art 0) y= +1, —1<z2<1) 
- 7 
ow 0. (r a y | | ). 
Oroy J 


In the deduction of the last condition of (39), use is made of the fact that, for instance, 


ow ow ; ow ow 
on the side r=1 


: - on the side y=1, 
Ondy Oroy on ot Oroy y 


and hence that the condition 


orn ; o-w ow 
lim lim means 0 at the vertex (1, 1). 
on of On Of Oroy 


! ~! 


Owing to the symmetry properties of (39) and those of the differential operator A’, His 
reduced by the following subspaces : 
W(o,0)=set of functions odd both in » and in y, 
Woe) =—set of functions odd in « and even in y, 
Wie.) set of functions even in « and odd in y, 
WMie,e)=set of functions even both in z and in y, 
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which are mutually orthogonal and whose direct sum is the whole space . We see immediately 
that //7 is also reduced by each of these subspaces, so that the inequality (37) still remains valid, 
if we regard 77 and // as operators in any of them. Let the eigenvalues of // in these subspaces 


be denoted “us 


Rat Dee Rex. . ; in Woo), 
Mo < M1 < pflo<.. either in Di(o.e) or in Wile), 
Vo <_ V1 <_ Vox. — In Wile 6) 


the union of which is the total spectrum of 77. On the other hand, the eigenfunction @,,, of 
with the eigenvalue X,,, belongs to W(o,o), Mo.e), Mle.o), or Wie.e), according as (m—odd, 
n=odd), (m=odd,n=even), (m=even, n=—odd), or (m=even, n—even). 

Because the zero point set Xt of // (the eigenspace with the eigenvalue 0) consists of fune- 
tions of the type ar+ by +c, it follows that 





Mio) LM, ‘ 
Woe) us uar 
\ 
a (40 
Wilea)yor Miu by 
WME eyom u;u=c}, J 
and hence that 
A, >, fo = Va =U. 
In consequence, the smallest positive eigenvalue of // is equal to min(Aygi.7). We can 


determine the order relations between Ay, u,, and », as follows. Because (40) means that 
is degenerate neither in Yi(o.e) nor in Wile,o), and v, is not degenerate in Wi(e.e), we note, 
according to (37), that Xyo( =») and X» are lower bounds of uw, and »,, respectively, and hence 


we have, by means of (33) and (35) 


“, = re dh» (] uu ks (] uw) X 31.28 : ) 
(41 
v, > dn =2(1 —p)ki— (1 —p) X 62.59. ) 


On the other hand, the substitution of the trial function «—asy 6 Wo) into (28) vields 
Ny < l pu) LS. (42) 


Combining (41) and (42), we conclude that the smallest positive eigenvalue of // is the smallest 
eigenvalue Ay in Wi(o.o). Moreover, owing to (42) and the relations A,,—0, Ayy=—~A (1 —yp) 
237.72 . . ., we obtain 


Nox A As), (433) 


which shows that \,, is the unique eigenvalue of 71 less than dX», and hence Xp is not degenerate. 


Also we have, by means of (37), AyyS A,, which enables us to take as 8 in theorem 1, with respect 
to \’=Ao, anv number such that 


B°s (1l—p) X237.72 . — $4) 


Furthermore, we divide Yt(o,0) into the following two mutually orthogonal subspaces 
W and WM’, 


VW u: “ue WMlowo), u(ay) =ulyya) 
We’ u: “we Woo), u(ry) u(y) 
Evidently, the spectre of 779% ind those of 7/99 ure lentical, and any eigenfunetion of oneof them is converted into an eigenftunctior 


of the other by permutation of. and 
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each of which is seen to reduce both of H and /7. 

From (43) it follows that A» is an eigenvalue of //7 not in WM’ but in M, beceuse the smallest 
eigenvalue of HT in W is ys with eigenfunctions Const -(@,,;— Ws,), though this is at the same time 
an eigenvalue of Hin M with eigenfunctions Const-(#,3+ Ws). In calculating bounds of Xp», 
to which we confine ourselves hereafter, we may therefore restrict /7 and T to M. 

3.3. Wego on to the construction of the trial elements '' u and vin theorem 1. Introducing 
l 


six elements u® of M, given by 


u?—P,(x)P,(y), u® =P (4) P3(y) + P3(2)Py(y), 
u P(x) P3(y), u =P (x) P;(y) +P; (2) Py), 
u P3(2)P5(y)+P;(2)P3(y), u® = P(r) P;(y), 


where 7’, represents Legendre’s polynomial of degree 7 (¢=1,3,5), we put, in accordance with 


(26), 

u=au'? + au t+ apu'® (45) 
with free parameters @,a,, .. . ,@ to be determined later. This u evidently belongs to 
MOT, because it is a polynomial and surely satisfies the regularity conditions of D. 
’;,2,03.04} is constructed as follows. We can easily verify that the 


-) 


The other element ¢ 
boundary conditions (24) of D* now reduce to the following forms: 


5 
Or; Or, Or, 
+ (1 = -+-2(1 0, r +1, lk y<cl) 
any HW)», ( Mh) dy (J y 
O/ Or, Or 
(1 2(1 ‘=(@ (<= ss 1) 
OV ad OY Hs, . 
4 (46) 
ue) (] u)ls 0, (7 +]. iz y: 1). 
ul, + (1l—p)r,=—0, (—l<r<l,y +1), 
0,.=9 (a ly 1) f 





In view of the symmetry of I, on the other hand, we can take as W* the subspace com- 


posed of elements ¢ with the following properties of symmetry: 


A 


ry, fo, tf, are odd both in « and in y, 


(47) 


r, is even both in s and in y, 
m(ry)=ry(yw), Po(L,y) Ms(Y,r), Mary) MY). 


Then the relations (27) obviously hold. 


A general element ¢ satisfying the conditions (47), and whose components are polynomials 


in v.y with the degrees shown below, is expressible as 


r= ary + ay(ry + ar°y) + agr’y’ + a,(ry’+2°y), 
wy bry | bory | bir yr bor y 7 b wy? 1 biry, ‘ (48) 
is bry t bry 1 bor yr boa yj bicy? | byry, 


; : : ; 4 4 
Ue=C, + C2(2° +") + e3ry? +e, (2+ y’). J 





we can construct the trial elements more easily (see 
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Moreover, for such a trial function ¢ the boundary conditions (46) are equivalent to 





ov, or, , Ors ~ 
T rs Q, (=I, l<y<.l), 
"Or O01 OY ; Y 
or, + Vo 0, (i a iz y: 1). r (49 
r,=0 y= 1) P 


where ¢=yn/(1—unz). 


Substituting (48) into (49), putting r=1 and equating the coefficients of y, y°, 7 to zero 
or putting s=y=1, we obtain the conditions for v given by (48) to satisfy (49) and hence to 


belong to D*. The results are the following six equalities: 





a(a,+3a,+ 5a,)+ (6,4 3b, ' 5b, )+ (4c, + 4¢,) =0, : 
a(@, + y+ d,) +- (6, 4 b; t b,) (), 
a(@y+3a4) + (b, +36.) t+ Se, 0), 
: > (50) 
a(y+ 3) ; (b. ' b. 0). 
al, bh, 0), 
( 2 ¢,+ 2¢,=0. 
J 
Bv inspection, we find the following 8 independent families, (1), (2), . (8), of 
parameters satisfying the system (50) of 6 equations in 14 unknowns: 
‘> |( 8 b C. 
(2) @=l, h, bh, oC. 
(3) bh, 1, 6é=—1, €:=1, Cs 1/2, 
(4) a | h “ 
(5) ¢,.=—¢ 1. ¢ ] 
(6 (ly | b,—b, a 
(7 h 1. 6=1 ( 1, ¢ l 
(S h 2.5 2 ( E, Gs 1/2 
The absent members of each family are equal to 0. 
In this wav, we have the following 8 elements r“ oo” ef. of r=1, 2, S 
contained in M@*—TD* 
; ry ory, a) | ) 
/ ry rey ory my) o(ry iy () 
/ 0 LY ary LY > LY ] (J y*) /2 
yp“ ry ory ory 0 
/ 0 0 0) (] r-)(] y-) 
/ ry rey a(1y ry) a(ry y) () 
/ 0 ry y y+ ry l—zr*y 
, 7® 0 (—sy' +r y 2(—xry+sr'y’) 1—(r*+y"*)/2 
2 Ther sills nent ng to zero the coefficients of n the first and second conditions of (49) coincides and gives tl { u) 





With these 8 elements as basis, we construct the trial function ¢ in accordance with (26): 


} 
r=Br+Br+ .. .+Be% (51) 
with free parameters 8,82, . . . ,Bs. 
Making use of the trial functions uw and ¢ given by (45) and (51) respectively, we get, after 
the necessary integrations, 


Tu? 2a Auciar (l—p) 5A jaa, (Aun=Ayi), 
u > | B ,.,,, (Ba= By), 

T*r? SS CuBib —p>IC7B Be,  (Ca=Ce), 
r2 >) Dubib (1—w>2D1,88; (Da=Dry), 

Tu,u)=(u,T*?) > 2p 8.=(1—w I Ene, 


with Ah, Bu, Ca, Dx, Fi given by the tables. 


A 
h = 
l 2 3 | ) 6 
l Ss 16 ba 16 16 8 
2 100e + 312 OH 11200 + 672 192 96 
: esse 
3 —, : sare 3600 + 96 2400 +816 288 
‘ ‘ 
| 16480 + 3336 1L008e + 336 240 
. 22248 37296 1680 17520 
> a” = 
( 8400 + 2640 
B 
l 2 3 } 5 6 
B.=0 (ik) 
B } S } s ba 4 
y 21 19 33 77 121 
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r 3 3 3 3 ) 7) ) 7 7 Ss 
k 3 5 7 8 5 7 8 7 8 8 For other values 
of i and Re 
Can =G. 
17$ 256 348 256 512 256 S704 1096 =17664 
Cu of 3 3 5 9 9 5 63 35 175 
D 
h . - 
l 2 3 } » t é e 
| 
' } 8 16 | a bs 32 32 
yg Tw Ty a5" 21 z 63° 7E° 
368 64 s 196 512 256 
2 —— hy —pt yj ‘t $y — _ — 
525 175 35 O45 O45 a2) 
144 14 128 32 992 2528 
] — — = i ~ 
30 175 15 135 1-9 25 
| S 32 32 
} T 0 T ao a 
19 1 225 245 
. 512 26 1024 
) = 0 . 
225 id 315 
H56 192 64 
F) T a a 
1617 53y 1+.9 
- 119744 LOSS 
‘ 17325 L759 
P 63424 
11025 
EB 
l l l l 2 2 
h 3 ° 7 s 7 & For other values of ¢ and 
| 32 ran 
, 6 32 ti4 32 64 o4 
k - > . 
] 9 9 ) 21 a) 
r l p 
Here o and + are constants given by ¢ a ; 


In the numerical computation we take «0.225, the value used by Ritz. Then the quan- 


titv 2e becomes an algebraic expression of the form 


Ag yi Yo wr’ bey - 
-€ 7 “/.)R.\'0 , =) 
r Br y Dy r’ Bar) y Dy 
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, 


where s, y are the column vectors with components | a;}, | Bk}, respectively, 2’, y’ are their 
transposed vectors, and A, B, C, D, EF are matrices with numerical elements. 

In order to minimize &, we choose the following procedure. First we minimize the first 
term on the right-hand side of (52). This gives rise to an eigenvalue problem of order 6, for 
which we determine the lowest eigenvalue m (which is the minimum value in question) and 
the associated eigenvector s=2,. Then we put s= 2, into (52); this converts (52) into a frac- 
tional form of which the numerator and denominator are known quadratic forms in y. We 
minimize this expression; this is equivalent to solving another eigenvalue problem of order 8 
and we determine the lowest eigenvalue m’ and the associated eigenvector y=y,. In the next 
step we should put y=y, into (52), converting (52) into a ratio of two quadratic forms in 2, 
and determine its minimum m’’. This procedure could be repeated indefinitely, and the 
sequence m’, m’’, . . . is certainly noninereasing. It is not clear whether these values con- 
verge to the true minimum of 2e. But this is not important, for actually we need not obtain 
this true minimum; it is sufficient for our purpose that a very small value of 2é is given by 
some set of r and Y. 

The actual computation was done by using SEAC, and it turned out that the value of 
m’ was sufficiently small, so that further computation was not necessary. ‘To avoid the uncon- 
trollable error arising from the use of the machine, the components of 2,,y, thus obtained have 


been rounded off to five significant figures, and the value of 2é€ for these arguments has been 


calculated anew on a desk calculator. The final values are 

ae 3.6100 000, 8 2.0362 000, 
ae 0.2071 000, B. 1.1853 000, 
Q3 0.0357 O00, B. 1.6269 000, 
a; 0.0335 000, 8, 0.3121 000, 
a 0.0012 000, 8. 0.5907 000, 
a 0.0017 O00, Bz 0.0112 000, 

8; 0.0176 000, 

By 0.1233 000, 
€ 0.009172, n 3.6266247. 


We can take a=0 and B=13.5 in (7) for Uo and L in conformity with (44). Thus we 


obtain 
(= 12.455433, L=12.430613., 


where the rounding-off errors are within 5 107° by the most conservative estimates. Thus , 
12.43805< Ay< 12.4555. 


For the upper bound we ean also use Up, given by (28). This gives 
PI Ry $ . 


U' p= 12.454620, 
which is somewhat sharper than 0°. In this wey we arrive et the final result: 
12.4306 < Ay< 12.4547. (53) 


If we adopt Ay= 12.4427, the average of upper and lower bounds, as the smallest positive 
eigenvalue of J//= 7*7, the relative error is less than 1/1000. 
The corresponding value of the fundamental frequency of a thin elastic square plate 
(u— 0.225) with free edges is 
\ \, Ie h wiuial Eh 
Vv 0.3420 ™ ’ 
2y aryl Ve \ pa \ pa 


with reletive error less than 1/2000. Here 2a= edge length, 24 thickness, p=density, and EF 


Young’s modulus. 
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A remark is necessary regarding the relation between the result of Ritz |2] and that of this 
study. Ritz gives the estimate \)S 12.43, and this appeared exceedingly sharp in view of the 
results shown here (53). Ritz used products of eigenfunctions of free bars as trial functions 
u”. This suggested that these trial functions might lead to a sharper upper bound than the 
one presented in this paper. Because the last figure in Ritz’ result is unreliable, the same 
calculation was carried out for purposes of comparison by retaining more significant figures. 

The result obtained gives only \y<12.488, a less sharp estimate than the one herein. This ] 
seems to show that, as far as the fundamental frequeney of a square plate is concerned, use of 
polynomial trial functions is preferable to products of eigenfunctions of free bars. 


4. Concluding Remarks 
4.1. This method is also applicable to a plate whose edge is not wholly free but subjected 


to conditions of a more general type. For instance, let us assume that C, the boundary of the 
plate under consideration, consists of three ares, (, Ci, and (>, and that the edge is free, sup- 


ported, and clamped along (, (,, and (, respectively. Among all angular points on C, we 
denote those on Cy by P;, Ps, ..., P, with s=s),8,...,8,. Then as in the case of the free plate 


discussed above, the vibration of this plate is characterized by the eigenvalue problem of a 


A . ~ ° A ° . . . 7 . A 
self-adjoint operator /7 in §. Again, // is the gradient operator of a quadratic functional J 
(uj, which is the restriction of J|u] to the class of functions w satisfying the boundary conditions 


0 =() on €- aa 
ou 54) 
0 on ¢’.. | 
on 


A A A 
Thus, following a course parallel to that of section 2, we can reduce // to 7* 7 and make theorem 


| applicable. The results of this decomposition are es follows. § and ’ being the same as 
A A 
before, Tu and 7*r are formally identical with Tu and 7*r, respectively, and now any wu be- 
A 
longing to the domain of 7 is subjected to the boundary conditions (54) and envy + belonging to 


A 
the domain of T* is subjected to the boundary conditions 


(V=0 on CC, +). 
d 
P+ j R=0 on (, 55 
«Is 
D)—0 =12 n 


where P, Q, R, D, are those defined in section 2. 

Applications of theorem 1 to the cases of such mixed boundary conditions will be treated 
elsewhere. 

In addition, if both ©, and ©, are empty, i. e., the whole edge is clamped, another .de- 
composition seemingly more convenient is possible. Namely, by means of (19) and the con- 
ditions u=0u/0n=0 on (, we have 


) drdy | ow drdy, 
D 


"Ou OU { ( Ou 
, u 
w/b Ou O07 * Ou Ou" 


,dardy 
Jpdor oy . 


and hence, in consideration of (10), we obtain 


Iu) Auydrdy. 


/ D 
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A 
Thus, taking 5='=12(D;1) we can decompose // into the desired form by using T, and 
T° such that 
Ou 


Tu =Au, (u=> 0 on), 
Cl 


T-v= Av. 


1.2. Before his derivation of theorem 1, one of the authors gave also the following 
theorem 2 [10], which stands in a close connection with the former. 

THeoreM 2. Let I be a self-adjoint operator of a Hilbert Space , and let (@,8) be an open 
interval containing at most a nondegenerate eigenvalue of Il but no other point of its spectrum. 
Let w be any vector of the domain Dy of HT with \w 1, and set n= (1w,w), é (11—n)w)? 

Hw ?—7. If &<(n—@) (B—n), then there is certainly an eigenvalue of H in (a,B) that satisfies 


the inequalities 


n <=)’ s7n-4 . (56) 


This theorem can be interpreted as follows. If we can construct a “trial function” w in 
accordance with the conditions stated above, we would obtain an estimate of the eigenvalue 
\’ in question by means of (56). Comparing theorem 1 with theorem 2, we note that, if a 
good choice of w is easy for the operator /7 under consideration, theorem 2 may be preferable 
inasmuch as it does not require the decomposition /7/=7*7, and also it involves fewer quan- 
tities to be calculated. However, this is not the case so far as (Pr.) is concerned, for the 
difficulties in obtaining a suitable w seem greater than those in obtaining a suitable pair of u 
and ¢ (essentially ¢) because of the increased complication of the boundary conditions of Dy. 
Furthermore, we have four components 7, 2, ?3, ? at our disposal in constructing the trial 
vector function ¢ belonging to T*, whereas we must make a single function w satisfy all boundary 
conditions in order that w@ Dy. For instance, let us consider the case of the square plate 


dealt with in section 3 and look for a polynomial trial function w contained in DyOM. For 
T 


any w belonging to ®, the boundary conditions (39) of Dy reduce to 


53 + (2 1) 0; gu=}, l<y 1). ] 
d t y 
ow, Ow 
arty (2=1, —1<y<), } (57) 
a yy 
o'w 0; (r=y 1). 
Oroy » 





On the other hand, a general form of a polynomial of degree 2N—2, belonging to I is given bv 
a z a) r Gr . 


Vo m+n 
- — ‘N 2m Lael? l j 27 1,,2m l a *) 
u A, (2 y J y ). (oS) 


Substituting (58) into (57), putting «= 1 and equating the coefficients of all powers of y to zero, 
or putting s=y=1, we derive the conditions for w given by (58) to satisfy (57). This gives 
rise to a system of 2N—3 equations in N?’/4 unknowns, or a system of 2N—4 equations in 
(V?—1)/4 unknowns, according as NV is even or odd. The smallest value of N permitting 
nontrivial solutions is seen to be 7. Therefore, in order to obtain a function of the type (58) 
in DyOoM, we have to solve at least a svstem of 10 equations in 12 unknowns, leaving two 
homogeneous parameters, 1. e., essentially only 1 free parameter in a 12-degree polynomial. 
1.3. Weinstein and Aronszajn have developed an approximating method [11],'* with a 
wide range of application, which also enables us to calculate upper and lower bounds of eigen- 


See |11] for a standard list of references concerning this method, 


185 








values of a certain kind of operators. Although the arguments used in establishing the con- 
vergence of their method are theoretically interesting, it appears that, from the viewpoint of 
practical applications, the essential parts of their method are the procedures leading from 
rough bounds to those that improve as sharply as possible. In this sense their method seems 
to have some similarity to the present one, and hence it will be worthwhile to reveal the mutual 
relations between these two methods and to compare their merits. However, we reserve 
this for some future occasion and remark here only that, with respect to lower bounds of 
eigenvalues in the case of the square plate treated in section 3, the eigenvalue problem of the 
operator 77 can play the role of the “auxiliary problem” in their method, of which the complete 
set of the exact eigenvalues and eigenvectors is required to be known. 
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A Thermochemical Study of the Reaction of Calcium 
Hydroxide, Silica Gel, and Water 


Edwin S. Newman 


The heat of reaction of calcium hydroxide, silica gel, and water was measured by a 


heat-of-solution method. 
solid ratio of about 1.0. 

Ca(OH),.+ SiO.-nH.O 
was found to be about 10.3 kilocalories. 


CaO-SiO,-nH,O+ mCa(OH).= (1 


is about 1.6 kilocalories per mole of added Ca(OH)». 
of reaction represents the heat of wetting of the reaction products. 


The reaction was carried out : 
The heat of reaction, 


t+ m)CaO-SiO,-(n-4 


t 25° C in pastes having a water 
AH at 25° C, for the reaction 


CaO-SiO.-nH,O 


The heat of the reaction 


m)H,O 


About half of the total measured heat 
Differential thermal 


analyses and drying experiments indicated that 3CaQ-2S8iO,-2H,O was formed in some of 


the pastes. 


l. Introduction 


The reactions of lime with silica and water are of 
fundamental importance in the study of the hydra- 
tion of portland cement, and are of particular interest 
in the study of the effect of pozzolans. The energy 
changes occurring dictate the direction and extent 
to which the reactions may proceed, and a knowl- 
edge of these changes will aid in elucidating the 
mechanism of the process. 

Measurements of the heat of reaction of lime with 
silica and water fall into three general classes: 
(1) Determination of the heat of hydration of port- 
land cement; (2) determination of the heat of hydra- 
tion of the individual compounds occurring in ce- 
ment; and (3) measurement of the heat of reaction 
between calcium hydroxide solutions and silica gel. 
Woods, Steinour, and Starke [1]! developed the 
heat-of-solution method of determining the heat of 
hydration of portland cement, which has evolved 
into an acceptance test [2,3] for use when it is im- 
portant to limit the heat evolved during the setting 
of concrete. Others [4,5,6] have measured the heat 
of hydration by direct methods, which are, however, 
limited to relatively short periods. The results of 
measurements falling in the first group are useful 
for control and have been extended to determining 
the contribution of the individual constituents of 
the cement to the heat of hydration. The direct 
measurements have been useful in explaining the 
important reactions that oecur in the early stages 
of the setting of concrete. 

The measurements of the heat of hydration of the 
individual compounds of portland cement have 
almost without exception been made by the heat-of- 
solution method. Lerch and Bogue [7], Brunauer, 
Haves, and Hass [8], Brisi [9], Thorvaldson, Brown, 
and Peaker [10], Cirilli [11], and others have deter- 
mined the heat of reaction of individual compounds 
with water. Many of these measurements were 
made with pastes of bydraulic material and water, 
approximating the proportions and consistency used 
in placing or testing portland cement. The heats of 





Figures in brackets indicate the literature references at the end of this paper 
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reaction obtained have been used to aid in the 
identification of the reaction products [8]. 

The third class of measurements, also made by the 
heat-of-solution method, have been made on the 
reaction products obtained [11, 12, 13] by the reac- 
tion between calcium hydroxide solution and silica 
gel performed in such a manner that a large quantity 
of solution came into contact with a small quantity 
of silica gel. These measurements also have been 
used to aid in identifying the reaction products 
obtained. 

The present paper is the report of the results of 
heat-of-solution measurements made at intervals over 
an extended period on pastes of calcium hydroxide, 
water, and silica gel. The constituents were finely 
divided and in intimate contact with each other and 
with the small amount of water, so that local deple- 
tion of the calcium hydroxide in solution might be 
avoided. These measurements should agree with 
those made in class three above, and together with 
those of both class two and class three, should pro- 
vide data for calculating the heats of formation of 
the reaction products, 


2. Apparatus, Materials, and Procedure 


2.1. Calorimetry 


The heat-of-solution method of measuring the heat 
of a reaction was used in this work. The procedure 
consists in measuring the heats of solution of the 
reactants and the products of the reaction under 
investigation. A solvent in which all the materials 
will dissolve should be chosen, or proper corrections 
made for the heat effects of any whim on substances. 
The final solution at the end of the heat-of-solution 
experiment made with the products must be identical 
in molar concentration and quantity with the final 
solution obtained at the end of the heat-of-solution 
experiment made with the reactants, or the heat 
effects of the reactions necessary to bring the two 
final solutions to identity must be obtained and used 
in the computations. Under these conditions, the 
heat of the reaction is the difference between the 








sums of the heats of solution of the products and the 
sums of the heats of solution of-equivalent amounts 
of the reactants in the same solvent. 

In the present work, 1- to 2-g samples of hydrated 
paste were dissolved in a mixture of nitric and 
hydrofluoric acids. A total weight of 600.0 ¢g of 
acid was used, of which 12.6 g (11.0 ml) was 48 
percent HF, the remainder being HNO,,26.03H,O 
(2.000 N at 25° [14]. The composition of the 
calorimeter acid charge was therefore 0.302HF,1.10 
HNO,,29.11H.O. Because some of the earlier work 
(9, 11, 12, 13] had been done with hydrochloric acid, 
some determinations were “yA using 600.0 ¢ of 
HC1,26.61H,O (2.000 N at 25° ) [14]. The heats 
of solution at 25° C were dete olive? in an isothermal- 
jacket, platinum-lined calorimeter which has been 
described previously [15]. At the end of determina- 
tions in which HCl was used, the calorimeter acid 
was filtered and the undissolved silica recovered and 
weighed. 

In determinations of the heat of solution using the 
nitric-hydrofluoric acid mixture, the quantity of 
sample dissolved was kept small enough to avoid the 
precipitation of CaF,. In this manner the necessity 
of correcting for the heat effect of the precipitation 
was avoided. If precipitation occurs, about 9.1 
keal will be evolved for each mole of CaF, so formed 
[16], and the heat of solution obtained will be too 
high. The data of Brunauer, Kantro, and Weise 
[17] showed a heat effect of about 60 cal/g of CaO 
(3.4 keal/mole) in excess of the normal value of the 
heat of solution of the CaO when the sample size 
was increased above about 0.37 g (in 420 g of mixed 
acids), where precipitation of CaF, was observed to 
begin. With calcium hydroxide in the present work 
with a larger calorimeter, the corresponding increase 
was about 20 cal/g of Ca(OH), (1.5 keal/mole) when 
the sample size was increased above about 1 g. It 
is evident that in both investigations only partial 
precipitation of CaF, occurred, as the full increase 
corresponding to 9.1 keal/mole was not obtained. 
It has been observed that when the CaO is combined 
in a silicate much larger samples can be dissolved 
without precipitating CaF,. With portland cement, 
for example, it is possible to dissolve completely 3-¢ 
samples in 425 g of mixed nitric-hydrofluoric acid 
[2] without pree =, of caleium fluoride, although 
there is approximately 1.5 g of CaO in the sample. 


2.2. Differential Thermal Analysis 


Differential-thermal-analysis curves were obtained 
with an apparatus that has been described previously 


[18]. The curves were recorded automatically, and 
the heating rate, controlled by a motor-driven 


variable transformer rotated at constant speed, was 
about 12.5 deg C/min. 


2.3. Preparation of the Pastes 


All materials were chemicals of analytical reagent 
quality. Single lots of calcium hydroxide and 
silicic acid were used in preparing all the pastes, in 
order to reduce possible variations in the results 


| years. 


caused by differences in surface or other properties of } 

the starting material. 
The calcium hydroxide and the silica gel were | 

heated over night at 1,200° C to obtain the losses on 


ignition. The residues were taken to be CaO and 
SiO,, respectively. Seventy-gram mixtures of the 


unignited materials were prepared. weighing the 

components to the nearest 0.01 g. The weight of | 
each material to be used was calculated from its 

ignition residue and the desired C/S ratio.” Fifteen | 
mixtures were prepared over the range of molar C/S 

ratios from 0.3 to 3.3 in such a way as to bracket and 

to include the simple ratios 0.5, 1.0, 1.5, 2.0, and 7. 

The weighed materials were tumbled for 2 hr j 

bottles with a few wooden balls to insure uniformity 
of composition. 


The individual mixtures were mixed with distilled | 


water and stored in wax-sealed glass vials in the man- 
ner described in Federal Specification SS—C—158¢ 
{2}. A water-solids ratio of about 1.0 was used in 
order to obtain pastes with sufficient fluidity to fill 
the vials easily. 


2.4. Testing the Hardened Pastes 


One vial of each mixture was opened for testing al 
the ages of 1, 7, and 28 days, i6 weeks, and 1 and 2% 


The samples were dried under vacuum for 
about 1 hr until enough water had been removed 


from the damp samples so that they could be crushed 
and passed through a No. 28 sieve without caking, 
The total residual water content of samples thus 
treated was about 50 percent. Crushing and sieving 
were done as expeditiously as possible to reduce 
reaction with CO, in the atmosphere. The samples 
were then placed in small rubber-stoppered Erlen- 
meyer flasks for ease in shaking. Portions were 
removed for heat-of-solution and loss-on-ignition 
measurements, CO, determinations, and chemical 
analysis. Only the first two measurements were 
made on every sample tested. The individual 
weighed portions of each paste were taken con- 
sec utively as quickly as possible, so that one loss-on- 
ignition value could be applied to all the samples of a 
single paste. The analytical procedures were started 
at once, and the heat-of-solution measurements were 
made thereafter, usually within 2 hr. The samples 
for the latter test were kept in weighing bottles until 
used, and no significant change of weight occurred in 
this period. 


3. Results and Discussion 


3.1. Heats of Solution of Calcium Hydroxide and 
Silica Gel 


Twelve determinations of the heat of solution of 
Ca(OH), in the HNO,-HF mixture were made using 
samples of 0.8 to 1.0 g. No precipitation of CaF, 
was Observed. Although the scatter was rather 
larger than desired, there appeared to be no correla- 


? The usual notation of C The 


molar ratio CaO/SiOn», for example, 


=Ca0O, , and H=H:20 will be used. 
will be written as C/s 


S=SiO2 
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tion with sample size over this short range. The 
individual values of the heats of solution were calcu- 
lated to a molar basis on the assumption that the 
ignition residue was CaQ. The average value, 
91.70+0.10 keal/mole,* was somewhat higher than 
that given by Brunauer, Hayes, and Hass [8] but 
agreed with the value 31.7 calculated from data 
available in NBS Circular 500 [16]. The experi- 
mental value 31.70 keal/mole was used for the heat 
of solution of calcium hydroxide in calculating the 
heat of reaction of calcium hydroxide and silica gel. 

Eight determinations of the heat of solution of the 
silica gel were made using samples of 1.0 to 2.0 g. 
The average value 34.78+0.13 keal/mole was ob- 
tained and used in calculations. Roth and Chall [21] 
obtained 34.49+0.02 keal/mole at 50° C for silica 
gel containing about 66.7 percent of water. Accord- 
ing to Mulert [22], the heat of solution of silica gel 
containing 10 percent of water is about 0.24 keal 
mole higher than that of gels containing large 
amounts of water. In calculations comparing heats 
of reaction obtained in this investigation with the 
results of others, correction has been made for this 
0.24 keal. 

After removal of the weighed portions for the 
above mentioned tests, the remainder of the crushed 
sample was returned to the evacuation chamber 
until the following day when differential thermal 
analyses were made and, in some cases, X-ray pat- 
terns were obtained. At 1 and 2% years the parts 
of the samples remaining after the usual first-day 
tests were evacuated through a dry-ice trap, accord- 
ing to the procedure of Copeland and Hayes [19], and 
the heat-of-solution measurements were repeated. 

In general, standard analytical procedures were 
used [2]. The sample for the determination of loss 
on ignition was dried for 2 hr at 110° C, placed in a 
cold furnace, and heated to 1,200° C over night. 
In some of the tests when the standard analysis was 
not made, the residue from the loss-on-ignition test 
was fumed twice with HF and heated to 850° C, and 
the residue weighed as CaF, [20]. From the data 
thus obtained, the C/S ratio of the sample was 
calculated. 

Some additional experiments were made to deter- 
mine the effect of sample size on the heats of solution 
of calcium hydroxide and silica gel, and also to 
determine the effects of either material on the heat 
of solution of the other. It was found that, within 
the precision of the present work, neither predis- 
solved Ca(OH), nor predissolved silica gel had any 
significant effect on the heat of solution of the other. 
Variation of the size of the sample between 0.5 and 
1.5 ¢ also had no significant effect on the heat of 
solution of either substance. These observations 
were true only if calcium fluoride was not precipi- 
tated. With larger samples, if the silica gel was 
dissolved before the caleium hydroxide was added to 
the solution, no precipitation of CaF, occurred, and 
the normal value was obtained for the heat of solu- 
tion of calcium hydroxide, even up to sample weights 


? The standard deviation of the average is given here and elsewhere in this 
paper, Heats unless otherwise indicated refer to ~AH at 25° C, 





of 1.5 g. However, if the calcium hydroxide was 
added first, a high value for its heat of solution was 
obtained, although seemingly insignificant amounts 
of precipitated CaF, were visible. If the silica gel 
was then added to the same solution in the calorim- 
eter, a low value was obtained for its heat of solution, 
and the calcium fluoride redissolved. In one experi- 
ment the sum of the high value obtained for the heat 
of solution of calcium hydroxide and the low value 
obtained for silica gel was equal, within 1.2 cal/g, to 
the sum of the normal value of the heat of solution 
of silica gel and the value obtained for the heat of 
solution of calcium hydroxide, when the presence of 
predissolved silica gel prevented the precipitation of 
calcium fluoride. 

The heat of solution of the calcium hydroxide in 
HCl, 26.61H,O (2.000 N at 25° C) was determined 
for use in calculations to be made using the heat of 
decomposition of the calcium hydroxide-silica gel 
pastesinthatacid. The heat effect of adding the silica 
gel was also determined. The values obtained were 
31.26 and 1.01 keal/mole for the calcium hydroxide 
and silica gel, respectively. Tests showed that the 
presence of either substance previously added to the 
calorimeter solution did not alter significantly the 
heat effect obtained for the other material. 


3.2. Heats of Solution of the Pastes 


The analyses and heats of solution of the pastes 
are given in table 1. The heats of solution were 
calculated to kilocalories per mole of silica in the 
paste. For this purpose, the percentage of silica in 
the paste, s, was calculated from the original molar 
(/S ratio and the percentage residue on ignition, 
R, by the equation s=1.071R/(C/S+1.071), in 
which 1.071 is the ratio of the molecular weights of 
SiO, and CaO, 60.06 and 56.08, respectively. 

The C/S ratios calculated from the analyses in 
table 1 were usually somewhat less than the propor- 
tions weighed in making the pastes. However, 
these calculated ratios are affected by any analytical 
errors and do not allow for MgO, Fe:O;, Al,O;, or 
other nonvolatile impurities in the original materials. 
Because the proportions of the original mixes and 
the molar heats of solution of the calcium hydroxide 
and silica gel were calculated on the basis of the 
ignited weights of the starting materials, the original 
(‘/S ratios were used in all computations. 


3.3. Heats of Reaction 
Based on 1 mole of silica, the equation representing 
the formation of a hydrated calcium silicate by the 
reaction of calcium hydroxide and silica gel in a 
paste may be written 
(C'/S)Ca(OH),+SiO,-pH,O + qH,O— 
r(aCaO-bSiO,-nH,O) + (C/S—aa)Ca(OH), 
+- (1—bx)SiO,-pH,O+ (q—naz+ pba+azx)H,0 (1) 


The heat of the reaction to form the hydrated 
calcium silicate may be found from the difference in 
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TABLE 1. Analyses and heats of solution of calcium hydroxide—silica gel pasles 
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112 (28. 30) (20.71 50.99 1.464 116.2 243.9 3. 55 | 
1.500 365 33.13 | 23.69 42.48 1.498 131.7 285, 0 3.90 
, 365P 50.65 | 36.77 2.07 11.54 1.475 180.3 
570 28.24 | 20.53 KO.90 1.473 242. 6 
570P 1.94 13.00 | 
846 1.38 | 44.46 273.8 | 
MOP 12.84 166.8 | 





1“P indicates that sample was kept at water vapor pressure of about 14 micron (ice at dry ice t mperature) until loss fell below 0.5 mg/g/day. 
2 Figures in parentheses were obtained by fuming the ignition-loss s imple with HF and weighing the CaF, formed, as described in the text 
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the sums of the heats of solution of the products 
and the reactants in the usual manner. If heats of 
dilution are neglected, the uncombined calcium 
hydroxide, silica gel, and water will contribute 
nothing to this difference, because their molar 
heats of solution will be the same after the reaction 
as before. Equation (1) can then be rewritten 


arCa(OH),+ b2SiO,-pH.O + (na —ar— 
pbrs)H,O 
2(aCaO-bSiO,-nH,O) + AH (2) 


where AH is the observed difference between the 
heat of solution of the paste containing the hydrated 
calcium silicate and the sum of the heats of solution 
of the reactants from which the paste was formed. 
If AH is calculated per mole of silica in the whole 
paste, its value increases with the amount of silica 
entering the reaction (i. e., 6z) until a final constant 
value is reached when all of the silica has been 
combined (62=1). This calculation does not depend 
on a knowledge of the composition of the hydrated 
calcium silicates formed. If the final values of AH 
are plotted against the C/S ratios of a series of 
pastes, the curve will consist of a series of straight 
lines intersecting at the compositions of the hydrated 
calcium silicates formed. The slope of each line 
will be proportional to the heat of the reaction 
involved. 

In the foregoing discussion a number of corrections 
have been neglected. These are: The heat effect of 
dissolving silica gel and calcium hydroxide separately 
rather than simultaneously, the effect of sample size 
on the value obtained for the heat of solution 
(excluding the effect of precipitated CaF,), and the 
heat effect of diluting the calorimeter acid with the 


| others. 


water formed in the hydration reaction or present 
in the paste. Experiments have already been 
described that showed the first two of these effects 
to be negligible in the present work. The heat of 
solution of the water contained in the paste is 
approximately equal to the partial molar heat content 
of water in HNO, 26-03H,O and in HF, 100H,0, 
(0.002 and 0.0004 keal/mole of H,O, respectively [16]. 
These quantities are also negligible. 

The heat-of-reaction technique was used by 
Cirilli [12] to determine the composition of the 
hydrated calcium silicate produced by the reaction 
between calcium hydroxide solution and silica gel. 
He determined the heat of decomposition of the 
reaction product in hydrochloric acid in calories per 
gram of CaO. The heat of solution calculated 
from these measurements remained constant up to 
C/S=1 and then increased with further increase 
in C/S. It was hoped that the procedure used in 
the present investigation would, by reason of the 
longer times of reaction and larger C/S ratios 
involved, give on the curve of AH versus C/S ratio 
a more definite change of slope indicating the forma- 
tion and composition of a compound containing 
more CaO than monocalcium silicate hydrate, such 
as the 3CaQ-2Si0O,-nH,O reported by Cirilli and 
Unfortunately this hope was not realized. 

In figure 1 are shown the heats of reaction of 
Ca(OH),, silica gel, and water (reaction 2) for each 
of the pastes of table 1, plotted as a function of the 
time. The major portion of the heat of reaction 
was released in 1 to 7 days. At C/S ratios larger 
than 0.5, the reaction proceeded further slowly. 
This behavior suggests that at any C/S ratio the 
monocalcium silicate hydrate reported by Cirilli was 
formed promptly, and that a more basic product 
was formed subsequently in pastes with C/S ratios 
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larger than 1. In figure 2 are shown the heats of 
reaction of the pastes at 1 year plotted against the 
C/S ratio. The two straight lines A and B shown in 
the figure, represent the least-squares equations 
calculated from the present data for C/S ratios, for 
A, of 0.3 to 1.0 inclusive, and, for B, of 1.0 to 3.33 
inclusive. The average value of the intersections 
of these two lines with C/S=1 is 10.3 keal/mole of 
silica. This value is taken as the heat of reaction 
of calcium hydroxide and silica gel to form the 
monocalcium silicate hydrate. The corresponding 
value calculated from Cirilli’s data [11,12,13] is 10.96 
keal/mole. The slope of the line between C/S=1 
and C/S=3.33 is 1.58. The value of 1.6 keal/mole 
of added Ca(OH), is therefore taken as the heat of 
reaction of calcium hydroxide and 1 mole of mono- 
calcium silicate at 1 year to form material of C/S 
ratios higher than one. This heat of reaction is so 
small that it easily might be heat of adsorption of 
calcium hydroxide on the monocalcium silicate 
hydrate. The precision of the data is not adequate 
to show definitely a change of slope indicating a 
change in the way the lime is combined, such as 
that reported by Cirilli at C/S=1.5. 

The data taken at 2% vears are similar to the 1 
year results. The heat of reaction and the slope of 
the line B were found to be 10.5 and 1.06 keal/mole, 
respectively. The slightly increased amounts of 
CO, found in the samples at the iater age may in 
part account for the differences in the heat effects. 
No correction has been made for the small CO, con- 
tent of any of the samples. It is more likely, how- 
ever, that surface or other changes in the pastes are 
responsible tor the differences. 

Brunauer, Hayes, and Hass [8] determined the 
heat of hydration of C,S and C,S to be 22.0 and 4.1 
keal/mole, respectively. From their results the 
heats of reaction of Ca(OH),, water, and silica gel 


15 


2 





5 
\ 


~AH, HEAT OF REACTION, KCGI / MOLE OF SiO 
wo 








0 1.0 2.0 3.0 


MOLAR RATIO C/s 
Figure 2.—Heat of reaction of calcium hydroride and silica ge] 
at 1 year. 


©, Newman; ©, Brunauer, Hayes, ard Hass [8]; @, Cirilli [11, 12, 13] 


' 


to form the hydrated calcium silicates produced by 
the hydration of C,8 and C,S have been calculated 
to be 11.96 and 12.13 keal/mole, respectively. These 
values are plotted in figure 2. Considering the un- 
certainties involved, the agreement is good, 

The heats of reaction calculated from Cirilli’s data 
are also plotted in figure 2. The agreement. is 
excellent until the C/S ratio becomes larger than 1,0, 
Cirilli’s results are based on the assumption that the 
silica gel contributed nothing to the heat of solution 
in HCl, being precipitated at the end of the heat-of- 
solution experiment in its original form. In the 
present work this assumption could not be made, 
because the heats of reaction as measured with 
nitric-hydrofluoric acid mixtures were not the same 
as those determined using hydrochloric acid. A 
considerable amount of the silica in the paste 
remained dissolved in the HCl at the end of each 
determination, and its heat of solution was included 
in the measured value. In figure 3 is shown the 
fraction of the silica dissolved in the HCl, plotted as 
a function of the C/S ratio of the paste. If the 
soluble silica were taken an indication of the 
reaction between calcium hydroxide, silica gel, and 
water, as was done by Shaw and Maclntire [23], it 
would appear that the reaction was not ended, even 
after the lapse of 2 years and in the presence of 
excess Ca(OH)... The measurement of the fraction 
of soluble silica does seem to represent some prop- 
erty of the paste, as it generally increases with in- 
creased age of the paste, and it would appear reason- 
able to correct the measured heat of reaction 
accordingly. However, if a correction is made, 
using the silica insoluble in HCl to represent unre- 
acted silica, large (at low C/S) and inconsistent 
values are obtained for the heat of reaction. The 
values shown in figures 1 and 2, therefore, have not 
been corrected. 

Thermochemical calculations show that the differ- 
ence in the heats of reaction obtained in the two 
solvents may be attributed to the difference in the 
heat content of the silica gel in its original state, and 
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as a component of the hydrochloric acid solutions 
resulting when the hydrated pastes were dissolved. 
The differences, after dividing by the fractions of the 
silica dissolved in the hydrochloric acid solutions, 
were scattered about an average value (excluding 
those determined at C/S ratios of less than 0.7) of 
AH=—4.8+0.25 keal/mole of dissolved silica. This | 
value may be considered as a measure of the heat | 
of gelation of silica gel from the low-molecular- 
weight form occurring in hydrocaloric acid solution, 
Tourky [24] has given AH 8.9 keal/mole as the 
heat of gelation from alkaline solution, and Green- 
berg [25] has estimated a value of —2.7 keal/mole 
for the heat of polymerization of silica gel. 


3.4. Heat of Adsorption of Water 


It has been concluded [8] that under the conditions 
of drying used in this work all of the water adsorbed 
on the hydrated calcium silicate in the paste would 
be removed. The difference in the heats of solution 
of the pastes before and after drying should, therefore, 
be equal to the heat of adsorption of the water. 
These heats of adsorption are shown in figure 4 as a 
function of the C/S ratio. The values found both 
at 1 vear and 2% years are given; the latter are 
generally the smaller. The heats of adsorption of 
water on the reactants were calculated from similar 
heats of solution determined before and after drying 
the calcium hydroxide and silica gel. The line repre- 
senting the reactants is shown in figure 4, calculated 
from the equation Hgio,+(C/S)Heao 
where Hgjo, is the heat of adsorption of water on 
the dried silica gel (of composition SiQ,,0.21H,O) 
and Heao is the heat of adsorption of water on the 
dried Ca(OH),.. The data indicate that the heat of 
adsorption is considerably larger for the pastes than 
for the reactants, and that it becomes slightly larger 
as the C/S ratio increases. Older samples show less 
heat of adsorption of water. These phenomena may 
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be caused by differences in surface area of the ma- 
terials. Brunauer et al. [8] have calculated the heat 
of adsorption of water on hydrated di- and tri- 
calcium silicate from their data combined with that 
of Verbeck and Foster [26] and found values corre- 
sponding to 4.8 keal/mole of SiO,. Their values, 
about 2 kcal/mole smaller than those found in the 
present work, are plotted in figure 4. Five values 
calculated from Cirilli’s data are also shown. 


3.5. Differential Thermal Analysis 


Differential thermal analyses (DTA) were obtained 
for each sample at each age of test. Partially dried 
samples were tested. Figures 5, 6, and 7 show some 
of the curves of particular interest. Curves obtained 
at 16 weeks are shown in figures 5 and 6. Figure 7 
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shows the curves obtained at six ages of test for the 
sample with a C/S ratio of 1.5. All the features 
observed for any of the samples are shown on one or 
another of these curves. The curves all show a 
broad endotherm in the range 80° to 350° C, which 
is attributed in part to the loss of adsorbed and free 
water and which resembles that found for silica gel. 
Two other endotherms found in the range 350° to 
500° C were assumed to be associated with the loss 
of chemically bound water from the hydrated calcium 
silicates formed in the pastes. 

The presence of calcium hydroxide is clearly shown 
in the curves in figure 6, representing samples with 
(/S ratios of 1.6 or more at the age of 16 weeks. 
The endotherm representing its decomposition in 
the DTA apparatus occurs between 480° and 600 
( with individual peak temperatures of 525° to 570 
(in the various samples. The increased amount of 
calcium hydroxide with increased C/S ratio is clearly 
indicated by the relative sizes of the endotherms. 
The continued reaction of Ca(OH), with silica gel, 
or with hydrated calcium silicate already formed 
with a C/S ratio less than 1.5, is clearly shown by the 
curves of figure 7, where the curves for a sample of 
that ratio are plotted for the different ages tested. 
The large amount of Ca(OH), present at 1 day is 
reduced by continued reaction, as indicated by the 
reduced size of the endotherm with increased age, 


until at 16 weeks no Ca(OH), was feund by DTA. 
Curves not presented showed no Ca(OH), at 1 day 
for samples having (/S ratios of 0.9 or less. DTA 
curves of all samples with C/S ratios greater than 
1.5 showed the presence of Ca(OH), at all ages 


tested. It is concluded from this part of the DTA 
curves that the calcium hydroxide reacts very 
rapidly with the silica gel up to a C/S ratio of about 


1, and that little further reaction occurs, within the 
time limit of these tests, after a silicate with the C/S 
ratio of 1.5 has been formed. 

Two exotherms representing reactions between 
CaO and SiO, in the DTA apparatus were found, 
one at about 830° C and one at about 900° C. The 
exact temperatures at which they are located on the 
DTA curves are dependent on many variables, 
including the degree of packing and the size of the 
samples, the placing of the thermocouples, and the 
fraction of the sample entering the reactions, so for 
convenience the exotherms will be identified merely 
by their relative positions at the higher or lower 
temperature. (The approximately linear time-tem- 
perature curves of the samples were violently 
distorted when these exotherms occurred, but the 
exotherms are plotted as though each time-tempera- 
ture curve had followed the straight line connecting 
its branches above and below the exotherms. This 
method of plotting is approximately equivalent to 
plotting the differential temperature against the 
temperature of the inert body over this portion of 
the DTA curve.) 

These two exotherms indicate that changes 
occurred in the pastes after the Ca(OH), had dis- 
appeared. For example, from the curves in figure 7, 
the exotherm at the higher temperature was the 





smaller at 1 day but grew, somewhat at the expense 
of the po hol as the age of the sample increased. 
The lower endotherms, howe ‘ver, never completely | 
disappeared from the curve for the sample with 4 


C/S ratio of 1.5. Similar behavior occurred at 
higher ratios, although the intensity of the upper 


exotherm was less than at 1.5. At lower C/S ratios 
except for the sample with a ratio of 1.2, there was 
practically no change of the two exotherms with | 
age of the pastes. There was, however, a change 
with change in composition. As the C/S ratio in- 
creased from 0.3 to 0.7 the lower exotherm increased 
at the expense of the upper until the latter dis. 
appeared. At ratios of 0.7 to 1.0 only the lower 
exotherm was found, but at higher ratios both were 

again present. Up to 4 weeks the sample with a 
C S ratio of 1.2 behaved like that with the 15 
ratio, but at later ages only the lower exotherm 
was present. 

The behavior of the exotherms indicates a change 
in the nature of the reacted paste between C/S ra- 
tios of 1.2 and 1.5. It must be emphasized that 
the exotherms represent reactions occurring in the 
furnace after the hydrated compounds have been 
destroyed, and they can not at present be related 
directly to the nature of the paste itself. X-ray 
patterns were made on samples obtained by inter- 
rupting DTA’s at temperatures below, between, and 
above the two exotherms, but the reactions causing 
the exotherms were not identified. Lines of CaO, 
of B and yC,.S, and @ and 8 CS were found in all the 
patterns. Kalousek [27] found similar exotherms 
at somewhat different temperatures for mixtures of 
CaO and SiO, that had undergone hydrothermal 
treatment, and also for a product of the composition 
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(.S,H, prepared at room temperature. Differences 
in apparatus and techniques could easily account for 
differences in his observed temperatures and those 
found here. 


3.6. Nonevaporable Water 


The nonevaporable water remaining in the dried 
pastes must be held in 1 of 3 ways. A portion is 
chemically bound in Ca(OH), in samples containing 
that material in excess either free or adsorbed. A 
portion is retained by the silica gel remaining in 
samples having a proportion of unreacted. silica. 
The remainder of the water is retained in the calcium 
silicate, and a study of the data in table 1 should 
indicate the amount of water so held in these pastes. 

Bernal [28] and Taylor [29] have given formulas 
corresponding to 3CaQ-2SiO,-3H,O and 3Ca0O-2Si0,- 
2H.O, respectively, for the composition of one of the 
hydrated calcium silicates formed at room tempera- 
tures. Accordingly, computations were made from 
the data in table 1 to see to which, if either, of these 
compositions the hydrated silicate formed in the 
pastes corresponded. The amounts of SiQ, and CaO 
in the dried samples were calculated from the original 
C/S ratios and the losses on ignition. The water in 
each paste was taken as the difference between the 
ignition loss and the CO, content. The total CaO 
was reduced by the amount assumed combined with 
the CO. to form CaCO,. It had been found that 
vacuum drying at room temperature would not re- 
duce the water content of the silica gel below an H/S 
ratio of 0.21, and it was assumed that this amount 
of water was retained by any unreacted silica in the 
pastes. It was assumed that in pastes of C/S ratios 
less than 1.00 the constituents would be SiO,, 0.21H,O 
and CSH,, and that in pastes of C/S ratios larger 
than 1.5 the constituents would be Ca(OH), and 
C,..5H,. The values of n were calculated on these 
assumptions from the CaO, SiOQ,, and H,O contents 
of the dried pastes. For the two samples with C/S 
ratios between 1.0 and 1.5, the values of n were cal- 
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culated on the assumption that C,,SH, was formed, 
where m is the corrected C/S ratio of the paste. 

The results of these calculations are shown in 
figure 8, where the values of n are plotted against the 
corrected C/S ratios. The value of n is less than the 


| C/S ratio at all values of C/S greater than 0.7, and 


is approximately 1 at C/S ratios greater than 1.5. 
The hydrated silicate formed in many of these pastes, 
therefore, appears to correspond in composition to 
Tavylor’s formula, 3CaQO-2Si0,-2H,0. 


4. Summary 


The heat of reaction of calcium hydroxide in pastes 
with silica gel and water was measured by a heat-of- 
solution method. It was found that the heat of the 
reaction 

Ca(OH), +Si0,.-nH,O =CaO-SiO,-nH,O 
is about 10.3 keal, of which about 5 keal is the heat 
of wetting of the reaction product. The heat of the 
reaction 


CaO-SiO,-7 HO + mCa(OH).= (1 
(n+-m)H,O 


+-m)CaQ-SiO0,- 


is about 1.6 keal per mole of added Ca(OH),. 

Calcium hydroxide reacts with silica gel and water 
in pastes at room temperature to form hydrated 
calcium silicate with a C/S ratio as high as 1.5. 
Above that ratio the Ca(OH), is less firmly bound so 
that its presence is revealed by differential thermal 
analysis, although the heat of reaction per mole of 
silica continues to increase. Continued changes in 
the paste with time and with C/S ratio are indicated 
by changes in the behavior of the heated pastes in 
the DTA furnace at temperatures above 800°C, 

Drying experiments indicated that the nonevapo- 
rable water remaining in the pastes was less than 
that equivalent to the CaO on a 1:1 molar basis. 
Above a C/S ratio of 1.5 the compositions of the dried 
pastes corresponded to the formula [3CaQ-2Si0,- 
2H.O} [Ca(OH Dole. 
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On Some Expansions for Bessel Integral Functions 
F. Oberhettinger' 


Investigations by van der Pol and Humbert concerning the Bessel integral function of 
order zero are extended to Bessel funetions of other kinds and to functions related to Bessel 


functions. 
The Bessel integral function of order zero, 


Jil) ft LJ (t)dt, (1) 
was introduced by van der Pol [8].2. Equivalent [3, 4, 6] with (1) is 
, Fao , 
Jio(x)= J Ci(x sin t)dt. (2) 
TJ0 


Various properties of this function have been investigated by Humbert [4]. A number of 
definite integrals, especially of the Fourier or Laplace transform type, can be reduced to (1). 
For this reason tables of Jio(#) have been published previously [5]. The numerical computa- 
tion of Jio(x) requires the transformation of the integral in (1) into an expression involving a 
series in either ascending powers (for not too large z) or descending powers (for large x) of the 
variable x. One of these expansions [4, 6] is 


(—1)"(2/2)?" 


. (3) 
2n(n!)? 


: we » 
Ji(r)=y+log gt 2 
- nm 
where one can write 
a (—1)"(a/2)*" 
pa 


—-— ss 2n(n!)? 


= [tse 1|dt. 


An asymptotic expansion of Jio(z) for large x [5, 7] is listed in (6). The preparation of a 
set of numerical tables by the Computation Laboratory made it necessary to derive expansions 
of types (3) and (6) for functions defined similarly to (1) but with Jo(t) replaced by another 
kind of Bessel or related function. The established results are listed below. (A list of nota- 
tions and auxiliary formulas is given at the end of this paper.) 


Sus: i & , l s\’ £ 
5V to() FI. t—'Y (t)dt if log 5) +y log 5 


2 — 2n 1 
+ =, 3. L)*(z/2) | vin+1)+3,—loe 5 (4) 


12 =) 2n(n!)? 


Ki(x) —| t~'Ky(t)dt —3 log ; )-7 log = sa. z +3: sd Ei n+1) + 5,—Hoe 5 | (5) 
Jz 2 


2 2 2 24° 47 2n(n!) 


Both these expressions correspond to (3) and are suitable for not too large x 
For large x we have the asymptotic expansions of the following functions: 


dt~ to eee IN te A ts in 


» - oz “4 ° < ° 
di t ) } al r) d J Jf L Jf Ir = 


(- ' Jt) aul r) 2 92.4 9?.42.6 J (x) | 92 92.4? 22.47.67 
t 
. ) 


Y (x) 


'‘ Guest worker, National Bureau of Standards, from The American University 
2 Figures in brackets indicate the literature references at the end of this paper 


197 











- ‘ - 2, 2°.4 , 2?.47.6 . 1 , 2? 2?.4°  2?.4?.6? 
f t-'Ko(t)dt~ -Ki(2) 3 +—+— +: | +a) +=+—-+—,-— +-:-]} (7) 
" ze" 2 x 2’ z' x r 
‘ea .  -. (2n)! P(2a)-*" 
t—"\H(t)—Y,(t)|dt~—>\(—1)4 —— (8 
| [M(t) > > n' OQn+1° 
“2 . 2 : 2 ae , (2n)!(2n—1)! ; 
J, [Ho(t)—Yo(t)] dt—_log (24)—™ ne ae -1) inl)? (2r)-?", (9) 
“2 2 2 2 = (2n)'(2n —1)! 
(t)—L,(t)| dt—=— log (22)——y~—— 3° : (Qa) ~2". (10 
J, Uo()—Lo(t)] d T °8 x! —p > (n!)? . 10) 
Furthermore, we derive the identity 
°@ ; l ee 
| t-?H(t)dt=(22)'H(x)-4 5 | t-'H(t)dt- (11) 
; 2), 


Note that with (6) we know the asymptotic expansion of | t~'H,(t)dt for large 2 by 
r 


virtue of (8). With this and the known asymptotic expansion of H,(x7) (see (21)), the asymp- 
totic expansion of (11) is likewise known. 
a. Proof of (4) and (5). We write 


%a *, 


| ! 'JAtidt= | ! 'JAt)dt— | t—.J.(t)dt, 


e Jt) 


and since (using term-by-term integration) 


( —1"(zr vA >- 2m 


fi: 'J,(t)dt=>D 


, 
a= m'T(v+m+1)(v+2m) 


we obtain by (23) 


lhe l (2/2) ; (—1)"(a2/2)t?™ 
tS (t)dt= |—— —>) = - . 
5 y r(l+yp) mai m'T(y-+m+1)(»+2m) 


This defines an integral function of v. If we now let »v tend to zero, using (13) we obtain the 
expansion (3). In order to prove (4), we differentiate the preceeding formula with respect to v. 
If we again let v tend to zero, using (12) and (13) we obtain (4). Finally, the relation 


- is 
K,(?) 5 ie 


Joi it) | iY ,(2t)], 





together with (3) and (4), yields (5). 


b. Proof of (6) and (7). We replace the Lommel functions in (24) by their asymptotic 
expansion (22), and the results in (6) follow immediately. Again, the relation 


» Bi ‘ rs 
K,(t) °S in| Jo(rt) ; iY ,(¢t)), 


together with (6), yields (7). (Note that the asymptotic expansion of 


Jt) 


| t“< Y(t) >dt 


can be obtained in the same way.) 


c. Proof of (8). From (21) we have, for v=0 and large ¢, 


4 254, ft T. . 
t (Hat) Yoo ~25(—1y2- | i aati 
Tn=0 


n! 
Integrating term by term, we obtain (8). 


d. Proof of (9). We use (19) and obtain, interchanging the order of integration, 


*, 9 ea 
| |H,(t)—Y,(t)|dt=— | v~(1+0?)-4(1—e-**)dv. 
0 us 


0 


Because it is necessary to represent the |. h. s. function of z as an integral of the Laplace 
transform type in order to apply Watson’s lemma [1] for the determination of its asymptotic 


expansion, we form, using (14), 
_ » 2 2° 2 
| (H,(t)—Y,(t)|dt—=— log r= | vo "[(1+0?)-4—« dot? [ e~**y—"[1—(1+0*)~4]de. 
0 T TIO TJ0 


The first integral on the r. h. s. is a constant equal to y+log 2. The second integral is of the 
Laplace transform type and can easily be developed asymptotically for large 2. We obtain 


o'[(1+e*)-#—e~"|do=lim | | v (14 rt) do | te-"do | 
vl 5b) F 3 
. ie 1+-(1-4 92) 16 
lim { i(—8)—log| ( : = | } 
6-0) D 3 


y+log 2, by (15). 


first 


The asymptotic expansion of the second integral for large z is as follows. 


We expand 
(—1)"(2n)12-2*y*"-* 


eo 'f1—(14+2°)72?|]=— : snmis 
; 2 (n!)? 


n=l 
The application of Watson’s lemma [1] gives (9). 


e. Proof of (10). We use (20) to obtain 


a ») *1 
[J(t)Lo(t)\dt =f v~(1—v*)-*2(1—e~™ dv, 
J 0 T 


and proceed precisely as before, arriving at the intermediate result 


. 9 9 9 1 
| [1,(t)—L,(t)|dt—= log (22) —= y=- { Ki r)+ J, eet 1a) } 
0 T 0 


We expand 


=] 
3 


v 1 (l- yp?) 4] = Si ——— ; 
n=1 (m?)° 
The application of Watson’s lemma [1] gives (10). (Note that the term —Ei(—z) 1s 
purely formal and may be omitted, because its asymptotic expansion consists of the product 
of an exponential function with negative exponent and a series in descending powers of z.) 
f. Proof of (11). By means of (18) we obtain 


%@ i) *1 
| t-°H(t)dt-= = (1—v*)’? si (vx)do. 
z 7 it] 
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By partia) integration, 


Ae - l , if ; . 
| t~*H,(t)dt=—s5 si (x)+5 27'A,(x)+ | v' sin~! p sin (vx)dze, 
rz & Wo 


but 


“1 : , l : l es 
| y~' sin! p sin (var)dv 5 WS (x)+5 | tH t)dt, 
{ z 


and (11) follows at once. 
Notations: 
J,(t) = Bessel function of the first kind. 
Y,(t)- - Bessel function of the second kind. 
1,(t)=modified Bessel function of the first kind. 
K,(t) modified Bessel function of the second kind. 
H(t) =Struve’s function. 
L,(t) ie */* A (it). 
All these functions are Bessel or related functions of the order v. 
I'(z)=Gamma function. 
¥(z)=T'’(z)/T(2). 
y= —wv(1)=0.577215 . . . . . . (Euler's constant). 


Su.v(z) = Lommel’s function. 


Od (2 T., 
| > | 5Yu(z) (9, p. 61, 62] (12) 
l T- =? 
1+ 2+( o— ) + XK 2° (13) 
rM1+z) , r~ 6/2 | ‘ 
(Taylor expansion of 1/f'(1+ 2) near z=0). 


Integral representations: 


log x | v~(e~*—e-* dp [2, p. 17] (14) 
0 
Ei(—é8) ~j v—'e-"dv (15) 
é 
Ci(6) ~{ v~! cos v dp (16) 
j 
si(6)- | v~' sin pv de. (17) 
5 


The last three functions are the exponential integral, the cosine integral, and the sine integral, 


respectively. 


on fi ve £ | _ l 
r(v+5 ALU) 2r ™5t) | (1—2’*)’-“@ sin (vt) de, Re »>—= (9, p. 328] (18) 


) 
2 ) 2 
. Fi 
H,(t)—Y,(t) | e~**(1+-n*)~""dp (9, p. 331] (19) 
TJo 
2 
I(t) —L,(t) | e~*"(1—v*)~— "de. 9, p. 338] (20) 
TJ 
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d isym ptotic expansions for large 2: 


ann te — = 
H,(z)—Y (z)~x-3” cos (wv)(5 7) ‘wif Ir ae n)(2n)! [9, p. 333] (21) 


n=0 


—])*—y? ( —])*—y?* —* = 
S, (z)~2 {i-% ) an u—1) le 3) ai_...\. (9, p. 351] (22) 





Integral formulas: 
oe or Pe _ 1 1 ' 3 ; ; 
| te" J, (t)dt= 2] ( su) /I (1 +57) — Re v< Re es (9, p. 391] (23) 


*2 J,(t) J (x) J (x) 
| g~ dt=2z S_» (2)—az S_1,0(2). 19, p. 350] (24) 
z Y(t) Y (x) Y (x) 


{1] Ek. T. Copson, Theory of a complex variable, p. 218 (Oxford Univ. Press, London, 1946). 

{2} A. Erdelyi, Higher transcendental functions, vol. 1, Bateman Manuscript Project, Calif. Inst. Technol. 
(MeGraw-Hill Book Co., Ine., New York, N. Y., 1953). 

(3) P. Humbert, Sur les fonetions de Bessel integrales, Compt. rend. 195, 854-855 (1932). 

[4] P. Humbert, Bessel integral functions, Proc. Edinburgh Math. Soc. 3, 276-285 (1933). 


[5] A. N. Lowan, G. Blanch, and M. Abramowitz, Table of Jio(x) J J,(t)/t dt and related functions, J. Math. 
I 
Physics 22, 57 (1948). 
[6] G. Petiau, La prt des fonetions de Bessel, p. 251-258 (Centre National de la Recherche Scientifique, 
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[7] V. G. Smith, An asymptotic expansion of Jig(z)=— J(t)/t dt, J. Math. Physies 22, 58-59 (1943). 
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[8] B. van der Pol, On the operational solution of linear differential equations and an investigation of the 
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WASHINGTON, February 27, 1957. 
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Table of the First Moment of Ranked Extremes 
Julius Lieblein* and Herbert E. Salzer** 


Let a sample in x independent random values from the extreme-value distribution, with 
c.d.f. &(y exp(—e-*), be arranged in decreasing order and denoted by y;, Yo, 
Ym,- - +», Un The table gives the expected values for all these order statistics for sample 
size not exceeding 25. For the larger samples, up to n=100, the expected values are given 
only for the first 26 largest values. 


1. Explanation and Use of Table _ However, there arise other situations in which it 
is important to know where these mean values /(y,,) 


are located, and with the great increase in applica- 
Kimball, gives the expected values, E(y,,), of order | tion of the extreme- value distribution this table 


statistics from the tvpe I distribution of ‘“‘reduced”’ should be a valuable asset. | 
largest values, Prob{ ¥ <y} =®(y)=exp(—e7"). This | A recent important use relates to the determina- 
distribution was discussed in considerable detail in a | tion of best linear unbiased estimators for the param- 
series of lectures by E. J. Gumbel [1].!. The entries | eters u, 8 of the distribution of largest values ? in 
were computed from the formula given by B. F. | general form 

Kimball {2}, 


Table 1, prepared at the suggestion of Bradford F. 


Prob {NX <r} =F (zr) =exp(—e7’), y=(xr—u)/B, 
. TE. . nt 

E(Y¥n)=C+ 2 | 1)'( P )a‘ In(n—t), or in general, for estimators of a given linear com- 
bination such as u+a8 with a known. The condi- 
tion of unbiasedness involves a knowledge of the 
means, /(y,,), of the order statistics in samples 
from the reduced distribution @(y). These are pro- 
vided by the present table. The second condition, 
— ; “best” 1. e., minimum variance subject to the first 
denotes the order statistics, that is, the observations | condition, involves in addition the variances and 
inasample of x arranged in decreasing order of magni- | eoyvariances of the order statistics (Ym), FT YmYp)- 
tude. The quantity C=0.57721566 . . . is Euler’s This problem has been treated by Lieblein [4], and 
constant. The values EY») are tabulated herein for numeric al results up to n 6 appear in a paper by 


Lieblein and Zelen [5]. 


where 


~ 
~ 
— 
~ 


YW2Y2o> ~~~ 2Ym2> ~~ » 2Yn-12Yr 


m= 1(1)min(n,26) 


n 1(1)10(5)60(10) LOO. 


In order to furnish values of /(y,,) for samples of 
sizes not included in the table an interpolation chart 
designed by B. F. Kimball is appended (fig. 1). 

The curves in the chart connect points having the 
same ranking. Thus, the lowest curve is for the 
expected values of the largest in a sample of size n; 
the second lowest curve relates to the second largest, 
ete. Each curve has two branches that meet at a 
vertex situated at the expected value of the median 
for an odd sample size. These vertices are connected 
by a nearly vertical line, and it is interesting to note 
how rapidly they approach the population median 


| 2. Interpolation Chart * for E(y,,,) 


The table of the first moment of ranked extremes 
is one of the many tables furnishing pertinent infor- 
mation concerning the extreme-value distribution of 
type I. At one time Kimball [2] advocated use of 
such values in plotting extreme-value data on a 
special-scale extreme-value plotting paper. Because 
in this method the /argest sample values are impor- 
tant, computations have not been carried beyond 
the 26th largest value (m—=26) for sample sizes n >26. 
Subsequent investigation [3] indicated other criteria 
that should be given consideration, such as the short- 
est confidence interval of fixed probability, which 


applies to the plotting position. In(—In 0.50) =0.3665 .. . as sample size n increases. 
*Present address, David T: ivylor Model Basin, W ashington, 1). 
**Present address, Convair-Astronautics, San Diego, Calif ? The scale parameter 8 is the same as the parameter 1/a used in reference [1]. 


' Figures in brackets indicate the literature references at the end of this paper. This section has been contributed by B. F. Kimball. 
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INTERPOLATION CHART FOR Ely.) 


where 


Fly) = exp (-e'") 


POSITION OF 
MEDIAN 





NUMBER IN SAMPLE 


Figure lL. 


The curves for m=constant rising to the left do not necessarily terminate in the positions shown in this 


figure However, values in the range beyond the 


curves as shown are rarely needed. In other words, the line joining these apparent end-points might better have been omitted 


For samples of n not in the table, the interpolation 
chart may be used mainly for finding the expected 
values for the ranks in the upper half of the sample. 
For example, let n=75. 
the horizontal line n=75 with the curves m=1, 
m=2,..., are read from the horizontal scale at the 
bottom as follows: 


n=75 
m E(Ym) 
l 4.90 
2 3.88 
3 3.38 
4 3.04 
26 O.SS 


Then the intersections of 


3. Method of Computation 


Logarithms of \ from 1 to 120 to 48 decimal 
places from J. Wolfram’s table [6] were punched on 
cards. The successive differences of the logarithms 
were then obtained with the aid of the IBM CPC-L. 
Sufficient differences were taken to attest to the 
accuracy of the functions and the individual differ- 
ences. Values of the binomial coefficients were 
obtained from Peters and Stein’s table {7|, and the 
computation of /(y,,) from eq (1) was then performed 
using desk calculators. The following partial check 
was then applied for n <25: 





lJ — 
S* Ei Yn)—C. 
IL m=1 ' 


For n>25 a differencing check was applied. The 
results are believed to be correct to within 2 units 
in the last significant figure. 


204 


Ra 


30 


oO 


the 


IN SAMPLE 


NUMBER 


Rank from 


top 


m 


19 


14 


P21 5th) 
i 
35 100 


0.8127 S817 


1080 S835 


Nido L2t 


5231 258 
4440 711 
S524 S26 
40065 935 


224 O42 


$153 099 


S459 S76 


S7US OOS 
S261 G5t 
2671 S822 
OSD SIS 
436 122 


0120 439 
2836 SUS 
5845 S81 
GSYUR 741 





2542 §57 


SORS 407 


6258 959 
1020 114 
2104 327 
O415 346 
S802 472 


7404 O44 


6192 453 
1963 746 
3791 O10 
2859 691 


1557 OOS 


O471 025 
0616 213 
Hag O55 
2812 545 


3068 504 


5195 330 
H536 OS!l 
8073 485 
YUS1 Sel 


JS82 SUS 


20) GAd 
2533 S28 


wo21 249 


3 QRA2 650 


— 


2503 728 
7132 872 
3403 534 


O47S 353 


SOIS STS 
S851 754 
872 S21 
2010 161 
0206 140 
1594 A79 
3458 O52 
$85 327 
7TSS3 747 


1326 S41 


UHLS HF 

4438 174 
23 291 
8237 41 

6044 Os 
$178 74 

2543 779 
1079 713 
74H 21 

S514 772 
734 333 
6278 773 
524 10 
$253 445 
3204 297 
2360) US1 
1443 471 
37 FA 
0364 O87 


1269 375 


2186 415 
3125 115 
Us LIZ 


9122 604 


224 14 

4 
1. 3838 782 
$ 3725 we 
2 SOLL O39 
2.5160 483 
2. 2540 897 


1. 6758 280 


) 4594 325 
0. 4036 136 
HN t) 


2. 3689 751 
1.2750 458 
0.6627 159 
+. ISS3 S853 

2545 345 


7772 937 


+ 
“ 
“4 
~ 
te 


3775 607 
-. 5304 220 
7022 321 
9119 529 
1, 2232 104 
n=35 


1325 637 


gag 

6030 713 
2.2544 241 
1. YSS7 OY 
1.7725 500 
l ed i2 
1.4293 167 
1, 2866 USS 
1.1574 147 
1.0386 471 
0, 9283 319 


6340 712 


5448 828 


1584 O15 
S755 427 
2942 &AS 
2146 659 


1362 310 
+ O585 513 
OLDSS 00S 
0963 110 
1744 5A7 


4.4802 387 
1791 033 
2 GOST G20 
2.6249 676 
2. S42 S68 


TABLE 


Rank from 
top 


m 


Who 


16 
17 
1s 
4 


20 
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+(-1)' () At In (n-4)., where C 


k 


Rat 


uler’s constant =0.57721566 


ik fron 


10 


V182 926 
7353 S848 
5779 Sil 
4379 596 
3114 804 


1956 921 
ORS5 487 
YSS5 O67 
8943 603 
SO51 S63 


7201 865 
6387 319 
5H02 S845 
4843 S15 
4106 182 


3386 325 
2680 442 


1986 51 

1301 405 

+. 0621 408 

OU55 967 
n=55 

$. 5845 489 

+. 5753 462 





3. 0659 996 





2. 7231 708 
2. 4635 211 
9¢ 

2 

1. 9249 664 
1. 7887 415 
1.6671 204 
1. 5564 180 
1.4546 060 
1. 3601 612 
1.2719 O91 
1. 1889 245 
1. 1104 647 
1. 0359 236 
0. 947 993 


S966 70S 
8311 S807 


7680 221 
7069 285 
476 664 
5900 293 
5338 321 


4789 079 
n=80 


4.9592 423 
3. 9529 307 
3.4465 700 
3. 1087 OR4 
» S502 897 


6437 091 
4703 881 
3208 O11 
IS8Y USS 
O709 G61 


1, 9640 281 
1. 8660 674 
7755 969 
6914 486 
6127 O44 


5386 292 
4686 252 
4021 905 
3389 408 
2785 021 


1, 2205 87S 
1. 1649 442 
1.1113 514 
0596 175 
1.0095 742 


0, 9610 726 


— me bo tS OS wmNwaw- 


O418 S41 
8627 537 
7071 623 
5691 446 
4447 183 


3310 904 


, 2262 243 


1285 871 
0359 939 
9505 O80 


8483 738 
7899 709 
7147 813 


.§423 661 


5723 476 


5043 959 
4382 186 
3735 525 
S101 567 
2475 O69 


. 1862 9O4 


n=60 


6715 602 
6631 331 
1545 855 
$125 S802 
5537 800 


3448 473 
1691 110 
O170 428 
R26 854 
7620 653 


6523 995 
5516 619 
4583 314 
3712 359 
2804 532 
2122 431 
1390 02° 
0692 342 
0025 197 
¥385 062 


S768 913 
8174 137 
TAOS 456 
7039 865 
6496 584 





l 


7651 672 


5926 280 
4438 397 
3128 486 
1956 S50 


0895 700 
9924 SIS 
9029 O38 
S106 690 
7418 600 


6687 427 
5997 200 
5343 001 
4720 725 
$126 O15 
3558 625 
3012 329 
2488 842 
1983 258 
1494 908 


1022 318 


2 
l 
l 
1 
] 


1 


1 
1 
l 
1 


1534 009 
9756 620 
8215 109 
6849 953 
5621 370 


4501 482 
3469 974 
2511 576 
1614 507 
0769 470 


0. 9968 991 


4 
3 
3 


) 
9 


bo to te bo 


1 
1 


l 
1 
1 
l 


l 
l 
] 
l 
] 


1 
1 
0 


4 
3 
3 
3 


9 
9 
2 
1 
l 


l 
1 
] 
l 
l 
1 
1 
1 
l 
! 


1 


9206 958 
78 204 
78 727 


_7104 609 






. §452 792 
5820 524 
. 5205 370 
4405 148 
.4017 883 


3441 759 
n=70 


8257 109 
8184 993 
3111 GOA 
9704 763 
7129 934 


5054 151 
3310 720 
1804 376 
0475 566 
9284 576 


8203 598 
7212 304 
6295 780 
5442 063 
4642 047 


3888 367 
3175 026 
2497 O78 
1850 392 
1231 476 


0637 3546 
0065 469 
9513 592 
8979 781 
. 8462 324 
7959 700 
n= 100 
1823 859 
1773 523 
6722 760 


3338 229 
O786 589 








8734 497 
7015 279 
5533 687 
4230 189 
3065 092 


2010 610 
1046 529 
O157 689 
9332 423 
8561 563 
7837 772 
7155 O86 
6508 589 
5804 186 
5308 422 


4748 362 
4211 483 
3695 608 
3198S S3Y 
2719 514 


2255 168 


Rank from 
top 


m 


to to to ty te 
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The calculations were done in the Computation 
Laboratory of the National Bureau of Standards. 
Computations were performed under the supervision 
of Herbert E. Salzer. The final checking was carried 
out under the direction of Lrene Stegun. 
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An Approximate Expression for Gamma-Ray Degradation 
Spectra 


U. Fano and Ann T. Nelms 


The equation for gamma-ray degradation spectra is reduced to a form suitable for a 


simple approximation. 


The approximation is valid in a wavelength range between two 


Compton wavelengths larger than the wavelength of the source and the region where photo- 


electric absorption becomes very intense; this range is wide for low Z materials. 
are provided for application to various materials and source energies. 


Input data 
The results of sample 


calculations are shown and compared with those obtained by standard numerical calculation. 


In the course of gamma-ray penetration, the spec- | 
trum of multiply scattered, degraded, secondaries 
approaches a shape that depends only on the prop- 
erties of the medium.' This limiting spectrum is 
that which when a gamma-ray source 
distributed uniformly throughout a medium and the 
only relevant factor is the accumulation at each point 
of secondaries of various energies. When degrada- 
tion spectra were first calculated numerically ,? M. H. 
Johnson pointed out, in a private discussion, that 
the main features of the degradation spectrum at low 
energies are determined by very simple arguments. 
A formalism of degradation theory, starting from 
this remark and from the corresponding neutron 
theory,’ was developed later.“ But this theory was 
not applied in actual gamma-ray calculations, 
because these were being carried out routinely by 
direct numerical computation.’ Lately, the exten- 
sion of calculations to degraded photon energies ~ 20 
kev has made it desirable to try to evaluate the simple 
approximate analytical formulas and compare the 
results with those obtained by numerical integration 
of the basic equations. These low photon energies 
are of interest, e. g., in LiH. This paper contains a 
description of the procedure, beginning with the 
qualitative analysis, together with a discussion of the 
results of the comparison. 

If successive steps of degradation by Compton 
scattering, characterized by wavelength — shifts 
t=\—X’, have approximately equal statistical dis- 
tribution (see footnote 4) with a mean (£), the prob- 
ability for each photon to attain a wavelength be- 
tween A and A+ Ad at the end of one of the steps of 
degradation is Ad/(é).° If the photon does “land” 
in this interval, it will ‘“‘spend” there a mean free 
path T/po(A), where poo is the probability of 
scattering per unit pathlength and coincides with the 
total u(A) when absorption is disregarded. For a 
source strength of S photons per unit volume and 


arises is 


(r 
= 





See, e. g., U. Fano, Nuceonics 11, Nos. 8 and 9, 55 (1953). 
P. R. Karr and J. C. Lamkin, Phys. Rev. 76, 1843 (1949). 
See, e. g.. R. Marshak, Rev. Mod. Phys. 19, 185 (1947), 
Placzek, Phys. Rev. 69, 423 (1946 

*U. Fano, Phys. Rev. 92, 328 (1953 

‘ Calculations of penetration by the moment method, as given by H. Goldstein 

ind J. E. Wilkins, Jr... NYO Report 3075 (1954), include, under the name of 


ind in particular, 
{ 


“‘zero-th moment”, the spectrum discussed here 


6 This problem is similar to the determinatien of the probability that a pedes- 
trian will step on a line drawn across a sidewalk. The probability equals the 
length of the pedestrian’s foot divided by his mean step length. 


unit time, we thus expect at each point a flux of 
photons in the assigned wavelength range 

I(d)AA=S{1/p(A)](Ad/(8)) (1) 
per unit area and per unit time. This approximate 
spectral distribution is expressed in the scale of 
wavelengths because it is in this scale that the 
statistical distribution of successive steps of degrada- 
tion is more nearly a constant. 

Equation (1) is subject to three main limitations: 
(a) It does not hold at photon energies > 500 kev, 
where the distribution of degradation steps varies 
rapidly from step to step; (b) it represents the 
spectrum for a “steady state’ of degration, which 
does not hold for photons that have seen scattered 
only a very few times because their spectral distri- 
bution depends on the source wavelength (‘‘tran- 
sient’ effect) (see footnote 4); (c) it does not take 
into account outright photon absorption, by photo- 
electric effect. The effect of absorption can easily 
be included in (1) as a progressive decay of intensity 
as \ increases, provided the probability of absorption 
at each step is small. 

Equation (1) is contained as a special case in the 
solution of degradation problems by successive ap- 
proximations (see footnote 4). In this solution the 
degradation spectrum is resolved into steady state 
and transient components. A numerical calculation 
of the first few terms of this solution for the X-ray 
degradation would be possible, but it was found 
more convenient to proceed as in the analogous 
electron degradation problem.’ One starts from an 
integral form of the degradation equation, namely 


«| A’ +2 
| dn’ T)| wave!) + f kw’ "yd" |=8 (2) 
r nN 


0 


which represents a sort of “balance sheet”’ of photons 

in the entire wavelength range from > to A. The 

source injects S photons per unit volume and unit 
ma 


ho. Of these, | dd’T(A’)mane(A’) are lost 


Xo 
through outright absorption, and the remainder by 
Compton scattering with final wavelength ’’>), 


time at 


’? L. V. Spencer and U. Fano, Phys. Rev. 93, 1172 (1954). 
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where k(X\’,\’")dd”" is the probability per unit path- 
length that a photon of wavelength ’ be scattered 
to a final wavelength between \’’ and \’’+dnX’’. 
Note that the wavelengths are expressed in Compton 
units of h/me=2.4X10~" em. 

Equation (2) should now be transformed so that 
it displays an approximate solution similar to (1) 
and so that the necessary corrections can be easily 
estimated and then calculated either by successive 
approximation or by numerical solution of the new 
equation. To minimize the difficulties (a) and (b) 
indicated above, we restrict our discussion to wave- 
lengths A> A9+-2; thus no photon from the source can 
reach \ in a single Compton scattering, because 
k(X’,d’’) =0 for \’’>)X’4+-2. As a first step we trans- 
form the Compton seattering term in (2), by setting 
\’’=)’+E and \’=A—?, and find 


| “aN TON’) | 


“y ei 


*\" 


kN) an” 


| dn dO n) | kO n,A-—n+ dé 
e ” 


ftr 
~~ 


Jae] ‘dln n)k(\—7,A— 9-4 


Because the procedure that is being developed is 
expected to apply when J and & vary slowly with », 
it Is convenient to represent (3) in a form suited to 
expansion in powers of 7. This is done by an 
operator symbol based on the power expansion 
formula 

>? 


f(\—n) FO) n~ fO)+5 Y Sef 


J 


e >F(r). (4) 


The integrand on the right of (3) can therefore be 
indicated in the form 

x 

T(A\— n)k(A— 9, AN— +) =e AT (NX) ROA, A+8). (5) 

We consider in particular the value ot (3) for 7 set 

equal to zero in the exponent of (5). This value 


Is 


10) | dé th(h, X42) =10)M,(0). 6 


The function 4,(A) defined here is the first moment 
of the wavelength shift and may be expressed as 
Useatt(=>, Which coincides with the denominator of (1). 
We add the right side of (6) to the expression (3) 
and subtract its left side from (3), utilize (5) and 
combine the result with the degradation equation 
(2), which becomes 


JOM.) + | dN’ T(X! ) tare (d’) 


+{ ds| os "A ] TOyka. A+-8)=S. (7 


The solution (1) is thereby obtained if the integrals 
in this equation can be disregarded by assuming 
Mans ~ 0, and 70/OA~0 in the exponent. 

Consider now the absorption, represented by the 
first integral in (7), taking into account that JM, 
should be effectively constant in the absence of 
absorption. The absorption term is then written 


\ *’ 
| dN par WLOX’) = | AN BIN) [TOODAMGO)], (8) 
Jy 


A 


showing that JAZ, is effectively reduced, as X in- 
| creases, by a fraction 


Mabs(A) 
MA) 


Maps! \) AX 


A\= - 
Mseatt(A) (f 


B(A) AX (9) 


| in every small interval AX. Notice that this fraction 
| has the dimensions of a pure number. This absorp- 
tion effect, taken by itself, would cause JM], to 


. ° "dr P 
| decay, as \ increases, according to exp - | B(r yan’ }. 
A . JA 


| Accordingly we replace the spectral distribution 
T(\) with a new unknown w(A)=1 defined by the 

expression 

| x. 

v Bir’) dn’ 


TX) Ma) J wi(X). (10) 


| When this expression for J(A) is substituted in (7), 
the source strength S drops out as a common factor. 
| In the first term we use the identity 


| 

*) Pp’ 

Bin’) dy’ *’ Bi") dd 

| e 2 1— | AN’ BIN’ Je 2 » (11) 


Ao 


| so that the second term from this expression may be 

combined with the absorption term in (7). In the 

| second integral of (7) the term exp [ -f"80) aN’ ] 
J, 


| may be factored out utilizing the operator formula 


after factoring, Mquation 


| (7) becomes now 


*A Bid") dr’ *¢ *¢ 
w (X) + ANB’) ¢ J { wn’) — | dé dn 
” , — BOA *( in © } P "€ 
Ls J oo ih wo b+{ ds | dn 


ALA) 
” Bir K(A,A-+-&) . » 
(. [. ] 1) MA) u dr) :. (13) 


(11) is applied again. 


This equation is fully equivalent to the initial 
eq (2), for A>A+2. Its zero-order approximation 
solution is w(A)=1, obtained by disregarding all the 
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integrals. To obtain a better approximation, the 
integrals may be evaluated by whatever procedure 
appears expedient. We proceed here by expanding 
into powers of the operator 0/O\ and of the rate of 
absorption B, to first order only. The integral 
over \’ vields no contribution to first order, because 
of the factor B(X’) and because the factor contained 
in the curly brackets vanishes in zero-th order. 
Where the zero order approximation 0 is poor, 
for AX Ao + 2, BCX’) is very small for most materials. 
The second integral is evaluated to first order by 


writ ing 


[oan( Ce) fan [2-009] 


bo © 
—5F | 5x8 | (14) 
and w(\)=1. Integration over & vields the first 
order estimate of the second integral 
ra) 1 AZ,(\) 
a Peni 5 2 f 
| ox Pp |aa¢0) ss 
where .V7,(\) dé &k(\,A+-£) indicates the second 


moment of the wavelength shift & Thus we find 


oO 1 4Z,()) 1 M,(\) 
IT, 2 Ma) 2 2 May? co 0 (16) 


, Di as 
10) S fie an 1 Oo 1.\M,(A) 
OA 2.M(A) 


1 M,(\) = 
8) sap tee} (17) 


The function (1/2).M5/.14, was determined analyti- 
cally using the Klein-Nishina formula for Compton 
scattering. A graphical differentiation was used to 
obtain (0/OX)(1/2)\2/\4,. These data, which are 
the same for all materials, are given in figure 1, 


. mr , 
The absorption exponent dd’ BX’), where only 
0 


photoelectric absorption is included in the tabulations 
of B(\’), is given for a few materials in figure 2. 
The computation of 8(\’) utilized total attenuation 
coefficients from G. Grodstein® and the total 
Compton cross section from the Klein-Nishina 
formu'a. The integration was carried out numeri- 
cally. Figure 3 gives values of ,7'(A). The solid 
curve is a dimensionless plot of the quantity 
[(ucompton/MTnompson)(¢)|~', Which is independent of 
atomic number. For a particular material, the 
nomograph on the left permits the conversion to 
units of g/em*®. If the source yields S photons 
g~'see!, S/M,(\) is in units em~*see~'(h/me)~*. 


With the data from figures 1 to 3 the zero-th order 


expression may be calculated, and the importance 
of the first order correction may be evaluated 





§G. W. Grodstein NBS Cireular 583 (1957 
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Figure 3. Spectral distribution of photon track length. 
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FIGurReE 4 


, Approximation, eq (17 ,A, exact numerical integratior 


Degradation spectra for LiH and H.O 


readily. This has been done for two cases, namely 
gamma rays with initial wavelength \y=1 in Lil 
and H,O, the source strength S being unity. The 
results are the curves in figure 4. The points in the 
figure show the results of direct numerical integration. 
The agreement is well within the expected accuracy 
of the approximation. The departures at increasing 
wavelength reflect the increasing strength of photo- 
absorption, which leads eventually to a breakdown 
of the validity of the assumption that the absorption 
changes little over a wavelength range (£) ~ 1. 

Of general interest is the absorption ‘‘cut-off”’ of 
the spectra. This point, which may be characterized 


— . N 
by a parameter dz satisfying the relation 4B(X") 
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Figure 5 Absorption "7 it-off”’ of the spectral distribution, 


d\’=1. It is the wavelength or energy at which 
the absorption reduces sharply the spectral distribu- 
tion. Figure 5 gives a curve of Az versus the atomic 
number Z. Again, the values for B(\’) take into 
account only photoelectric absorption. 


The authors thank L. V. Spencer and M. Berger 
for many helpful discussions, and J.C. Lamkin, who 
contributed the results of the direct numerical 
integration. 


WaAsHINGTON, March 12, 1957. 
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Some Infrared Bands of Deuterium Sulfide’ 
Harry C. Allen, Jr., Earle K. Plyler, and L. R. Blaine 


Two regions of the D.S spectrum have been observed, 5 and 3.6 x. 


An analysis of the 


absorption in the 3.6 w region attributes it to two overlapping bands with excited states 


(4 ,02,2 (1,1,0) and (0,1,1). 


The band center of the 


(O,1,1) band is found to be some 


50 em”! lower than the previously accepted value while the (1,1,0) band center is found 


some 60 em~! higher than the previous determination. 
arise from the (0,0,1) band centered near 1906 em~!, 
makes a detailed analysis impossible at this time. 

The excited state constants for the two bands near 3.6 uw are 


lower than previously reported. 


The absorption near 5 yw seems to 
Overlapping atmospheric absorption 
This band center is almost 100 ecm~! 


1,1,0) »9 2742.77 A*=5.478 B*=4.512 C*=2.395 
O,1,1) » 2754.44 A*=5.531 B*=4.521 C*=2.405. 
1. Introduction 3. Analysis 


Although the spectrum of hydrogen sulfide is now 
fairly well understood, very little work has been 
done on the totally deuterated compound. The 
spectrum of D.S was investigated with low resolu- 
tion by Bailey, Thompson, and Hale [1]? and with 
somewhat higher resolution by Nielsen and Niel- 
sen [2]. In general these two groups agreed as to 
the regions of absorption, the only real discrepancy 
arising in the region near 3.6 uw. Although Nielsen 
and Nielsen resolved the rotational structure in 
several regions they did not present a rotational 
analysis of any of their data. The vibrational 
assignments of the two groups were completely 
different. With this considerable confusion existing 
it was thought advisable to reinvestigate the spec- 
trum of D,S. 

One D,S band has been successfully analyzed [3} 
previously, and in the present work the analysis of 
the rotational structure in the region of 2750 em™! 
is reported. The fundamental region near 1900 em 
was also investigated although the atmospheric ab- 
sorption in this region makes a complete rotational 
analysis quite impossible at this time. 


l 


2. Experimental Procedure 


The intensity of absorption of D,S is weak, thus 
it was necessary to use several centimeters pressure 
of the gas and a path length of 10 m to observe the 
absorption in the 3.6 uw region. The absorption at 
5 w could be observed satisfactorily with a 1 m cell 
and several centimeters pressure of the gas. 

The spectra were recorded with the grating instru- 
ments of the Radiometry Section, using a 7,500 
lines/in. grating with a ruled surface of about 6 
by 8 in. and a cooled PbTe cell as a detector. The 
resolution was about 0.2 em~! The wavelengths of 
the lines were measured by employing a Fabry-Perot 
interferometer as previously described [4]. 


The work 
Commission 
? Figures in brackets indicate the literature references at the end of this paper 


reported herein was supported by the U. 8S. Atomic Energy 


Figure 1 is a reproduction of a recorder trace of 
the absorption near 2750 cm™'. The absorption 
shows a well-defined P-, Q-, R-branch type of struc- 
ture. The appearance of a collected Q-branch indi- 
cates that the change in electric moment is along 
the least inertial axis, hence A-type selection rules 
must be invoked [5]. 

Assignments of the strong lines in the P- and 
R-branches of the band could be made by using the 
previously determined ground state [3] of D.S. From 
these assignments estimates were made of the excited 
state inertial parameters from which a trial spectrum 
was calculated using published tables of the reduced 
energy |6], E(K). 

Relative intensities of the transitions were calcu- 
lated by combining the Boltzman factor and ap- 
propriate nuclear spin statistics with the published 
values of the line strength [5]. Line strength values 
were used for K=0.5. These were found to be close 
enough for our purposes. This trial spectrum was 
compared with the observed spectrum and _adjust- 
ments were made to improve the agreement by 
means of the partial derivatives of the energy with 
respect to the various inertial constants [3]. The 
final adjustment of the parameters was obtained by 
carrying out a least squares treatment [3] in which 
some 90 assignments were used. The agreement 
between the observed absorption and calculated 
transitions for these 90 assignments was excellent. 
The average deviation between the observed and 
calculated values was 0.07 em™. 

After assigning the strong transitions of the A-type 
band, it was found that there was a sizable number 
of strong absorption peaks still unaccounted for. 
Closer examination revealed that these peaks arose 
from a second transition obeying the selection rules 
for a B-type band. Proceeding in the same manner 
as for the A-type band it was possible to assign 
unambiguously about 70 transitions to these ab- 
sorption peaks with an average deviation of 0.06 
em~'. Strong transitions of these two bands ac- 
count for almost all the absorption in this region 
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Comparison of observed and calculated spectrum of DoS near 3.6 p 


The numbers on the peaks correspond to the numbers in column 1 of table 1 


The remaining peaks are due to a piling up of weaker 
transitions or to absorption by HDS which was 
present to a few percent in the sample. Because 
the moments of inertia in D,S are relatively large, 
the resulting line spacing in the vibrational rotational 
bands is quite close, resulting in severe overlapping of 
the absorptions due to neighboring transitions both 


Figure 1 compares the observed absorption with 
the two calculated bands whereas table 1 compares 
the observed and calculated frequencies of some of 
the stronger transitions. 

The excited state constants leading to the fit shown 
are 


in the same band and in two closely lying bands. | 4-type »—=2754.44 em™! A*t=5.531 B*=4.521 
This situation leads to some difficulty in making an C*—2 403 

analysis for it cuts down the number of assignments | piyye y»,—274977 em-! A*—5478 B*—4.512 
that can be made, thus introducing more uncertainty NNT li coitinteaie (* —9 295 aT oe 
into the determination of the inertial constants — 

than one would like. The strongest transitions in | -p),, previously determined ground [3] state was 


the P- and R-branches of both bands arise from 
levels that are strongly dependent on C*. These 
levels could be readily identified and assigned, 
hence the value of C* for each band should be very 
good. The strong Q-branch transitions of the A- 
type band arise from levels strongly dependent on 
A*. Although the calculated @Q-branch reproduces 
the observed absorption in remarkable detail, indi- 
cating a good value for A*, this may be somewhat 
fortuitous because many transitions of the B-type 
band are calculated to fall in this same region. The 
transitions determining 2* most directly in both 
bands as well as A* in the B-type band are spread 
throughout the band and in general are weaker 
transitions. As a result these inertial constants are 
probably somewhat more uncertain than the others 
although the general excellence of the fit for the 
transitions in the two bands that could be unam- 
biguously assigned indicates that the constants 
determined should be fairly reliable. 


found to be adequate for both bands. 

The absorption in the 5 uw region is too complicated 
by atmospheric absorption to enable a quantitative 
analysis. However, the D,S absorption seems to 
have a Q-branch located near 1906 cm! indicating 
an A-type band. The main transitions of the P- 
and R-branches that are free from contamination 
by water vapor absorption indicate that the center 
of the band is 1906 +1 em™!. 


4. Vibrational Assignment 


The A-type band at 2754.44 cm™! must have 2; 
odd. The most likely transition with +; odd that 
could occur in this frequency region is the transition 
with the upper state (7,72,73;)=(0,1,1). The B-type 
band must have v; even and hence the most likely 
assignment on the basis of frequency is for the 
excited state to be (1.1.0). 
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} TABLE | Comparison of observed and calculated spectra * 
—_ 
Calculated Caleulated 
Line Observed rransition Intensity Line Observed rransition Intensity 
calculated calculated 
,1,1 1,1,0 0,1,1 1,10 
cm cm cm cm 
a asi G..— 10 2087.71 0 6S 
5 l 2080.49 } 9 10 2687.71 2 6o 
2 MSS, ST 70 » $ 2733. 86 18 
{ Lt Qu 208U, 20 3.9 71 
‘ 3 28% 14 4 89 72 l—1 2734.8 1.9 
} 2OS4 x4 73 
i), 32 74 15--4 1.6 
6 200 54 = (3: —3 1.7 
. #91. 40 , \2 2 2.1 
S OL. 74 76 ; 3 19 
q 2692. 50 77 2 $ 2737. 33 1.3 
10 2593, 13 73 
3 ast K. 9 2603, 3S 79 20-2 738, 48 8.7 
7 il LOIS. J | s rT) 24593. 38 27 x) 
- 12 2604. 12 SI 2 2 739. 72 1.5 
- } 7 it 2094, 83 2 0 SS 
° i4 miO5. 22 Ss $ 3 1.1 
7) 15 2606, 10 SS l 2 3.1 
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* Because so many absorption peaks are composites of several transitions 


numbers correspond to numbering of peaks in figure 1. 


The band center of the (0,1,1) band is some 40 
em™' lower than the Bailey, Thompson, and Hale 
value whereas it is some 70 em™ higher than the 
center of the absorption observed by Nielsen and 
Nielsen. 

The A-type band at 1906 em~? is in the region 
where one would normally expect the stretching 
fundamentals. The band type demands that this 
absorption be assigned to the (0,0,1) band. Such 
an assignment places the band center of v, nearly 
100 em™ lower than that reported by Nielsen and 
Nielsen. 


5. Discussion 


With the help of these new band centers it is 
possible to make some estimates of the (0,1,0) band 
center. The difference between the (0,1,1) and 
(0,0,1) frequencies is ~848 em™'. This difference is 
the band center of the (0,1,0) band minus X2,, the 
interaction between v. and »,. This interaction can 
be estimated from the corresponding interaction in 
HS as about 10 cm™ leading to a value of 858 
em™' as the frequency of (0,1,0). Using the (1,1,0) 
frequency, and as », the Raman [7] line at 1892 
em one finds 851 em~! which includes Yj). 

This frequency can also be estimated from the 
spectrum of crystalline [8] H.S and D,S which leads 
to a value of 855 em“. 

The existing data in this region by Nielsen and 
Nielsen runs out at 853 em™', so this absorption 
seems to represent only the R-branch of the bending 
fundamental. Attempts to analyze their data as 
an R-branch were unsuccessful. Their data show 
no regular series that can be identified with the main 
strong R-branch series. However, the three pieces 


, only those for which definite assignments could be made 


ire included in the table. Ling 


of evidence offered above strongly argue for a » 
band center close to 855 em! fully 80 em~'! lower 
than the value previously suggested [2]. 

These findings argue strongly for more work on 
the D.S spectrum, which will be undertaken shortly 
at the Bureau. 


In conclusion the authors thank R. E. Florin for 
preparing the D.S and Carroll Dannemiller for 
doing the least squares calculations on SEAC. 
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a Some Effects of Low Temperatures and Notch Depth on 
the Mechanical Behavior of an Annealed 
Commerciall re Titanium’ 

y Pu t 
Glenn W. Geil and Nesbit L. Carwile 
Unnotched and notched specimens (60° notch angle, 0.05-inch root radius and various 
notch depths) of an annealed commercially pure titanium were slowly strained to fracture in 
tension at +-150° or + 100° to —196° C, to reveal the combined effects of temperature and 
notch geometry on the tensile behavior of the metal. Impact tests were made on Charpy 
V-notch specimens at + 300° to — 196° C for a determination of the impact notch-toughness 
of the titanium. True stress-true strain relations were determined for the titanium in 
tension and a study was made of the effects of test temperature, stress system, and inter- 
stitial content on the mechanism of deformation and work-hardening characteristics of the 
metal. 
l. Introduction | comprehensive program at the Bureau to evaluate 
| the fundamental factors affecting the deformation of 
Several investigators [1 to 5]? have published | metals. The present paper summarizes the data 
- results of their studies to determine the mechanisms | obtained at low temperatures in Charpy V-notched 
involved in the tensile deformation of unnotched | impact and in tension on both unnotched and notched 
lV specimens of titanium and some of its alloys. The | specimens of selected notch depths of this heat of 
wer data obtained on single crystals [4,5] and on coarse- | titanium. 
grained polycrystalline specimens [1,2,3] indicate 
on that the tensile properties and the mechanisms in- 2. Material and Test Procedures 
rtly volved in the tensile deformation of alpha titanium 
are quite temperature sensitive and complex. <A 2.1. Material 
large proportion of the , deformation mn vests at All test specimens were prepared from bars proc- 
oa —196° C was by twinning [1,2]; the amount of | gcced from a single heat of commercially pure 
ra twinning, in general, decreased with increase iN | titanium (ASTM Designation B265-52T, grade 2). 
temperature. Single and duplex slip on the pris- | ‘The bars were supplied in the form of 1-in. rounds in 
matic planes also occurred in the deformation at | the hot-rolled and annealed condition; the bars were 
sae a ah mrense in temperature | fine ‘grained (ig. ‘The impurities aumeria 
eae Ti Me x Be senegal Cwty weed * “9) | Values are in percent by weight) determined * in the 
” also on pyramidal planes [1,4] and the basal plane [2]. titanium were as follows: Oxygen, 0.21; nitrogen, 
The amount of slip on pyramidal planes increased 0.04; hydrogen, 0.012; carbon, 0.04; silicon, 0.02; 
77 considerably as the temperature was raised to +500° iron, W; chromium, W; tin, W; aluminum, VW; 
J = : vr tal rd Qos es I! ep ee magnesium, VW; manganese, VW; nickel, VW; 
v7 ery Tew data are avaliable on the combinec effect copper, T: tungsten, not determined. ; 
ch of notches and low temperatures on the tensile 
behavior of titanium and its alloys; most of the 2.2. Specimens 
m. published data were obtained on titanium alloys a ; 
either in tests at different temperatures on specimens The forms and dimensions of the tensile specimens 
™ with a constant notch geometry [6,7,8], or in tests at | &re shown in figure 2. lhe reduced section (D=0.40 
24) room temperature on specimens of selected notch | 1.) of the unnotched specimens was slightly tapered. 
depths [9]. Results obtained at the Bureau in The notched specimens were machined with a con- 
n tension tests at different temperatures on a heat of | Stant root radius of 0.05 in. and a constant minimum 
annealed commercially pure titanium [10] contained diameter, d, of 0.35 in., thereby minimizing the size 
some data on notched specimens of selected notch | effect. The notch depth * was varied by changing 
depths. However, the quantity of titanium avail- | the diameter, D, of the cylindrical portion of the 
able for these tests was relatively small, and only an | Specimen adjacent to the notch. P 
exploratory investigation of the effect of low temper- ASTM standard (E23-47T) Charpy V-notch 
atures on the properties of notched specimens was | Specimens were prepared from one bar of the tita- 
possible. nium, 

An additional heat of commercially pure titanium 3 Chemical analysis determinations for nitrogen, carbon, and silicon were made 
was procured for use in the present study of the | py {he Anatytiel Chemistry, Secu of NBS: Determinations for een and 
combined effects of low temperatures and notch modified vacuum fusion method [11]. The spectrographic-analysis determina- 

. = tions designated by letters were made by the Spectrochemistry Section of NBS, 
the tensile behavior, as a part of a In general, W denotes weak (0.01 to 0.1%), VW denotes very weak (0.001 to 


geometry on 


This investigation was partially supported by the Materials Laboratory, 
Directorate of Research, Wright Air Development Center, USAF 
? Figures in brackets indicate the literature references at the end of this paper 


} 


0.01%), and T' denotes trace (0.0001 to 0.001%). The spectrographic method 
used was not sensitive for determination of tungsten, 

‘ Notch depth is expressed as the percent of cross-sectional area removed in 
machining the notch in the specimen and it is equal to 100 (.D?—d?)/ D?, 
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FIGuRE 2. 


Dimensions and forms of notched and 


Keo 


lensile specimens. 


innotched 


Microstructure of the annealed commercially pure titanium. 


1 part HF, 10 parts HCl, 10 parts HNO 


ind SO parts water 


2.3. Test Procedures 


The Charpy V-notch tests were made at tempera- 
tures ranging from +300° to —196° C in a machine 
of 224.1 ft-lb capacity and a striking velocity of 
16.85 ft/sec. The specimens, except those tested at 
room temperature (about +25° C), were immersed 
in a bath at the selected temperature for a minimum 
period of 30 min, and then quickly transferred to the 
impact machine and broken. The total time elapsing 
between the removal from the bath and the breaking 
of the specimens ranged from 3 to 4 sec and there was 
no significant change in the temperature of the 
specimens during this period. 

A detailed description of the tensile test equipment 
and the method of maintaining the selected tempera- 
ture is given in previous publications [12, 13, 14]. 
The tensile tests were made at temperatures ranging 
from + 150° or + 100° C to —196° C with notched 
and unnotched specimens in a pendulum hydraulic 
testing machine. The specimens were strained 
slowly with the deformation rate maintained at 
about 0.5 to 1.0 percent reduction in area per min- 
ute during the deformation beyond initial yielding 
The specimens, except those tested at room tempera- 
ture, were strained to fracture while completely im- 
mersed in an appropriate. bath maintained at the 
selected temperature. The specimens were im- 
mersed in the bath for a period of 30 min prior to the 
application of load. Simultaneous load and diam- 
eter (minimum diameter of unnotched specimens, 
diameter at base of the notch of the notched speci- 
mens) measurements were made throughout each 
tension test. 
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3. Results and Discussion 


3.1. Impact Tests 


The results obtained in the impact tests are sum- 
marized in figure 3. A transition temperature of 
about 110° © was indicated by the criteria of 
either (1) the steepest slope of the energy-tempera- 
ture curve or (2) the mean value of the energy 
absorbed at the lowest temperature, and at the high- 
est temperature at which complete fracture occurred. 
These impact data differed appreciably from those 
obtained previously [10] on another heat of com- 
mercially pure titanium with relatively small 
amounts of interstitials, which had a_ transition 
temperature about 40° C higher, and a considerable 
degree of notch toughness (10 ft-lb, energy absorbed) 
at —196° CC. This difference in the notched-bar 
impact behavior of these two heats of annealed 
commercially pure titanium may be attributed to 
the effects of the interstitial elements (carbon, nitro- 
gen, oxygen, and hydrogen) as the other impurities 
did not vary significantly. The observed lower 
notch toughness (energy absorbed) at temperatures 
below the transition temperature range of the present 
heat of titanium with the higher interstitials, and 
thus with the higher hardness and strengths, con- 
forms to the normal behavior for metals. However, 
the observed lower transition temperature for this 
heat of titanium is surprising, and no satisfactory 
explanation for this behavior can be given at the 
present time. As the grain size of the previous heat 
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fracturing the Charpy V-notch specimens 


of titanium with low interstitials was slightly finer 
than that of the present heat, the difference in the 
transition temperatures cannot be readily attributed 
to a grain size effect. 


3.2. Tensile Tests 
a. Flow Curves 


The true stress-true strain relations obtained in 
tension tests on unnotched specimens at different 
temperatures are presented in figure 4. The true 
stress °-true strain ® values that were obtained in 
duplicate tests at + 100°, +25°, and —196° C were 
fairly consistent and showed only slight deviations 
from the curves representing the average values at 
each temperature. The marked effect of tempera- 
ture on the strength of the titanium is shown by the 
increase in height of the true stress-true strain curves 
as the temperature is lowered. The initial resistance 
of the metal to plastic flow at — 196° C is more than 
three times that at +150° C. 

Logarithmic graphs of the true stress-true strain 
data (not shown) were sigmoidal. These data, like 
those obtained previously on another heat of titanium 
|10], and on other metals [13, 14], do not conform to 
the often postulated relationship, c=ké", for tensile 
deformation, where o and 6 are the true stress and 
true strain, respectively, and & and m are constants. 

The true stress-true strain relations obtained in 
tension tests at + 100°, +25°, —78°, and —196° C 
on circumferentially notched specimens are presented 


’ True stress, o, was determined by dividing the current load by the current 
minimum cross-sectional area of the specimen. 

® True strain, 6, was determined as the natural logarithm of the ratio of the 
original minimum cross-sectional area of the specimen Ao, to the current mini- 
mum cross-sectional area, A, 
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Figure 4. True stress-true strain relations obtained in the 
tension tests at different temperatures on unnotched specimens. 
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depth was increased above 70 percent. This is to | * ¢-smitian FRACTURE 











be expected, because the increase in triaxiality was 4 (1 
also relatively small as the depth was increased P tests 
from 70 to 87 percent. ‘ mari 
Both the stress concentration and the triaxiality ; mini 
induced in a notched specimen in a tensile test vary offse 
with the geometry of the notch. The variation of 1 was 
these factors with the notch depth for the titanium (a . “* How 
specimens of circular cross section containing a cir- ee ee ee eer a Se) . ined + obtained in th appl 
cumferential notch of 60° angle and 0.05-in. root ” inlininea~Conuatioaaniae:. °° (aa 
radius, is shown in figure 9. The values for the plas’ 
stress concentration factor A,, were derived from be h 
published stress concentration design charts [15]. valu 
The triaxiality, determined as the ratio S,/S;, anv a a ee high 
where S, and S; are the transverse and longitudinal stra 
stresses, respectively, is based on the method pro- + lowe 
posed by Sachs and Lubahn [16]; that is, S,/S ---*, dept 
1-(1/f), where R is the ratio of notch strength ’ to the F \ T 
tensile strength. The triaxialitv increases at a PE om ndigr ook Saeaet ium 
nearly constant rate with increase in notch depth to | +B yore sent 
about 20 percent, and then increases at a continu- | 5 [ ¢ (0.02 
ously decreasing rate with further increase in notch | & 9 | sligl 
depth. The stress concentration factor increases | @ LY? M-MAXIMUM LOAC , the 
very rapidly with increase in notch depth to about FE F- INITIAL FRACTURE | (70 
10 percent, less rapidly with increase in depth from ie and 
10 to 30 percent and then increases only slightly iy corr 
with further increase in the depth. s crea 
en Se. ae oer Pere See note 

’ Notch strength was determined by dividing the highest load attained on the : ae sia taden ' — , tem 
notched tensile specimen by the original minimum cross-sectional area at the NUE siete Je “o/* of 1 
notch In determining the values of #, the notch strength values are restricted FIGURE 8 ine sirens tree thenin wthitene aitained te the aes 


to those in which the specimen had sufficient ductility to attain a normal maxi . 
mum load condition, thet is, the slope of the load-extension curves becomes zero fension tests alt 196° C on notched specimens fact 
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Figure 9. Effect of notch depth on the stress concentration 
factor and the triaxiality of the tensile specimens use d. 


60° notch angle, 0.05-in. root radius, 


notch 


b. Effects of Temperature and Notch Geometry 


(1) Strength Indices. The results of the tensile 
tests on unnotched and notched specimens are sum- 
marized in tables 1 and 2, respectively. A deter- 
mination of vield strength, as normally based on 
offset values from autographic load-extension curves, 
was not feasible with the equipment emploved. 
However, an indication of the at which 
appreciable vielding has taken place is given by the 
values of true stress at a true strain (elastic plus 
plastic) of 0.02. These ‘0.02 stress’? values would 
be higher than the 0.2 percent offset vield strength 
values for unnotched specimens, particularly at the 
higher temperatures because the amount of elastic 
strain prior to yielding, in general, increased with 
lowering of temperature and with increase in notch 
depth. 

The relations between the 0.02 stress of the titan- 
ium and the notch depth and temperature are repre- 
sented by the surface delineated in figure 10. The 
0.02 stress at each test temperature increased, at a 
slightly decreasing rate, with increase in notch depth; 
the values for the specimens with the deep notches 
(70 and 87% notch depths) were approximately one 
and one-half times those for unnotched specimens at 
corresponding temperatures. A relatively large in- 
crease in the 0.02 stress was observed for both the 
notched and unnotched specimens with decrease in 
temperature from + 100° to —196° C; the strengths 
of the specimens at 196° C were greater than 
those of corresponding specimens at +100° C by 
factors of 2.5 for unnotched specimens and 2.7 to 


stresses 


( 
2 


$VV4S— HT 5 


ind 0.350-in. diameter at the base of the 


TABLE |. Tensile properties at + 150°, 
78°, and 196° © of unnotched 
annealed commercially pure titanium 


+ 100°, 4 
specimens 


fo 
25°, 


30°, 
of initially 


True stress True strain > 


Speci- __| Re- | Elon- 
men * Tem- Tensile duec- | ga- 
desiv- wra- strength At At At in- At Atin-| tion | tion « 
nation ture true | maxi-| itial | maxi-| itial of in 
strain | mum | frac- | mum | frac- | area | 2 in 
of 0.02) load ture load ture 
Cc iin?  tblin2 | thlin2 | Dline q% Qq 
B10 +150 60,000 50,000 69,009) 95,000 0.121 | 0.741 &2 44 
Bx +100 75,000 64,000) 85, 000/115, 000 123 774 57 40 
B6 +100 75,000 63,000) 86,000,117,000 137 770 87 42 
Bl +25 04,000 84,000) 106, 000'144,000, . 123 708 53 30 
Bg +25 94,000 84,000) 108, 000 145, 000 136 703 51 30 
B5 30 (111,000 100, 000 128, 000) 145, 000 140 602 5 26 
B2 78 124,000 113, 000) 146, 000) 179, 000 160 522 12 rai) 
B7 —~196 (169,009 164, 000/202, 000' 227, 009 176 330 29 20 
BA -196 169,000 157, 000) 198, 000/229, 000 159 . 358 31 22 


* All specimens were prepared from the same bar, B. 

> Includes clastic and plastic deformation. 

¢ Measurements based on punch marks in fillets 
section of 2.0 in. length 


Specimens had a reduced 


2.2 for the notched specimens (table 3). No simple 
relation between 0.02 stress and temperature is indi- 
cated as the curves (fig. 10) representing the varia- 
tion of the stress with temperature for specimens of 
the same notch depth have an inflection point near 
0° C. The inflection point may be an indication of 
a significant increase in the number of active slip 
plane systems and accompanying decrease in the 
stress for activation of slip with increase in temper- 
ature above 0° C, 

The observed relations between the notch strength 
of the titanium and the notch depth and tempera- 
ture are represented by the surface delineated in 
figure 11. The shape of this surface is nearly similar 
to that shown for the 0.02 stress (fig. 10). As shown 
in figure 12, the percentage increase in notch strength 
is approximately equal to the notch depth for 
depths up to 50 pereent. However, nearly all of 
the values for notch depths of 70 and 87 percent lie 
considerably below the dotted line representing the 
above relationship. This deviation may be attrib- 
uted partially to the relatively low rate of increase 
in triaxiality with increase in notch depths for 
depths greater than 50 percent (fig. 9). A consid- 
erable portion of the deviation of the values for the 
deep-notched specimens at —78° and —196° C 
may be attributed to the very small ductility of 
these specimens (table 2). 

The combined effects of temperature and notch 
depth on the true stresses at maximum load and at 
initial fracture are illustrated by the two surfaces 
delineated in figure 13. These two surfaces coincide 
at temperatures of —78° C and below, for notch 
depths of 10 percent or greater, and for all notch 
depths at —196° C. This coincidence is associated 
with the fracture of these specimens while the loads 
were still increasing. 

The true stresses at maximum load and at initial 
fracture increased greatly with decrease in tempera- 
ture (fig. 13). However, the relative increase in these 
values with decrease in temperature from +25° to 

78° or to — 196° C (table 3) is significantly smaller 
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Fiaure 10 Variation of the true stress at a true strain of O02 Figure 11. Variation of the notch strength with the notch 
with notch de pth and te m perature of the tensile Specimens de pth and te mperature of the tensile specimens, 
Tasie 2. Tensile properties at 100 25°, 78°, and 196° © of notched specimens (60° notch angle, 0.05-in. root radius 
of initially annealed, commercially pure titanium 
Ratio of true stress of Ratio of true 
rue stress at notched to unnotched rue strain ! strain of notched 
Speci Ratio of specimens at at to unnotched Re 
men * rem- Notch Notch notel specimens at tuc 
desig pera depth D d rengtl strengtl tio 
nation ture to tensile if 
strenetl Prue Max Initial rrue Max Initial Maxi Initia! Max Initial rea 
tra mum fracture strain mum iri mum ire num irae 
of 0.02 load of 0.02 load ture load ture load ture 
( in nN in f Ly lhiin 
(i4 +100 0.360 0.35) 77, 000 lf 65, 000 Sv. 000 119. 000 1.02 1o2 1.03 124 0. 631 OF 82 4 
gil +100 10 $70 3m) 81.000 1.08 64, 000 43, 000 121, 000 Lol Lo 1.04 136 4 1 O48 wy 42 
qi th i th 70 +) 8? O00) Lo TO. Oye 44 OOD 123, 000 Lo 111 om 132 2? 1.02 is 4) 
ais +1%) wv) 418 $f) a O00 1.30 a O00 116, 000 130, 000 130 1. 3 1.12 173 $i 5] 17 4 
Gl +100 ” $05 +) TOs, O00 1.44 43, 000 126, 000 142. 000 1. 46 1.47 1. 22 1 $s lv 17 ) 
in +100 ") $05 +) 109, 000 1.45 41, 000 124. 000 142. 000 1.43 1. Al 1.22 ied 332 1.2 4 2s 
Gis +100 70 ny +0) 114. 000 1 TOO. O00 136, 000 14. 000 | 159 133 135 4 rot ‘1 $2 
17 +1) 70 639 3%) 12. 000 1“ 11, OOO 138. 000 154. 000 1 oS Ll 133 130 359 ] ; 17 vw) 
Ell +100 x7 (Hy) tu) 125, 000 147 107. 000 143, 000 157, 000 1. 67 1. 47 l . 35 s1s 14 41 n» 
D4 +2 it) su) 89, 000 105 42, 000 112. 000 140, 000 1 10 0). 97 11S 1sS Tu “iu $! 
Dis +2 0 370 +) 105, 000 1.12 4), 000 1 000) «141.000 1.14 11 is 1} 446 117 i) $2 
(i +25 10 $70 $3) 103 000 1.10 41. O00 120. 000 139. 000 1. ON 1.1 ) 155 $41 i is 30 
FS 2 0 $70 $a) 1034, 000 LW SA OD 123, 000 144. 000 105 1. te Loo 172 4] ; 5 5) 
D7 +2 uM) 11S 30) 125, 000 1. 33 114, 000 146, OO) 156, O00 1. 36 1. 37 1. US 151 Zt 1.7 2 ee) 
qi 1) 11s 3M) 1), O00 1. Js 14, O00 145. 000 154, 000 1. 24 1 36 1.07 142 Di) 1. 4s $ a4 
Fi + wu) 118 $M) 123, OOM +] 106,000 «148,000 156,000 1. 2h 1.39 1. OS Is? 2%) 1.41 4 
14 + 1) $05 $0) 141, 000 1M 130. 000 159, 000 177. 000 1. 55 1.44 1. 23 122 2m 4 ¢ 3 
15 1) 105 350 140, 000 1.49 120,000 164, O0F 172, 000 1.43 1. M4 1.19 15S 218 22 $] a 
Fa +2 ” vv 3m) 142, 000 1 5} TTS, O00 166, 000 174, 000) 1.40 1 120 15s »~? 1. 22 l 23 
D6 +2 70 Hay 3m) 182, 000 1.72 138. 000 183. 000 195, 000 1.4 1.72 1. 35 122 27 O.u4 hal 21 
F21 +2 70 139 wh) 1 000 1. 126, 000 178. OOO 192, 000 1”) 16 1.33 130 2 1 oOo tf 44 
Div 2 7 u4O 3M) 167, 000 1. 78 135, 000 TSA. OM) LL 141 1.77 1. 37 11s ISS 0.41 2 1S 
As 2 x7 Gn) uM) 155. 000 1. 65 134. 000 ISO, O00 PSS, 000 Low 169 120 146 Jt 1 is $1 X 
N6 +25 x7 uu) 3) 160, 000 2 132, 000 IST. O00 145, 000 1.457 1.70 l 12 222 yy l Zl 
IS 7s 5 win) ") 128, 000 1.08 117, 000 154, 000 162. 000 1.04 106 0.91 17s 252 lL. 11 is 24 
Is 7s 10 {70 3%) 134, 000 1. OS 128), 000 14. 000 164, 000 106 1.12 v2 203 D3 1. 27 4 1Y 
12 7s 0) 418 i) 14, 000 1. 26 142. 000 177. 000 1. 26 ay 132 2 12 
li 7s ») $05 i) 177, 000 1.43 146, OOO 200, 000 1.28 1.12 119 23 ia 
4 7s 70 aay +0) 199. 000 1.41 143. Oo 773 OO) 1. 35 1. 25 115 22 12 
Al2 7s a7 won 350 200, 000 1. 61 173, 000 217, 000 1.53 1. 21 ON6 17 It 
N7 7s 87 wun) 3m) 200, 000 1. #1 150, 000 72? 000 133 1. 24 102 AY 10 
Fl 16 5 3) 300 177. 000 1.0 164. 000 199. 000 1.05 O87 x 11s 4 ll 
F2 1 10 $71 350 183, 000 1. ON 182, 000 199, 000 1.13 x7 Os! 24 s 
F3 106 30 418 4) 213, 000 1. 26 210, 000 221. 000 1.31 “7 “ O39 a 5 
Fi7 106 ”) 195 tw) "249, 000 1.47 215. 000 263. 000 1.34 1.15 « On It t 
All 196) ”) 195 3M) 243, 000 1.44 225, 000 255, 000 1. 40 1.12 « O52 15 6 
F4 106) 70 Hae 0 258, 000 1. 53 240, 000 270, 000 1. 43 1. 18 ‘ O46 13 5 
49 196 a7 4) tO 257, 000 1. 52 235, 000 266, 000 1. 46 1.17 « 037 11 t 
N& 106 7 a $50) 256, 000 1. 52 239, 000 24, 000 1. 48 1.16 “ O31 ov i 





* All specimens prepared from the same bar are designated by the same letter 
b+ Diameter of cylindrical portion of the specimen adjacent to the notch 
¢ Diameter of specimen at the base of the notch 
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e Determined by dividing the highest load attained by the original minimum 
ross-sectional area 

f Includes elastic and plastic deformation 

« Fractured before the slope of the load-extension curves became zero 
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Ctracture at T2 


Effect of test temperature and notch depth on relative values of some strength indices for initially annealed commercially 


7 » 2 
TABLE ». 
oi~o.02 at T 7Maximum load at T2 
a pat ‘I Maximum load at T; 
Notch 
depth 
Ti. ( +25 +25 +25 +100 +25 +25 +25 
H °< +100 78 196, 196, +100 78 — 196 
0 0. 76 1.35 1.91 2. 53 0. 80 1. 36 1. 87 
5 71 1.27 1. 84 2. 60 78 1. 38 1.78 
10 73 1.31 1.99 2.72 77 1.35 1. 64 
20) 77 1.31 1.94 2. 53 79 (>) (>) 
50) 75 1.19 1.79 2.39 78 (t (>) 
70 76 1. 16 1.74 2. 29 76 (>) (>) 
7 80 1.21 1.77 2. 21 78 (t (>) 


» Ratio of tensile strengths. 








Notch strength at T, 


Tiracture at T; Notch strength at T; 


+100 +25 +25 +25 +100 +25 +25 +25 +100 
— 196 +100 —78 196 — 196 + 100 —78 —196 —196 
2 0. 80 1.24 1. 5S 1.97 * (0. 80 * 1.32 * 1.80 ® 2.25 
2. 29 85 1. 16 1.42 1. 67 78 1.29 1.79 2.30 
2.13 6 1.16 1. 41 1.63 78 1.29 1.76 2. 24 
(>) M4 1.14 1. 42 1. 70 80 1.27 1.73 2.17 
(t 81 1.15 1.49 1. 82 77 1. 26 1.7 2. 27 
(t 80) 1.15 1. 40 1.75 75 1.25 1. 62 2.16 
(t SI 1.13 1.37 1.69 78 1.24 1.59 2. 08 

















» Fractured at —78° and —196 C before the slope of the load-extension curves became zero 
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than the corresponding relative increase in the 0.02 
values; this feature can also be attributed 
mainly to the decrease in ductility of the specimens 
with decrease in temperature. 

The true stress at initial fracture at temperatures 
of +25° C and below, decreased with increase in 
notch depth up to 10 percent, then increased with 
increase in depth through the range of 10 to 70 per- 
cent, and thereafter changed only slightly with 
further increase in depth to 87 percent. The de- 
crease for 5 and 10 percent depths can be attributed 
directly to the decrease in the true strain at initial 
fracture associated with the stress concentration at 
the root of the notch. The decrease in the true 
stress at fracture due to this factor is greater than 
the small increase in the true stress associated with 
the induced triaxiality for these shallow notches. 
For notch depths of 30 percent and greater, the 


stress 


relative effect of the above two factors is reversed, 
and the triaxiality is the predominant one, as the 
stress concentration factor and the ductility of the 
specimens changed only slightly, whereas the tri- 
axiality and flow stress increased considerably, with 
the increase in notch depth. 

(2) Duetility. The effects of temperature and 
notch depth on the true strains at maximum load 
and at initial fracture are shown in figure 14. The 
true strains at maximum load for the specimens that 
attained a normal maximum load condition did not 
vary significantly either with temperature or notch 
depth. A slight, although probably insignificant, 
maximum is exhibited in the curves showing the 
effect of notch depth on the true strain at maximum 
load at temperatures of +25° and +100° C. The 
true strain at initial fracture decreased continuously 
with decrease in temperature below + 100° C. 
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Figure 14. Variation of the true strains at maximum load 


and at initial fracture with the notch depth and temperature of 


the tensile specimens. 


Moreover, the true strain at initial fracture decreased 
very rapidly with increase in notch depth to about 
10 percent, less rapidly with increase to 30 percent, 
changed only slightly with further increase in depths 
to 50 and 70 percent, and then decreased slightly 
with increase in depth to 87 percent. The pro- 
nounced dependence of the ductility of this metal on 
the notch depth can be attributed to variations in 
the stress concentration factor and triaxiality (fig. 9). 
The stress concentration factor increases very rapidly 
with increase in notch depth to about 10 percent, 
and less rapidly with increase in depth from 10 to 30 
percent, whereas the triexiality for these shallow- 
notched specimens is relatively smell. Therefore, it 
can be postulated that the large decrease in ductility 
of these titenium specimens (fig. 14) with increese in 
notch depth to 5, 10, and 30 percent is due meinly 
to the detrimental effect of the stress concentration 
and gccompenying stress and strain gradients in the 
region near the base of the notch. The triaxiality 
of the stress system also reduces the ductility of 
specimens but its effect for these shallow-notched 
specimens is relatively small. However, for the 
deep-notched specimens, the effect of the relatively 
large triaxiality on the ductility is significant as is 
shown in figure 14 by an appreciable decrease in the 
true strain at initial frecture et temperatures of 
+ 100°, +25°, and —78° C as the notch depth is 
increased from 70 to 87 percent; the change in the 
stress concentration factor with notch depth in this 
range is very small, and the triaxiality apparently is 
the predommant factor involved in the initiation of 
* the fracture at these temperatures of the specimens 
with notch depths of 70 and 87 percent. 

Some corroboratory evidence in support of the 
above views is provided in the photographs of the 


fracture surfaces of the notched specimens pre- 
sented in figure 15. The fractures at + 100°, +959 
and —78° C of specimens with notch depths of 5 
10, 30, and 50 percent, and the fractures of all the 
notched specimens at —196° ( were initiated at 
points of high-stress concentration at, or near, the 
root of the notch. Conversely, the fractures aj 

100°, +25°, and —78° C of specimens with 
notch depths of 70 and 87 percent were initiated 
in the region of high triaxial stresses near the axis 
The total deformation of the deep-notched speci- 
mens at these temperatures was considerable: jt 
apparently was great enough to reduce the stress 
concentration near the root of the notch so that 
the effect of the high triaxial stresses was the pre- 
dominant factor in the initiation of these fractures. 
At —196° C, the total deformation of the deep- 
notched specimens was very small. It was not 
great enough to reduce, to any significant extent, 
the stress concentration near the root of the notch, 
and the stress concentration remained as the pre- 
dominant factor in the initiation of the fractures. 


c. Work Hardening Characteristics 


The rate of work hardening,* do/dé, of the tita- 
nium was affected by the temperature and noteh 
depth. This is indicated qualitatively in figures 4 
to 8 by an increase in the slopes of the curves at a 
given true strain with either a decrease in tempera- 
ture or an increase in notch depth. 

The relations observed between the rates of work 
hardening and true strain of unnotched specimens 
of the titanium at different temperatures are sum- 
marized in figure 16. The rates of work hardening 
at constant true strains, except at the small true 
strains in the region of initial yielding, increased 
continuously with decrease in temperature over the 
range of +150° to —196° C. The rates of work 
hardening of the specimens decreased greatly with 
increase in true strains for strains less than those at 
the maximum loads. The decrease in the rate of 
work hardening with increase in strain, above that 
at maximum load was still considerable at the lower 
temperatures, but it was very slight at + 100° and 
+ 150° C, 

The effect of notch depth and the accompanying 
induced triaxiality on the rate of work hardening of 
the titanium is illustrated quantitatively in figure 
17 for the tension tests at +25° C. The rate of 
work hardening of the metal increased with increase 
in the notch depth and decreased with increase in 
true strain; the rate of decrease was relatively small 
with strains in excess of about 0.3. Similar results 
(not shown) were observed for the tests at 100°, 

78°, and —196° C; however, at —78° and —196° 
(* the true strain range was much smaller due to 
the limited ductility of the notched specimens at 
these temperatures. In the tests at 100° ©, 
considerable scatter was observed in the rates of 
work hardening of the notched specimens at strains 

* The rate of work hardening, determined as do/dé, at any true strain is the slope 
of the true stress-true strain curve at that strain. For strains of an unnotched 
specimen beyond that at maximum load, the increase in the flow stress resulting 


from the triaxiality induced by the necking of the specimen is included in the 
value of do/dé 
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rk Ficure 15. Effect of test temperature and notch depth on the appearance of the fracture surfaces of the tensile specimens. 
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above that at maximum load, and the effect of 
notch depth on the rates of work hardening was 
considerably less than that observed in the tests 
at +25° C. 

The dependence of the rates of work hardening of 
the titanium on the stress system may be explained 
by the hypothesis, suggested in a previous paper 
{10}, as follows: Deformation in tension under a 
multiaxial stress system, such as the triaxial stresses 
induced in these notched specimens, occurs by the 
movement of dislocations on e# greater number of 
slip planes within the individual crystals of the poly- 
crystalline specimen than is the case for deformation 
under uniaxial tension, and more duplex or cross 
slip will also occur with deformation under the multi- 
axial stress system. The number of intersections of 
dislocations, with the accompanying formation of 
jogs in the dislocations and point defects in the metal, 
increases with increase in the triaxiality and accom- 
panying increase in cross slip. Energy is required 
to move one dislocation across another dislocation. 
Moreover, the energy required to move a dislocation 
with a jog is greater than that required to move the 
dislocation prior to the formation of the jog. The 
point defects formed may also act as barriers to the 
movement of dislocations. These factors tend to 
increase the rate of work hardening of the metal 
(17, 18, 19}. 


d. Some Effects of Interstitials 


Recent investigations [20, 21, 22] have shown that 
the mechanical properties of high-purity alpha tita- 
nium are affected greatly by the interstitial elements 
(carbon, nitrogen, oxygen, and hydrogen). Carbon, 
nitrogen, and oxygen strengthen the metal and elso 
tend to decrease its ductility. Hydrogen is generally 
believed to be one of the main causes of embrittlement 
of the metal. The main difference in the chemical 
composition of the present heat of titanium and the 
of the heat used in the previous study [10] at the 
Bureau is in the amounts of the interstitial elements: 
the interstitial contents (percent by weight) of thes» 
two heats of the titanium were as follows: 


Carhon Oxryoen Nitrogen Tlydroge 
Previous heat 0. 02 0. 07 0. 02 0. 006 
Present heat 04 21 . 04 012 


A comparison of some of the strength and duc- 
tility indices of unnotched specimens of the two 
heats of titanium studied at the Bureau is given in 
table 4. 

The 0.02 stress at 100°, +25° 78°, and 

196° © (table 4) of the present heat of titanium 
ranged from 12 to 30 percent greater than that of 
the previous heat of titanium with its lower inter- 
stitial contents. It may be postulated that this 
difference in the 0.02 stress values of the two heats 
is due mainly to the greater density of “Cottrell 
clouds” (clusters of interstitial atoms around the 
dislocations) [17] in the annealed metal with the 
higher interstitial contents; more energy is required 
for the initiation of the movement of a dislocation 
with a surrounding cluster of interstitial atoms than 
for that of a similar dislocation not partially anchored 


TABLE 4. Some comparative strength and ductility values of 
the unnotched tensile specimens of the present and previous 
heats [10] of initially annealed, commercially pure titanium 


rem perature +100° C +25 196° C 


Ratio of true stress at a true 
strain of 0.02 for the present 
heat to that for the previous 
heat 1.12 1. 22 1. 27 1.30 

Ratio of true 
stress ata true 
strain of 0.15 


+» {Present heat 1. 32 1.3 1.27 1. 2¢ 

to that at 0.02.) previous heat 1.31 1. 32 1.31 129 
Ratio of true strain et initial 
fracture for the present heat 

to that for the previous heat 1.25 1.03 0.70 0.44 


Ratio of true stress at a true 
strain of 0.02 at ~196° t0 that present heat... 1.90 
” _ | Previous heat 1.75 
. ne _ . . . 
by interstitial atoms. The binding effect of the 


surrounding interstitial atoms on a_ dislocation 
increases with decrease in the rates of diffusion of the 
interstitial atoms such as accompanies a decrease 
in temperature of the metal. This feature is indi- 
cated in table 4 by the significantly higher ratio of 
the 0.02 stresses at — 196° and + 25° C, respectively, 
for the present heat of titanium with its higher 
interstitial contents. 

The interstitial atoms also greatly affect the 
resistance of the metal to the continued movement 
of dislocations and thus affects its work hardening 
characteristics. Evidence of this effect is provided 
indirectly in table 4 by the ratios of the true stresses 
at true strains of 0.15 and 0.02, respectively, for 
these two heats of titanium.’ These ratios, which 
may be considered as measures of the relative work 
hardening during this range of deformation, are 
nearly the same for both heats of the metal. How- 
ever, as the 0.02 stress was greater for the heat with 
the higher interstitial contents, the increase in the 
resistance to flow during the deformation of these 
specimens was greater than that of specimens of the 
other heat of titanium. The slight increase ex- 
hibited in these ratios for both heats of titanium as 
the temperature was increased may also be inter- 
preted as indirect evidence of slip along more pris- 
matic planes, and including some pyramidal planes 
and possibly basal planes at the higher temperatures. 
This would be accompanied by an increase in the 
amount of cross slip and the number of intersections 
of dislocations, thereby increasing the work harden- 
ing of the metal. 

The effect of the interstitials on the embrittlement 
of the titanium at low temperatures is indicated in 
table 4 by the ratios of the true strain at initial 
fracture of the specimens. The ductility at room 
temperature of specimens of both heats was about 
the same. However, the low temperature embrittle- 
ment of the titanium by the interstitials was very 
pronounced as shown by the large decrease in these 
ratios as the temperature was lowered to —78° and 

196° © 


rhe upper range of the strain for these ratios was limited to 0.15 in order to 
eliminate any effects that might be attributed to triaxiality accompanying the 
necking of the specimens 
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4. Summary and Conclusions 


A study was made of the behavior at low temper- 
atures of unnotched and notched tensile specimens | 
0° noteh angle, 0.05-in. root radius, 0.35-in. 
diame ter at root of notch, and selected notch de »pths) 


| 
if 
| 


and Charpy V-notch impact specimens of a heat of | 


annealed commercially pure titanium. 

The transition from ductile to brittle behavior in 
the Charpy impact tests occurred within the range 
of +125° to +80° C; only a small amount of notch 
toughness (3 to 4 ft-lb energy absorbed) was retained 
by the metal at 196° C 

The resistance to deformation of unnotched and 
notched tensile specimens increased continuously as 
the temperature was lowered within the range of 

-150° to —196° C. The strength of the titanium 
was also affee ted gre atly by the stress system induced 
by a notch; stre ‘ngth indices, such as true stress at a 
true strain of 0.02 notch strength, and true stresses 
at maximum load and at initial fracture, increased 
greatly with increase in the induced triaxiality 
accompanying an increase in the notch depth. 

The ability of the titanium to deform under 
—— stress apparently was not affected greatly 
by lowering of temperature, as the true strain at 
maximum load in the tension tests on unnotched 


specimens increased only slightly as the temper- 
ature was lowered within the range of +150° to 


196° C 

The ductility of the titanium under multiaxial 
stresses, however, was very dependent on temper- 
ature. The true strain at initial fracture of notched 
specimens of constant notch depth decreased con- 
tinuously with lowering of temperature. Moreover, 
the ductility of notched specimens decreased very 
rapidly at all test temperatures with increase in 
notch depth to 10 percent, decreased less rapidly 
with increase in depth to 30 percent, changed very 
little with further inerease in depth to 70 percent 
and then decreased slightly with increase in depth 
87 percent. The concentration factor 
increases very rapidly with increase in notch depth 
to about 10 percent, less rapidly with increase in 
depth from 10 to 30 percent, and then rather slowly 
with further increase in notch depth, whereas, the 
induced triaxiality is relatively small for the shallow- 
notched specimens but increases to a fairly large 
value for the deep notch specimens. Thus, it may 
be concluded that the decrease in tensile ductility 
of the notched specimens apparently depends mainly 
on the and strain gradients accompanying 
the stress concentration near the root of the notch, 
and only to a minor extent, except perhaps for very 
deep notches, on the induced triaxiality. This 
conclusion, however, should be considered as ten- 
tative, as the change in notch geometry in this 
study was limited to a variation of the noteh depth 


to stress 


stress 


of notches with constant notch angle, root radius, 
and minimum diameter at the root of the notch. 


Additional experimental data involving other varia- 
tions in the geometry of the notch, such as notch 
angle and root radius, are needed to definitely estab- 


lish the relative effect of the stress concentration 
factor and triaxiality on the ductility of the titanium. 

The rates of work hardening of the titanium speci- 
mens under tension depended on the temperature 
and the stress system; the rate of work hardening 
increased with decrease in temperature and with 
increase in the triaxiality of the stress system. 

A comparison of the tensile data obtained pre- 
viously on another heat of annealed commercially 
pure titanium of low interstitial contents, with 
those obtained in the present investigation shows 
that the strength and work hardening increase and 
the ductility at temperatures below + 25° C de- 
creases with increase in the interstitial contents of 
titanium, 
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Determination and Use of the Sag Point as a Reference 
Point in the Heating of Glasses 
Sam Spinner, Given W. Cleek, and Edgar H. Hamilton 


when heated, undergo 
over a rather wide temperature range. 


Glasses, 


a gradual continuous change from solids to liquids 
As a consequence, it has been found necessary to 


define certain temperature reference points in this transition range so that different glasses 


may be compared with each other. 
stages in the annealing or forming process. 
point, is deseribed. 


These points are chosen to correspond to established 
A new reference point, designated as the sag 
rhe experimental determination of this point is rapid and simple. 


Also, it gives useful information concerning the annealing temperatures and ease of formation 


of glasses of widely varying compositions. 


The sag point is defined as the temperature at which a glass fiber 0.5 to 0.8 millimeter 
in diameter, horizontally supported at 's-inch intervals, will sag under its own weight in 


25 +5 minutes 


1. General 


One of the more widely accepted definitions of 
glass is that it is a material that has cooled from the 
liquid state without undergoing a change in phase 
(ij! Though glass may be regarded either as a 
supercooled liquid or an amorphous solid, the latter 
designation is more appropriate to stress the fact 
that ostensibly glass behaves like a solid, i. e., it 
retains its shape and reacts elastically to shear 
and tensile The process of transition 
from liquid to glass on cooling or from glass to liquid 
on heating is a continuous one. There is no sharply 
defined temperature, corresponding to the melting 
point of a crystal, at which it can be said that glass 
melts. Glass technologists have found it necessary, 
therefore, as a practical matter, to define certain 
reference pomts in the heating or cooling process 
so that reasonable comparisons can be made between 
glasses of different compositions, similar to the com- 
parison of crystals by their melting points. Because 
the reference points for glass are not sharply defined, 
their selection is, in a arbitrary. Their 
choice, however, is governed by established pro- 
cedures in the annealing, forming, or other stages 
in the production process. 


stresses. 


sense, 


Glasses of different compositions have different 
temperature-viscosity relationships. Hence, glasses 
whose reference points might correspond in one tem- 
perature region, such as the annealing range, will, 
in all probability, not have corresponding points in 
another region such as the forming range. Consider- 
ations of this nature show why more than one refer- 
ence point is necessary. 

Littleton and Roberts [2] and Littleton [3] made 
the original contributions to the establishment of 
reference points by defining successively the anneal- 
ing and softening points. Significant refinements 
and additions have been made by Lillie [4]. The 
American Society for Testing Materials [5] defines 
several such points now in use. These points are 
usually defined in terms of a specific phenomenon 
associated with a set of experimental operations, and 


Figures in brackets indicate the literature references at the end of this paper 


also in terms of the viscosity of the glass at the 
temperature of this phenomenon. For instance, the 
deformation point is defined as ‘“‘the temperature 
observed during the measurement of expansivity by 
the interferometer method at which viscous flow 
exactly counteracts thermal expansion. The defor- 
mation point generally corresponds to a viscosity in 
the range from 10' to 10" poises.”” The point 
defined by an operational procedure is preferable 
because the experimental procedure can usually be 
reproduced more accurately than the viscosity can 
either be calculated from the experimental condi- 
tions, or measured at the same temperature by some 
independent method. Furthermore, the viscosity 
frequently is not the only determining parameter. 
In the case of the softening point, for instance, 
the density of the glass must also be considered. 
The kinematic viscosity, then, would be the relevant 
factor. The situation is further complicated by the 
fact that at a lower temperature a considerable time 
may be required before the glass comes to the equili- 
brium viscosity for that temperature [6, 7]. This 
fact has led Lillie to regard the strain point, which 
occurs in such a low temperature region, to be of 
limited significance [4b]. It may also be mentioned, 
in this connection, that Adams, in his well-known 
paper on the annealing of glass [8], doubted whether 
the viscosity, which is designated as a megoscopic 
property, was really related to the release of stress, 
which involves only microscopic processes. 

In the Bureau’s Glass Section laboratory, it 
became necessary to devise a reference point that 
would serve as a reliable guide for annealing a large 
number of glasses of widely varying compositions. 
In order to be serviceable, it was also necessary for 
the determination of this point to be rapid and 
reproducible. 

It seems desirable to summarize some of the 
requirements, which a good reference point should 
satisfy, as follows: 

(1) The apparatus and procedure should be simple 
so that they may be easily duplicated by other 
laboratories. Methods that require elaborate speci- 
men preparation, for instance, would be undesirable. 
In general, the fewer qualifying specifications that 
need be made, the better. 
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The method should yield reproducible results 
In order to fulfill this most important requirement, 
the effect sought for should vary sharply with tem- 
perature and, conversely, should be relatively in- 
sensitive to extraneous factors. If, for instance, the 
specimen measured is in the fiber form, the thickness 
of the fiber should not be critical. 

(3) The determination should be rapid. This 
condition is especially important when the specimen 
is used as a control test for some portion of the 
manufacturing process. 


(4) The value of such a point is obviously en- 
hanced, if significant additional information is 
> 
obtained. 


A reference point satisfying these requirements was 
suggested by C. A. Faick, formerly of this Bureau, in 
1948, and has been utilized successfully ever since. 
This point is hereinafter designated as the sag point 
It is about 40° C higher than the deformation point, 
usually regarded as the upper limit of the annealing 
region. In the technique of determining the sag 
point, it was found that additional significant in- 
formation could be gained as to the feasibility of 
making sizable melts of hitherto untried com positions. 
This will be discussed later. 

It is the purpose of this paper to define the sag 
point, to describe the experimental equipment and 
procedure for its determination, and to show some 
typical results with glasses of different compositions. 
The sag points are compared with the deformation 
and softening points. 


2. Experimental Technique 


The technique about to be described for deter- 
mining the sag point bears some similarity to that 
originated by Silverman [9] and refined by Grauer 
and Hamilton [10] for measuring the liquidus tem- 
peratures of glasses. 


2.1. Equipment 


The basic equipment consists of a temperature 
gradient furnace, shown in figure 1, and a platinum 
specimen holder, figure 2. The platinum holder 
is 6 in. long with supports for the fibers spaced 
approximately }y in. apart. 

The furnace is essentially the same as that de- 
scribed by Grauer and Hamilton [10] except that the 
heating coils are of Nichrome rather than platinum 
wire. Figure 3 shows a typical curve of the tempera- 
ture distribution along the length of the furnace 
The temperature range, as well as the gradient, 
can be altered considerably by proper adjustment 
of the power controls. 


2.2. Procedure 


A specimen is drawn in the form of a fiber, 0.5 to 
0.8mm in diameter. Suitable fibers can be prepared 
.by drawing from the lip of the crucible in which r 
‘glass was melted. This is a convenient me ‘thod « 
obtaining fibers directly from the melt, even in ios 
cases where the devitrify readily. The 
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specimen is placed lengthwise on the holder (see 
fig. 2A) and the holder is inserted into the core of 
the furnace. After 25 min the holder is withdrawn. 
If the furnace has ioe properly adjusted, that see- 
tion of the fiber that was in the hotter region will 
have sagge “l between the supports, whereas the sec- 
tion of fiber that remained in the cooler region will 
not have sagged. The temperature at the first 
compartment, measured from the open end of the 
furnace, in which the glass sagged between the sup- 
ports (fig. 2B) is the temperature of the sag point. 
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The furnace is properly adjusted when the sag 
int occurs in the linear gradient region of the tem- 
perature distribution curve, shown between points 
4 and B in figure 3. A probe thermocouple is then 
inserted into the furnace up to this sag point distance 
and the temperature 1s measured. 


2.3. Definition 


Following the above procedure, the sag point is 
defined as the temperature at which a fiber, 0.5 to 
08 mm in diameter, horizontally supported at 
vin. intervals, will sag under its own weight in 
95 +5 min. 

2.4. Precision of Results 
a. General Considerations 


Under the conditions given in figure 3, repeated 
measurements with different specimens of the same 
glasses, but with the position of the holder varied 
slightly in the furnace, indicated a reproducibility 
of the order of +7° C. This is the limiting preci- 
sion for the determination of the sag point and fol- 
lows from the condition that 's in. of furnace length 
corresponds to a temperature change of about 15° C. 
Increased precision of results may be obtained by 
decreasing the over-all temperature gradient and 
hence the gradient per compartment. In this case, 
the following factors must be given greater consid- 
eration: Maintenance of a uniform distance be- 
tween supports of the holder, the fiber diameter, 
and measurement of the temperature. It may also 
be noted here that it was not found necessary to 
reinsert the platinum holder in the furnace when 
measuring the temperature of the sag point with 
the probe thermocouple. Within the experimental 
precision of the measurement, no difference was 
observed whether the holder was in or out of the 
furnace. The apparatus and method described here 
were found to give the speed, convenience, and pre- 
cision sufficient for the purpose of the determination. 


b. Effect of Fiber Diameter 


As might be expected, the temperature of the sag 
pomt was found to increase somewhat with increased 
fiber diameter. Within the precision used here, no 
significant variation was found within the given fiber 
diameter tolerance. The following table for one 
glass, BSC 517, shows typical data relating fiber 
diameter to sag point. 


Diameter Sag point 


mm ( 
0.50 to 0.84 650 
52 to .69 650 
50 to .60 657 
55 to .70 657 
62 to .76 650 


It may also be noted that the diameter for any single 
fiber may vary over the specified range without 
affecting the reproducibility of the result. 


c. Effect of Time 


The sag point was found to be relatively insensitive 
to time in the furnace within rather broad limits. 
Although longer times may be used safely, less than 
20 min may not give reproducible results. The 
time is not critical within the specified 25 +5 min. 
One of the conveniences of the method consists in 
the fact that it is not necessary for the operator to 
observe the specimen during the time of test, as is 
usually the case with other reference point determina- 
tions [11]. 


3. Use of Sag Point in Problems Pertaining 
to Glass Formation 


As was mentioned earlier, significant information 
could be obtained in the process of determining the 
sag point as to the feasibility of cooling sizable melts 
of new compositions to form glasses. In considering 
the problem of the glass-forming possibilities of a 
melt of a new composition, one is interested not only 
in whether the melt will cool as a glass, but whether 
glass will form in sufficiently large quantity so that 
sizable articles can be fabricated from it. Once the 
glass has formed, one is further concerned if such 
articles can be annealed without danger of devitri- 
fying. It is necessary to recall that the liquidus 
temperature is the temperature above which crystals 
cannot form. The temperature of maximum rate of 
crystal growth is usually 50 to 100° C below the 
liquidus temperature [10]. This region of maximum 
crystal growth is of the greatest concern in the 
practical problem of cooling a melt to glass. As 
the temperature is decreased below this region, the 
rate of crystal growth decreases rapidly because the 
increased viscosity associated with lower tempera- 
tures counteracts the tendency toward crystal 
growth. For practical purposes, this rate becomes 
zero at sufficiently low temperatures. Many compo- 
sitions that will form glasses in small quantities, or 
when the melt is poured on a cold plate, will devitrify 
to a considerable extent when an attempt is’ made 
to make these same compositions in large quantities 
for melts of about 500 g or more, because the center of 
the cooling piece must remain at temperatures in this 
region of maximum crystal growth for considerable 
periods of time. Thus, when specifying whether a 
particular composition will form a glass, the experi- 
mental conditions under which glass formation 
occurs must be specified [12}. 

In the technique of determining the sag point, a 
portion of the specimen is usually held in the tem- 
perature region where crystal growth may occur. 
The presence of erystals on withdrawal from the 
furnace, and their relation to the sag point, gives an 
indication of the tendency of the particular glass to 
devitrify in this critical region. As a consequence, 
the feasibility of cooling large melts of a given com- 
position as glass can be predicted. 

The one extreme would be a specimen that emerges 
without the sag point having developed, and yet. 
shows the presence of crystals along the section that 
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was in the hotter region of the furnace. Such a 
composition would present the least favorable condi- 
tion for glass formation. The opposite, and most 
favorable condition for glass formation, would be a 
specimen that emerges with a well-defined sag point 
and shows no evidence of crystallization anywhere 
along the fiber. 

The fact that crystals do not form in the tempera- 
ture region above the sag point does not necessarily 
indicate that the temperature of maximum crystalli- 
zation is higher than the hottest portion of the speci- 
men. It may indicate that the rate of crystallization 
is so slow that erystals do not form in the time that 
the fiber was in the furnace. This shows that the 
melt can be held for considerable time in this dan- 
gerous region with little fear of devitrification. 

In general, the difficulty of glass formation in- 
creases as crystals appear closer to the sag point. If 
crystals form on the cold side of the sag point, it 
indicates that even if glasses are formed, the problem 
of crystallization during the process of annealing 
would still remain. 

If a glass crystallizes in the annealing temperature 
range, the difficulties involved are not necessarily 
insurmountable. Fused silica, for instance, erystal- 
lizes in the annealing region, the crystal growth 
being initiated at the surface. It would appear, 
however, that if the criterion suggested by this pro- 
cedure, and by what has been said before, is adopted, 
then fused silica is far from being the ideal glass it 
is frequently described to be. If a glass is regarded 
as the opposite of a crystal, then those glasses that 
devitrify least readily come closest to being “ideal” 
from a theoretical as well as a practical point of 
view. Thus, glasses that contain “network modifi- 
ers’”’ and “intermediates” in addition to the so-called 
basic silica network, but which crystallize fairly 
slowly in the region of maximum crystal growth, 
come closer to being “ideal” than does pure silica. 
On the other hand, pure B,O,; can be crystallized 
only with the greatest difficulty and is more nearly 
an ideal glass in spite of certain (usually) undesir- 
able properties (i. e., high hygroscopicity, low melt- 
ing point, high coefficient of expansion). 


4. Results 


Table 1 lists several optical glasses * comparing 
the temperatures of the deformation point, sag 
point, and softening point. The fiber diameter 
chosen for the sag point determination (0.50 to 
0.80 mm) is practically the same as that used for 
the softening point (0.55 to 0.75 mm) [14] so that 
no correction for this dimension is necessary in 
comparing these two points. The last two columns 
give successively the difference between the soften- 
ing point and deformation point, and the difference 
between the sag point and the deformation point. 
For the various glasses, these reference points vary, 
the differences between them vary, and the sag 
points lie between the deformation and the softening 


? The compositions as well as the meaning of the symbols designating these 
glasses may be found in reference [13] 


TABLE | Temperatures of reference points of several Optica ' 
glasse s 
’ 
Temperature of 
CGilass type * '—a C~—a 
Deforma Softening 
tion Sag point peint! 
point * 
a (h 
( Cc ( 
BSC Sll 80) HSA Tol 7 12) 
BSC 517 604 653 720 14 125 
LC 423 suo 629 702 y 112 
BaC 441 640 650 737 Ww “7 
BaC A725 tit Hu7 ris) s xy 
BaC 574 H70 718 702 ts s? 
BaC 617 Hu2 731 7R5 at} 8 
CF 529 195 ‘ti 41 14 
F 620 Yr 515 aed s 1S 
* Compositions are given in [13], Leo Shartsis and Sam Spinner, Viscosity an 
density of molten optical glasses, J. Research N BS 46, 176-194 (1951) R P2199 
b ASTM Standards, Standard definitions of terms relating to glass and glass 
products, C 162-52 (1955 
e ASTM Standards, Tentative method of test for softening point of glass 


C333-MT (1955 


points. This is an example of the fact, mentioned 
earlier, that the temperature-viscosity relations of 
glasses of different compositions can vary widely, 
and consequently of the necessity for using several 


points to specify more fully the change in properties | 


of different glasses as they are heated. The par- 
ticular reference point chosen is governed by its 
ease of determination and proximity to the temper- 
ature range of the particular process involved. 

The sag point is not offered to replace either the 
detormation or the softening point. 
tion point is obtained from the determination of the 
thermal expansion by the interferometric method 
and is part of the data that are essential to establish- 
ing the correct temperature to use for the fine an- 
nealing of optical glasses. The sag point is deter- 
mined in a matter of minutes as compared with a day 
or more for the deformation point, and can be used 
to establish a maximum temperature for the com- 
parable annealing of experimental glasses. The 
temperature of the sag point is nearer to the anneal- 
ing range than is the softening point and _ hence 
devitrification troubles encountered in the latter 
determination are avoided. 


5. Summary 


(1) A new reference point, the sag point, in the 
heating (or cooling) of glasses has been defined and 
the method for experimentally determining this 
point described. The determination is rapid and 
simple. 

(2) The sag point has proved itself to be of con- 
siderable value in determming the maximum ab- 
nealing temperature of new glasses of widely varying 
compositions. 

(3) A method is also shown for obtaining addi- 
tional information concerning the devitrification 
tendencies of experimental melts and the feasibility 
of cooling them as homogeneous glasses. 

(4) Data have been presented in which the sag 
point is compared to previously defined reference 
points in the heating of glasses. 
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